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The interpretation of spectral band shapes in terms 
of the reorientational motion of the absorbing species 
has proven to be a powerful method for studying inter- 
molecular interactions and reorientational motions of 
molecules in the liquid and the gaseous state. In this 
thesis the infrared band shapes of gaseous and liquid 
mixtures of methane and the ESR linewidths of vanadyl com- 
plexes are interpreted in terms of the reorientational 
motion of the absorbing molecules. 

In Chapter II an extended rotational diffusion in- 
terpretation of Gordon's semiclassical band shape theory 
is developed and this modified theory is used to inter- 
pret the infrared band shape due to the V3 mode of 
methane. The angular momentum correlation times, Tye 
obtained from this analysis were in very close agreement 
with those obtained from an analysis using the classical 
model. This result indicates that the commonly used 
procedure of studying reorientational motion in terms 
of the Fourier cosine transform of the band shape is 
justified. This theory proves to be of limited applic- 
ability due to the details of the calculations. 

In Chapter III a "true" semiclassical extended 
rotational diffusion theory is developed. The memory 
function approach is used tosimplify the computation 
of reorientational correlation functions and spectral 


densities. This semiclassical theory is used to inter- 
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pret the infrared band shapes of methane due to V3 and 


v The angular momentum correlation times obtained 


Ae 
from the two different infrared bands and those obtained 
from an independent nuclear magnetic resonance T) study 
on closely related systems show a similar density de- 


pendence indicating that the values of T_ assigned by 


J 
the extended rotational diffusion theories have physical 


Significance and that T_ is not a meaningless parameter. 


J 
The ESR spectra of a series of vanadyl complexes 
in various solvents were obtained. These spectra are in- 
terpreted in terms of the Kivelson linewidth theory. A 
procedure for fitting the M-dependent part of the line- 
width using the Kivelson theory is introduced. This 
procedure is more representative of the Kivelson theory 
than the previous procedure of considering the mM? con- 
tribution to the linewidth because each of the parameters 
of the Kivelson model are considered according to their 
contribution to the linewidth. It is demonstrated that 
the Kivelson theory gives better agreement with experi- 
ment if each spectrum is studied separately rather than 
if all the spectra for a given complex/solvent system 
are studied simultaneously and the Stokes-Einstein re- 
lation is used to relate the reorientational correlation 
time, Toys to the temperature and viscosity of the solvent. 
This result indicates that possibly the Stokes-Einstein 


relation is the cause of the disagreement between pre- 


dicted and experimental results. 
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CHAPTER I 

The conventional spectroscopic method of studying 
a system is to assign the observed spectral frequencies 
to transitions between the various quantum states of the 
system. Because of its emphasis on energy levels, this 
approach might be referred to as the Schrodinger : pic- 
ture of spectroscopy. The Schrodinger picture, with its 
emphasis on transition frequencies, yields information 
about the environment of the absorbing species. This 
information includes such things as bond lengths, bond 
strengths, electronic environment and chemical inter- 
actions with solvent species and can in general be des- 
cribed as dealing with the time-independent or time- 
averaged characteristics of the systems. If a time- 
dependent parameter of a system causes the absorption 
frequency of a given molecule to vary with a frequency 
We, between two values separated by Aw, then by time- 
independent we mean ae << Aw, while by time averaged 
we mean i >> Aw. When ue << Aw the two lines will be 
resolved and the two values of the parameter can be 
calculated but when w >> Aw only a single line is 


observed and an average value of the parameter can be 


Ie 


determined. In the region when a, Aw the spectral 


line will be very broad and little can be said about 


the parameter. In most liquids and dense gases, the 


sctugunl to heeding yatidionessiiate! ac ER 
eninanoger® iste “bible i. 5 28/09 ee 
sits 3¢ eadnye eS aneezey isty! aaawe wrsd ; 
akitst pat wage 1G etenilan 234 1 
et * tom bdadbe: bent iad oo + boxset x ‘ath 


Jeot> bees 
ak a ia pire 


Seth ; 


-#ob sl 2 kewonae sean te ai we 


| (aanwpag? % ahs we oF aleaaton dehy 5 perpen 
Recuad nt poe 108 Ve avenacpe nbulay oi? scented | 4 

: oats ve! etide nh > q* Feet aw. sreheade beh 

od In Roath omy add. wh > eine ria a apeul- 
oye ou! eae 2sdenureem lt jo eautay' cvs ait bas bovkoges 
"gi wR okenie « yLao.wh: << at Feit ‘ed hesahwo as 
sd ree redann aay, itd Io oetie distur ns Fite Bowater 


Schrodinger picture does not present the most effective 
description of the observed spectra because the individual 
spectral lines are blended together to form a continuous 
band. This broadening and coalescence of the spectral 
lines is due to time-dependent interactions within the 
system. 1’? 
The Heisenberg formulation of quantum mechanics 
allows us to rephrase the Schrodinger expressions for 
spectral absorption so as to emphasize the time develop- 


ee The Heisenberg formulation of 


ment of the system. 
quantum mechanics places the time dependence of the sys- 
tem into the operators and the wave functions are time 
independent. The Schrodinger formulation on the other 
hand uses time dependent wave functions and time inde- 
pendent operators. The Heisenberg picture of spectros- 
copy with its emphasis on the time development of the 
quantum states has proved to be a powerful tool in the 
interpretation of spectral band shapes in which the in- 
dividual transition lines are not resolved.” 1° In 
addition the Heisenberg approach to spectroscopy pro- 
vides an easily visualized picture of how time dependent 
interactions influence observed band shapes. + af 9 'ate 

In the frequency region above 1000 cm +, where the 
influence of stimulated emission is negligible, the spect- 
tal band shape in the Schrodinger picture of spectros- 


pale 
copy is expressed by o 
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oA 
where the sum over i and f includes all relevant quantum 
states of the system, i indicating the initial and f the 
final state involved inthe transition. Ps is the Boltzmann 
factor for the quantum state i> and E, is its energy; 6 
is the Dirac delta function and w is the frequency of the 
applied radiation. The operator, Y, represents the inter- 
action being observed. For example, in infrared absorp- 
tion spectroscopy, Y is €.m” which represents the inter- 
action of the electric field of the incident radiation, 
é, with the vibrational transition dipole moment of a 
molecule, m’. The Heisenberg equivalent of Eq. (I-1) can 


be obtained by replacing the Dirac delta function by the 


equivalent Fourier integral, 


6(X) = — dte P (I-=2) 


where t is the Fourier variable and ~ = vy-l. Making this 


substitution, we obtain 


Figo ois o<il¥lO<ele]D f at exp[i((E, - E,)/M-w)t] 
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The quantity, e fi Y e f Thy, is the Heisenberg time de- 


pendent operator which we will denote as Y(t). Then Eq. 


(I-4) becomes 


2 
al 2,6 


I(w) = Hh fac ett) 5 <i] ¥ (0) [£><E]¥() [4d oe (t-5) 


Eq. (I-5) expreses the spectral density I() in the form 
of the Fourier transform of a correlation function G(t) 


defined by 


G(t) =>) p <ily (0) |£><E[¥(t) |4) >2, Pee (OK ist (t) 
if Let 


(356) 


where the * indicates the complex conjugate and Y,¢(t) 
is the matrix element (il¥(t) |£). 

The function, G(t), is the autocorrelation func- 
tion, for, the, quantity, Y,, represented by the operator Y. 
Let us consider the meaning of the correlation function. 
First of all we. must note that the correlation function 
"gives a description of the ensemble as a whole and not of 


an individual member of the ensemble. From Eq. (1-6) 
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we see that the correlation function contains information 
about the ensemble average of the individual products of 
ee tt) and Y5_¢ (0). Let us consider a hypothetical ex- 
periment in which all of the quantities Yi are the same 
for every molecule at time zero. The correlation function 
would then simply describe the result of the dephasing 

as time progressed. The dephasing is due to the differ- 
ing time dependence of the matrix elements, Yet), from 
molecule to molecule in the ensemble. In magnetic reson- 
ance spectroscopy, such an experiment can be performed. 

A pulse of incident radiation is used to give all tran- 
sition magnetic dipoles of the ensemble the same phase at 
time zero. The magnetization vector, whch is the sum of 
the transition magnetic dipoles, is detected. The magnet- 
ization decays in time and the decay curve is the cor- 
relation function for the transition magnetic dipole. 

The only purpose of the pulse in this experiment is to 
give all the transition magnetic dipoles the same phase 
at time zero (the time of the pulse) so as to produce 

a physically observable quantity. Discussion of cor- 
relation functions is still meaningful if the individual 
matrix elements, Yager are not in phase because the 
product considered does not depend on a phase which is 


universal to the system but depends only on the change 
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in phase of Yi¢ relative to its phase at time zero. 

The correlation function is the vehicle through 
which we can introduce the influence of molecular inter- 
actions to the band shape. In all systems studied in 
this thesis, the time dependence of the operator, Y(t), 
is due to the reorientational motion of the absorbing 
molecules. The correlation functions for Y in these 
systems will be referred to as reorientational correla- 
tion functions. 

A reorientation correlation function is character- 
ized by a reorientational correlation time, Tos which 


is defined as 


Cae ne} fat ro) : (I-7) 


O 
Tp is a measure of the time over which the molecule 
"remembers" its orientation in space i.e. the time re- 
quired for the molecule to reorient by ~1l radian. A 
second correlation time of importance in describing 
reorientational motion is the angular momentum correla- 


tion time, which may be interpreted as the average 


Ty: 
time over which a molecule remains in the same rotation- 
al state. The two correlation times are related, but 
the exact details of this relationship depend on the 


nature of the reorientation process. 


Reorientational motion affects various forms of 
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spectroscopy in different ways. In magnetic resonance 
spectroscopy the energy levels of a spin system depend 
on the orientation of the molecule relative to the 
applied magnetic field. Reorientation of the molecule 
thus modulates these energy levels and the resulting 
fluctuating energy levels give rise to linewidth contri- 
butions. In vibrational-rotational spectroscopy re- 
orientation of the molecule changes only the magnitude 
of the interaction between the incident radiation and 
the molecule, but does not alter the transition fre- 
quency. Rapid reorientation of the molecule causes the 
reorientational correlation function for the dipole- 
electric field interaction to decay rapidly. In terms 
of the band shape this means that the individual lines 
will be broadened. 

In order to further explore the influence of re- 
orientational motion on the band shape it is necessary 
to have a model to describe the reorientation of a 
molecule. There are two basic views on molecular re- 
orientation in liquids. One view suggests that a 
liquid resembles a loosely structured solid in which 
the individual molecules are restricted from moving 
freely because of interactions between species. In 
this theory reorientation comes about because of very 


strong torques exerted on a molecule when it is in- 
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volved in a collision. This approach has been developed 


b2-15 


by O'Reilly but it has not been widely exploited. 


Similar approaches have been suggested by Brown et een’ , 


and by Kivelson et aap! 18; 


as possible solutions to 
problems dealing with hydrogen bonding solvents. The 
second view pictures a liquid to be very similar to a 
dense gas. Reorientation models based on this picture 
assume that reorientation takes place by rotational 
steps inthe periods between collisions. These models 
are known as rotational diffusion models. The word col- 
lision here and in subsequent usage will be interpreted 
as an event in which the angular momentum of a molecule 
is altered. The models considered in this thesis are 
based on the dense gas picture for liquids. 

The first rotational diffusion model was developed 
by Debye ee to interpret dielectric relaxation experi- 
ments. In his model Debye assumed that reorientation 
took place by a series of short (<<l radian) rotational 
steps. The requirement of short rotational steps has 
limited the application of the Debye theory to polar 
liquias. 7° Another reorientation model based on 
rotational steps is the perturbed free rotor model. 

In this model reorientation is assumed to take place by 
long (>>1 radian) rotational steps. Again the length 
of the rotational step has restricted the application 


of this model and it has been used only in studies of 
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low pressure gases. 7? 

Recently Gordon 2 has developed a theory of molecu- 
lar reorientation which is not restricted by the length 
of the rotational steps. In this approach, which is 
referred to as the classical extended rotational diffusion 
theory, the classical motion of the molecules is followed 
in detail in the individual rotational diffusion steps 
and the reorientational correlation function is obtained 
by taking an ensemble average over these individual 
steps. Gordon's theory which was limited to linear 
molecules has been extended by McClung to the more compli- 
cated systems of spherical < and symmetric : top mole- 
cules. The classical extended rotational diffusion 
theories have been successfully used to interpret a 


[eee As 26-20 
‘ 


considerable number of infrared Raman 


and nmr 78-33 pana shapes. 

A limitation of the classical extended rotational 
diffusion theory is that its application is limited to 
those systems which can be described by classical 
mechanics and thus the use of this theory to interpret 
asymmetric band shapes is often questionable or impos- 


2,34,35 proposed a semiclassical 


sible. Gordon et al, 
theory to predict reorientational correlation functions. 
This theory makes use of semiclassical scattering theory 


to describe collisions and has been used to interpret 
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the infrared spectra of CO/He mixtures = and HCl/Ar 


mixcuLces. 


interpretation of Gordon's semiclassical theory. 


37 


In Chapter II we develop an extended diffusion 
2,34,35 


We have opted to pursue this approach rather than to use 


semiclassical scattering theory for several reasons. 


(1) 


(2) 


(3) 


(4) 


The Gordon approach using semiclassical Bea veri ng 
theory requires a good deal of preliminary 
mathematics in order to calculate the collision 
cross sections. These calculations are costly 

and in this sense make the theory of limited 


use. 


Semiclassical scattering theory requires the 
knowledge of an intermolecular potential. This 
in effect introduces another variable to the 


model. 


The semiclassical scattering calculations have 


been seriously questioned in the literature. 7° 


The success of the classical extended rotational 


diffusion theory Ue: suggests this picture 


of a collision to be a reasonable approximation. 


The band shapes predicted by this extended diffusion 


interpretation of Gordon's semiclassical theory are 
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compared with the infrared band shape of CH, due to v, in 


4 3 
the pure liquid ao in solutions of CH, in liguid inert 
gases 38 and in gaseous mixtures of CH,/He and CH,/N, at 


densities between 100 and 1000 amagats.>? 


Application of this semiclassical theory was re- 
stricted to band shapes which could be described by a 
small number of transition frequencies (<40). For ene 
reason it was not possible to study the infrared band of 
CH, due to V4 because of the large Coriolis coupling con- 
stant, C4,’ which splits each transition line into several 
components. 7° 

In Chapter III we developa "true" semiclassical 
version of the classical extended rotational diffusion 
model. In this model the number of transitions is not a 
limiting factor. The semiclassical extended rotational 
diffusion model is then used to interpret the infrared 
band shapes of CH, due to Vy in CH, /He and CH,/N. gas 
mixtures at densities between 100 and 800 amagats. >” Ln 
addition, the gaseous systems studied in Chapter II are 
reanalyzed. 

As pointed out above, reorientational motion of 
molecules in liquids can give rise to linewidth contri- 
butions in magnetic resonance spectra because the 
“rotational modulation of anisotropic magnetic inter- 


actions gives rise to relaxation of the spin system. In 


Chapter IV, we describe the results of an experimental 
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study of the ESR linewidths of several vanady1l-f-diketo- 


nate complexes inavariety of solvents. Kivelson et al, 


a have developed the theoretical relation between line- 


widths and the correlation times for molecular reorienta- 


tion. Spin-rotational coupling also makes a contribu- 


42,43 


tion to T, and T Kivelson et al, have developed 


2° 
theoretical expressions for this contribution. Recent 
investigations at low temperatures < and in hydrogen 


17,44,45 suggest that the expression 


bonding solvents 
for the spin-rotational contribution to Tos and possibly, 
the entire Kivelson theory are not adequate. The study 
of vanadyl complexes described in Chapter IV was designed 


to further test the Kivelson theories in non-hydrogen 


bonding solvents. 


-ssmiaieaiadie soem 20 aia: Womb wy 
he 28 wonkavies _agnevitoe erect’ 
| wenn! eesenrtend aoitehes sates ald 
et ae os aia neal vot Reel ne | 


hem —! ‘gtieal: 


“yhuse nt sua eee feth acts 


Par =) ) = i) 7 
bey) / | 
; ‘ 
) : a 
: 7 
3 : in + 
& i 
oie 
am Y I 
” 7 
a 
‘~ ee 
‘_ ' 
a 


SEMICLASSICAL BAND SHAPE THEORY 


CHAPTER II 
A. INTRODUCTION 


The reorientational motion of molecules can be 
studied by analysis of the infrared and Raman band 


2,4-6,8,9,46 In order to extract this informa- 


shapes. 
tion it is often necessary to make reference to a model 
for the molecular motion in the fluid. A number of 


5,6,8,9,18,47-49 
a 


theories based on the hypothesis that 


a liquid resembles a dense gas, have proven quite suc- 


cessful in accounting for experimental results.>7/110, 


aa. These models are generally referred to as rota- 
tional diffusion models. Rotational diffusion models 
consider reorientational motion to consist of a series 
of rotational steps which are terminated by "collisions". 
The reorientation during the rotational step is des- 
cribed by either free rotation of the molecule or by 
rotation executed under the influence of a retarding 
torque. The latter description has been employed by 
Debye 4 in his theory of dielectric relaxation. In 

the Debye model, the molecular reorientation is assumed 


OFA 


to follow a rotational diffusion equation 
rotational diffusion coefficient in this equation is 
related to the bulk shear viscosity of the liquid by a 
Stokes-Einstein hydrodynamic relationship.?> The 


Debye theory is limited to systems in which the diffu- 
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Sive steps are very short (<<l radian). This is a par- 
ticularly restrictive condition for systems of spherical 
molecules where the experimental results suggest that 
rotational steps areoften of the order of 1 radian.>* 

A perturbed free rotor model has been suggested for 
low pressure gaseous systems. 7? This approach, which 
assumes that the free rotation steps are long (>>1 radian), 
is a lifetime broadening interpretation of the band 
shape. This approximation neglects all but the first 
collision 8 and the theory has proved of limited applic- 
ability even in low pressure gases. 7° 

Recently Gordon ? has suggested a rotational dif- 
fusion model which is not restricted by the length of 
the rotational stép. In this approach, which,is, referred 
to as the classical extended rotational diffusion theory, 
the classical motion of the molecule is followed in 
detail in the individual rotational steps and the re- 
orientational correlation function is obtained by taking 
an ensemble average over all possible sequences of 
rotationasteps: ijinmsthisimodelyvaricollisiontassyiewed 
as an instantaneous event in which the angular momentum 
of .a molecule is changed. Two limiting cases are con- 
sidered: M-diffusion and J-diffusion. In the J-diffu- 
sion model, the magnitude of the angular momentum 


vector is randomized onto a Boltzmann distribution and 


the orientation of the angular momentum vector is 
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completely randomized at every collision. Each colli- 
sion randomizes only the orientation of the angular 
momentum vector in the M-diffusion limit, and the magni- 
tude of the angular momentum does not change. A more 
sophisticated approach a7 permits one to describe the 
reorientation by any desired combination of M- and J- 
diffusion processes. Gordon's theory which was limited 
to linear molecules has been extended to the more com- 
plicated systems of spherical e and symmetric har? top 
molecules. The classical extended rotational diffusion 
theories have been successfully used to interpret a 


geo 2 ana Baek, 26228 


28-33 
band shapes and nuclear relaxation times. 


considerable number of infrared 


A disadvantage of the classical extended rotation- 
al diffuston theory is that *it is limited to ‘those 
systems whose rotational motion can be described by 
classical mechanics. The classical approach predicts 
symmetrical band contours and the application of clas- 
sical extended rotational diffusion theory is often 
questionable or impossible in systems with asymmetric 
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band shapes. Gordon has proposed a semi- 


classical theory to predict reorientational correlation 
functions. This theory, which in many ways parallels 
the modified Bloch equation treatment of chemical 
exchange used in magnetic resonance a2y makes use of 


: 34,36 : 
semiclassical scattering theory : to describe 
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collisions. Semiclassical scattering theory is a pro- 
cedure which approximates quantum mechanical scattering 
theory by using classical scattering theory and the cor- 
respondence principle. In this procedure a particular 
intermolecular potential is adopted and classical scat- 
tering theory is used to predict the trajectories of 
particles before and after they have been involved in a 
collision. The classical angular momentum of the particle 
before and after the collision is determined and the 
semiclassical scattering cross section is calculated by 
rounding off the classical angular momenta to the near- 


86737 have used this 


est multiples of fi. Gordon et al, 
procedure to study reorientation in CO/He and HC1/Ar mix- 
tures at low pressures, 

In this chapter, we will develop an extended ro- 
tational diffusion interpretation of Gordon's semi- 
classical theory. In section II-B this semiclassical 
band shape theory will be developed, and in section 
II-C the theory will be applied to a study of the V3 
infrared band of CH, in liquid and gaseous mixtures.°°" 39 
Section II-D will be devoted to a comparison of these 
results with those obtained by application of the 
classical theory 3 and with the results of an independ- 


56 


ent NMR-T, study on CH, in the dense gas phase. In 


z 
addition, the limitations of the theory will be dis- 


cussed. 
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B. DEVELOPMENT OF SEMICLASSICAL BAND SHAPE THEORY 
se The Gordon Semiclassical Theory 


The Fourier transform of the infrared absorption 
band shape is the dipole reorientational correlation 


function defined by a 


ae 4 ; | 
G(t) = > o.<ilm%o)]/£><E|m’ (t) [i> = > o,d, (0)d, (t) 
a! kk k 
Bae k 

(II-1) 
where Os is the Boltzmann factor for the rotational 
state |i) and the sum over i and f is over the initial 
and final states respectively. The subscript k used 
here and in the following equations replaces the pair 
of subscripts i,f. The Heisenberg operator, m’ (t), rs 


the dipole operator at time t, and d, (t) is 


dj (t) = (£[m’(t) |i). ar) 


Equation (II-1) follows from Eq. (I-6) by replacing Y 
by the dipole operator, mY’. In order to calculate thi 
correlation function, it is necessary to know the time 
dependence of dy (t). If the molecule is not involved 
in collisions which cause changes in its angular mom- 
“ entum (and thus its reorientational motion) dy (t) is 
readily calculated. Replacing the Heisenberg operator 


m (t) by the equivalent expression a 
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tHt iy) -itt/fA 
e m’ (O)e , Eq. (II-2) becomes 
Sy : tw t tw t 
di ft) = <fim" (0) |i) .e = 4 (0)e Caren 


where Oy, is defined in terms of the rotational energies 


Ey and Er by 


W. = (Ep - E,)/ * (II-4) 


In the absence of collisions the transition dipole 
4u.-t 
moment dy (t) oscillates as e with amplitude d,, (0) 


and its time-derivative is given by 


d[d, (t)] 
————— = tw,d,(t) . (EL-5) 
dt 


Eqs. (II-3) and (II-5) are exact provided that the 
molecule is not involved in collisions. If the molecule 
is involved in collisions then the reorientational 
motion of the molecule in the other angular momentum 
states must also be included in the calculation of d, (t). 
One approach to this problem would be to follow the re- 
orientational motion in each of the rotational states 
which the molecule samples between time 0 and t then 

dy (t) could be calculated as the statistical average 

of all possible sequences of sampled states. This is 


the procedure employed in extended rotational diffusion 
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theory and will be considered in Chapter III. Gordon 
has suggested an alternative procedure which is similar 
to the "jumping spin" treatment of chemical exchange in 
magnetic resonance spectroscopy.” Gordon suggests that, 
instead of calculating the product in Eq. (II-1) for 
each individual transition dipole moment, the product 


be calculated at time 't for the sum of the transition 


dipole moments ofa particular form. He replaces Eq. 


(II-1) by 
a 
G(t) = dy d, (t) ave (II-6) 
k 
where dir the intrinsic line amplitude, is the quantity 
Saye : 
(£ |m | 2) and d,(t),,,. is the sum over all molecules, 


at time t, of the transition dipolemoments of the form 
([m” |3) - The subscript, Ave., has been introduced 
to emphasize the "averaged character" of q(t) ave, and 
to distinguish it from dy (t) which refers to the transi- 
tion dipole moment on a particular molecule. In order 
to clarify the wording in the development to follow, 
Qe) eee will be referred to as the line amplitude of 
the k-th spectral line. 


; ; P F £ 
Since, q(t) ave, is simply the sum at time t o 


We ; 
matrix elements of the form (é |m |i) the time de- 


pendence of a can readily be described simply 


by adding transfer terms to the equations of motion 


Let 
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[Eq. (II-5)]. These transfer terms would account for 
the collisional transfer of line amplitude from the k-th 
spectral line to,all other lines and for the transfer of 
line amplitude to thek-th spectral line from all other 
lines. The time dependence of d, (t) can then be 


Ave. 


expressed by, 


d ime 
pach rvs. ST nee oe ae 
i SD el CEN es (II-7) 
nea 


where the first term on the right represents the change 
in dt) ave. due to the free motion, while the second 
and third terms represent the transfer of line ampli- 
tude out of thek-th line and the transfer of line ampli- 
tude into thek-th line respectively. Thus yy gives 

the rate at which a unit amplitude in line k would de- 
crease due to collisions, while ~The gives the rate at 
which collisions would transfer unit amplitude in line l 
into line k. The set of simultaneous differential equa- 


tions represented by Eq. (II-7) can be expressed in the 


matrix form 


a pa . a 
a Aves iat (70, LONE) ae ae 
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where A(t) ro, is a column matrix whose k-thelement is 


d, (t) and Wo is a diagonal matrix with kk-th element 


Ave. 
W,-+ The matrix [] is defined above. 

Eq. (I1I-8) can be solved subject to the boundary 
condition that we start with an equilibrium distribution 
Ofyiline amplitudes, ji.e:; d, (0) nue. = a, Pie where Py is 
the equilibrium population factor for the k-th line and 
dy is the intrinsic line amplitude defined earlier. 
SolvingEq<| 411-8) ) subject to this initial condition, we 
obtain 


a(t) (exp [(¢w, -11)t)).P.4 (II-9) 


Ave. 
where P is a diagonal matrix with kk=th element PL and d 
is a column matrix with k-th element d,- 


The reorientational correlation function defined 


in Eq. (II-6) can now be expressed by 


Il 


+ 
asd) 


G(t) Spt 
k + A 
d -exp [(tw, -f])t].-P. 


2 Qu 


where at is the transpose of the matrix d. The Ccorres= 


ponding expression for the spectral density is obtained 


by the Fourier transformation of Eq. (II-10): 
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Toh) = Im}a" . (wl rie ae )-t.p.a} (II-11) 
where w is the frequency of the incident radiation and 1 
is a unit matrix; Im{} implies the imaginary part of 
the quantity within the braces and ( ie implies the 


inverse of the matrix within the brackets. 


ae Extended Diffusion Interpretation of the Gordon 
Semiclassical Theory 


3655 7 


Gordon et al. have determined the matrix, I], 


for the CO/He and HC1/Ar systems using semiclassical 
scattering theory. Because semiclassical scattering 


theory calculations are very costly and possibly un- 


36 


reliable , and because they require a detailed know- 


ledge of the intermolecular potentials in each particu- 


lar system, we have opted to interpret [] in terms of 
the classical extended rotational diffusion mode1l.>'° 
This might appear to be an arbitrary decision, but we 


feel that it is justifiable in terms of the success 


achieved in the application of the classical extended 


rotational diffusion models. /*227° 


In the classical rotational diffusion models 576, 


See in which the rotations of the molecules are 


treated using classical mechanics, collisional events 


are assumed to occur at an average rate Ts , where Ts 
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is the mean time between events. Since each collisional 
event destroys all correlation in the angular momentum 
of the molecule,the angular momentum correlation time 
is TS: The "collisions" in the extended diffusion models 
are assumed to randomize the direction of the rotational 
angular momentum, and may randomize its magnitude as 
well. Analogous semiclassical rotational diffusion 
models may be constructed so that the discreteness of 
the rotational levels may be taken into account, but 
the basic intuitive simplicity of the diffusion picture 
is preserved. 

The rotational states of a spherical or symmetric 
top molecule are described by Lanes quantum numbers: 
the total angular momentum quantum number J and the 
projection quantum numbers on molecule-fixed and space- 
fixed directions, K and M respectively. The semi- 
classical analogs of the classical rotational diffusion 
models are models in which collisional events randomize 
the projection quantum numbers K and M in both M- and 
J-diffusion limits, and randomize the quantum number J 
onto a Boltzmann distribution in the limit of J-dif- 
fusion. Since the rotational states of a spherical 
top molecule within a J-manifold are degenerate, energy 
“transfer can occur at the collisional events in J- 


diffusion only. In the semiclassical J-diffusion 


model, the probability that a molecule initially in 
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the state |5,K,M) is put into the state ary : aaa m') 
by a collisional event depends only on the Boltzmann 
factor for the final state. In the semiclassical M- 
diffusion model, the quantum number J is not changed 
by collision, but the quantum numbers M and K are 
randomized. 

The transfer matrix I] does not describe the rate 
of transfer of molecules between rotational states, it 
describes the rate at which intensity amplitude (tran- 
sition dipole moment) is transferred from one spectral 
line to another. One would expect, however, that the 
rates of these processes are related. A molecule in 
the rotation state |5,K,M) of the ground vibrational 
state can absorb infrared radiation at as many as three 
distinct ~cequencies, Le. it may give rise to absorp— 
tion lines in each of the P, Q and R branches of the 
spectrum. Consequently, the determination of the 
rate of transfer of line amplitude from one spectral 
line to another is somewhat more complicated than the 
determination of the rate of transfer of systems from 
one rotational level to another. In order to relate 
the rates of transfer of line amplitude between 
spectral lines to the rate of transfer of molecules 
between rotational levels, we make the argument that, 
although the individual molecules in the J-th rotational 
level are associated with more than one transition 


moment (i.e. can absorb at more than one frequency), 
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we can divide the molecules in the J-th rotational level 
into groups associated with each of the possible transi- 
tions. This assumption implies that if one molecule is 
transferred out of the J-th rotational level, the P-, 
Q- and R-branch lines, which are derived from this 

(3d 


INP Ie J(Q) 
and ds (p)/3 respectively, where ds is the intrinsic line 


rotational level,lose intensity amplitudes d 73 
amplitude of the i-th line in the spectrum. The inten- 
sity amplitude which is lost by the i-th line when one 
molecule is transferred out of that line is dj, and the 
intensity amplitude gained by the j-th line to which 
transfer occurs is ore In general, d; varies from line 

to line in the vibration-rotation spectrum. Since the 
fraction of collisional events which transfer a molecule 
from the J-th rotational level to the J'-th level, in 

the J-diffusion limit, is given by the Boltzmann factor, 


re) for the J'-th level, the fraction which transfer 


gt 
line amplitude from the i-th spectral line, which is 
associated with the J-th level, to the i'-th spectral 
line, which is associated with the J'-th level, is P 5/3. 
By analogy to the classical rotational diffusion models, 
we take the average rate of collisional events to be 
Eee Then the transfer matrix, IY, which describes the 


~ rotational behaviour of molecules undergoing J-diffusion 


has elements 
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(1-P;)/t, = 
Ni; = (ers 13) 
Ean Solids te ifj 
where P. is equal to one-third of the Boltzmann factor 
for the rotational level with which the i-th transition 
is associated. 

The prescription for the elements of the transfer 
Matrix in the semiclassical M-diffusion model is some- 
what different because the collisional events are 
assumed to randomize the direction of the angular momen- 
tum vector, but not change its magnitude. [In this 
limit, line amplitude is transferred only between lines 
in the different branches of the spectrum which origin- 
ate in the same J-level. In the absence of collisions, 
the equations of motion for the line amplitudes of the 
transitions originating in a common J-state are not 
coupled, so there is no reason to expect line amplitude 
transfer between these lines unless collisions occur. 

We therefore assume that the transition matrix, ae in the 


M-diffusion model has elements 
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where Ty is the mean time between collisional events 
just as in the J-diffusion model. 

A further subdivision of the molecules in an en- 
semble of spherical molecules like methane, which 
contain several “sperns oan nuclei with nuclear spins, 
is required because the molecules exist in different 
nuclear spin isomers, analogous to the ortho and para 
isomers of hydrogen. 7° Since the rate of interconver- 
sion of the nuclear spin isomers is several orders of 
magnitude slower than the rate of rotational energy 
transfer in methane au one must view the infrared 
spectrum of methane as a superposition of subspectra 
arising from ortho-CH, (symmetry F), para-CH, (symmetry 


E) and meta-CH, (symmetry A). The band shape of methane, 


4 
therefore, should be calculated in three separate com- 
putations - one for each symmetry species - which allow 
transfer of line amplitude only between the lines of a 
Single nuclear spin isomer. 
Ce APPLICATION OF THE THEORY TO CH, BAND SHAPES 
ifs Calculations 

The vibration-rotation spectrum of the V3 mode of 


CH, in the low pressure gas phase is well resolved and 


4 
exhibits some splitting of the lines arising from the 


rotational levels with high J, by Coriolis interactions.° 
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The observed absorption frequencies of CH, have been 
Keportedwby:Plyer et atey and formulae for the 
intrinsic line amplitudes, d., have been given by 
Childs and Jann. 92 The Boltzmann factors for the 


rotational levels of the various nuclear spin isomers 


of CH, were computed using the weights given by Wilson 


4 
and Herzberg. 79" 61 


The calculation of I(w) and G(t) for the V3 band 


of CH, with a given [J] matrix was accomplished using 


the diagonalization algorithm suggested by Gordon 


et al. 22184 


This algorithm improves the efficiency 
of the calculation of I(w) because it takes advantage 
of the fact that the inverse of a diagonal matrix, C, 
is the diagonal matrix whose elements are the inverses 
of the elements of C. The G(t) calculation is simpli- 
fied by noting that the exponential of a diagonal 
matrix C is a diagonal matrix whose kk-th element is 
exp[C,,]. Using these algorithms Eq. (II-10) and 


Eq. (II-11) become 
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where Re{} implies the real part of the quantity within 
the braces and the complex matrices § and A are defined 
by 

=] ; 

S : (tw-IT) ~ <5. =k (II-16) 


where A is a diagonal matrix and Ses is the inverse of 
S. The calculations were performed on an IBM 360/67 
digital computer using complex matrix inversion and 
diagonalization procedures throughout. A listing of 
the Fortran program is given in Appendix II-A. It was 
found to be more efficient to calculate G(t) by numeri- 
cal Fourier transformation of I(w) calculated from 
Eq. (11-15) then to calculate both G(t) and I(w) inde- 
pendently from Eq. (11-14) and Eq. (LI-i5).. “This 
numerical Fourier transformation procedure was used 
throughout and employed trapezoidal rule integration. 
In a set of preliminary calculations the band 
shape for the F-species was calculated in two different 
ways. In the first calculation, all 44 spectral lines 
associated with the rotational states in the range 
0 < J< 7 were included. In the second approach, the 
transition frequencies associated with a given J-state 
in a given branch were averaged and the averaged fre- 
quency was used in the band shape calculation.- This 
method resulted in a band shape calculated from only 


21 lines. The nature of the averaging process used is 
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indicated in Fig. 1-:(a)-(b). The resulting spectra for a 
is 


Ty aye Choe ame ae sec, representative of a CH, /He mixture 
at 100 amagats, were identical. Since the Coriolis split- 
tings of the frequencies averaged in this latter calcu- 
lations are representative of those encountered in the 
range 0 < J < 12 and because the correlation time is 
representative of the least dense sample studied, it 

is apparent that neglect of Coriolis coupling is ee 
justified in these calculations. 

In a second set of calculations the overall 
methane band shape was calculated by superposition of 
the component ortho-; para- and meta-CH , subspectra cal- 
culated from frequencies averaged over the Coriolis 
components of each as described above. In order to 
assess the importance of the separation into subspectra, 
a series of band shapes were computed by relaxing the 
restriction that line amplitude be transferred only 
between lines of the same nuclear spin isomer. In these 
latter calculations, the Boltzmann factors for the 
rotational levels were taken to be the sums of the 
Boltzmann factors for the nuclear spin isomers, and the 
absorption frequencies of the lines for each nuclear 
spin isomer were used to construct an average absorp- 
tion frequency [see Fig. Te eee Methane band shapes 
‘calculated using this procedure are given in Fig. 2 


as the dashed curves. Also included in Fig. 2 are the 


corresponding band shapes for the individual nuclear 
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FIGURE 1 


Construction of average frequencies for the P(9) 


lines of CH, 


(a) All multiplet lines in P(9). 


(b) Average of Coriolis multiplets for each 


nuclear spin isomer. 


(c) Average of nuclear spin isomer frequencies. 
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FIGURE 2 


V3 band shapes of CH, computed from the semiclassical 


J-diffusion model 
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——— - Individual nuclear spin isomer spectra, 

----- - spectra calculated using line frequencies 
obtained by averaging all Coriolis multiplets of 

all nuclear spin isomers. The vertical bars under 
the computed band contours indicate the positions 

and relative intensities of the line frequencies © 
used in each calculation. The solid curves in the 
lowest section of the figure represent the normalized 
sum [(5A + 2E + 9F)/16] of the three nuclear spin 


isomer spectra. 
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spin isomers of methane and the band shape calculated by 
taking the weighted sum of these isomer band shapes. 

The band shapes in Fig. 2 (a) and (b) approximate the 
observed spectra of a CH, /He mixture at 100 amagats 

and the observed spectra of liquid methane respectively. 
It is clear that the separation of the band shape com- 
putation into separate calculations for each nuclear 
Spin isomer is unnecessary, in the range of Ts of 
interest, since the differences in the spectra calculated 
by the two methods are insignificant. We have based 
subsequent calculations on a total methane picture which 
allows collisional transfer of intensity between lines 
belonging to different nuclear spin isomers, since this 
reduces computation time. 

In Figs. 3 and 4, the band shapes calculated 
with the J- and M-diffusion limits of the semiclassical 
diffusion model are given for a range of values of Ts: 
It is clear that the two limits of the diffusion model 


give similar results when tT. is long, but give signifi- 


J 
cantly different band shapes in the limit of motional 


narrowing where Ts is short. Motional broadening of 


the individual vibration-rotation lines in the dilute 


gas limit, where t. is long, is more efficient in the 


J 


J-diffusion model, for a given value of Tye thanran 


the M-diffusion model. This observation is readily 
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FIGURE 3 


V3 band shapes of CH, computed with the semiclassical 


J-diffusion model at 295K 


The angular momentum correlation times, Tz, are 
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FIGURE 4 


V3 band shapes of CH 


4 computed with the semiclassical 


M-diffusion model at 295K 


The angular momentum correlation times, T 5, are 
(a) BN be < eos sec 
(b) 7.0) % 10), sec 
(c) 4.2 x yas sec 
(a) 2.0 x 10°77 sec 
(e) 1.0 x Vom sec 
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understood when one considers the lines between which 
intensity amplitude is transferred in the two models. 
In the J-diffusion model, line amplitude can be trans- 
ferred between lines within a branch (spacings of 9 cm”? 
LOmLE cay while the M-diffusion model allows only 
those transfers which are between lines associated with 
the different transitions from a single rotational level 


to 250 en 


(spacings of 8 cm. ). The amount of broaden- 
ing of the lines is related to the spacing between them 
as well as the magnitude of the off-diagonal elements of 
I] which connect them - the smaller the spacing and the 
larger the off-diagonal element, the broader the ine. ° 
In the limit of motional narrowing, the semiclassical M- 
diffusion model predicts a band shape with a sharp spike 
in the region of maximum absorption just as the classical 


6,8 The J-diffusion 


extended M-diffusion model does. 
models (both classical and semiclassical versions) give 


a less rapidly varying spectral density in this region. 


Z< Comparison with Experiment 


In comparing experimental results with those 
predicted by the semiclassical rotational diffusion 
models, one may choose to compare either the observed 
band contours with spectral densities computed from 


Eq. (II-15), or one may compare the reorientational cor- 
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relation functions obtained by Fourier transformation of 
the band contours with theoretical correlation functions. 
We have found that the correlation functions obtained in 
our calculations were unreliable because we truncated the 
rotation-vibration spectrumatthose lines arising from 
the J=12 rotational level. The neglect of the lines 
associated with J>12 was necessary to maintain camedes” 
tional efficiency by keeping the order of I] as small as 
possible. It did however lead to some unreliability in 
the spectral amplitudes in the wings of the spectral 
band and this unreliability is transmitted and amplified 
when the correlation function is obtained by Fourier 
transforming the computed band shape. We have chosen 
therefore to compare experimental and computed band 
shapes rather than their Fourier transforms. 

Examination of Eq. (II-10) explains this depend- 
ence of G(t), at small t, on the wings of the spectral 
density. The contributions to the correlation function 
are determined by two factors: a damping term exp [-[It] 
and a term exp[tw t] which oscillates with frequency 
Wo° At short times the oscillating terms are nearly in 
phase and the contributions to G(t) reinforce each 
other. As time progresses these terms get out of phase; 
’ theckarger «the tdifferencelim the sinequenezes: ofoscil- 


lation, the greater the dephasing. Therefore the con- 
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tributions from the transitions in the wings (large w) 
are dephased much faster than the contributions from the 
transitions near the center of the band (small w). Be- 
cause of this dephasing the contributions to G(t) from 
the wings of the spectral density soon cancel each other 
as time increases and their influence becomes negligible. 
The fitting of observed band shapes with dhexe 
calculated from the semiclassical rotational diffusion 
models involves comparing experimental spectra with a 
series of theoretical ones calculated with different 
values of Ty in both J- and M-diffusion models. In pre- 
paring the calculated spectra for comparison with 
observed ones, one must choose a procedure for normaliza- 
tion of the spectral amplitudes. One can require the 
calculated and observed spectra to match at the point 
of maximum absorption (intensity normalization) or one 
can require them to be normalized such that the areas 
under the calculated and observed curves are equal (area 
normalization). Intensity normalization tends to weight 
the central section of the band very strongly, while 
area normalization allows a more realistic weighting 
over the whole spectral band and is the preferred method. 
As we have indicated above, the wings of our calculated 
spectra are not too reliable because we mays not in- 


cluded all of the weak transitions which would contribute 


rAd Ah it i ‘i eon 
«58 mira Ss Aas 
: 


ta ae plat Fe eas fo 


a> mont sRob sida ae ah 
ae’. Go) Diege babe ss, to weirs sas Shen. 
mor. (h)'p dd ei Perpren ete abds. a8 
illite, iesiada dentin * Foon aiken pasta itt to 
atdipilt gen sepisgelad Perera hesnaticeal rapes 
sant? chi, aaa tial Rermado. fe en idas? aaol 
abe tae Sener s a om fedeantoinoe est ret 
bbw atone Le aca grlbtaqaoe. malo: : 
derma PLES agile ins 9 Lepshint were ippkvo roast? Foi 
—sKqy at @Seisbin rode, aoe he gE! cok . ‘+ bias 
29 iw ove agai Sol ere hod shwokge wale 9 
“ost isnton oa% Hebecotey i SRE Ce ee 
att shape sige. 60 | Seabat igen ferseqs ed3 30:4 
trkeg Sis) hk hoi gi af mapbeme peveoue bis 
HHS: RO {nozsswie BIO 

cee ald edd. owe, ophtanven od ad mots waluipor. ee i 
BOTH) - heupe bas) a ai bevsesicte Onb eae Eine acto a a | 
tian | shila, mised l tenor urkenednt 
lide. onkwned®: yer Baie “esti hid noideae bwidqan’ ot 
ninktctes, pitaiiees’ Sic x yeti Es Oks Eieamt nos 
lsouttieh bostotaty ont: ae oe Bived Cemboasté Stoo ot seve an 
podatuatan prt tdi Nrecie odd (OVDER “sotoanite aie ean | 
«hon ‘gv Bil sal peracid bhaaeos cs . ¥ We: om eudoaqe 
i aie aie 1a, ang Me aad tabeet 


43. 


to the amplitude in the wings. We have chosen to take 
the spectral area in a 120 em? region symmetrical about 
the frequency of maximum absorption as the area used in 
the normalization. The best fits of observed spectra 
with those calculated from the semiclassical diffusion 
models were judged by visual comparison. 

The infrared data for methane in high pressire 
mixtures with inert gases were obtained from Armstrong >”, 
and those for liquid CH, and CH, in liquid noble gases 
were obtained from Cabana. 78 Equally satisfactory fits 
of the observed high pressure gas spectra were obtained 
using the M- or J-diffusion model. Representative 
examples of the agreement between calculated and observed 
band shapes are given in Fig. 5. The values of the 
angular momentum correlation time which were obtained 
from the comparisons of calculated and observed spectra 
are given in Table 1. It should be noted that the 
spectra were calculated using frequencies appropriate 
to the dilute gas phase and agreement between calculated 
spectra and those observed in the dense gas can only be 
achieved if the whole of the calculated spectrum is 


- to higher frequencies for the 


shifted by 3-10 cm- 
gaseous systems and 5-13 om + to lower frequencies for 
the liquid systems. This is not too unexpected because 


we have only considered the broadening effects of 
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FIGURE 5 


Comparison of the observed V3 band shapes of CH, in dense 
gas mixtures at 295K with spectral band shapes calculated 


with the semiclassical J-diffusion model 


The band shapes are area normalized over a 120 cm + region 
about the center of the band. The intensity scale numbers 


are for comparison purposes only and are relevant only 


within ‘the figure. Observed :— , Calculated - ==-=== ° 
(a) CH,/He, 96 amagats; tT .5 MAE pte 1006p ise 
(b) CH,/N,, ° 186 amagats; 9) v7 =)4- 25% 10°73 sec 
(c) CH,/He, 543 amagats; Tae Zeon NOs sec 
(da) CH,/N,, 463 amagats; ci Pia x 10°? sec 
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molecular collisions. Frequency shifts can also be pro- 
duced by intermolecular interactions e2. but we have 
made no attempt to include them. Our neglect of these 
frequency shifts may account for the poorer quality of 
the fits of the highest pressure gas spectra. 


The spectra of CH, in the liquid phase and in 


4 
dilute solutions in liquid rare gases were Moneta with 
calculated spectra, and it was concluded that the best 
agreement was obtained with the J-diffusion model for 
the neat liquid, and with the M-diffusion model for the 
liquid rare gas solutions. This is in accord with the 
conclusions drawn in an earlier analysis of these spectra 
using the classical extended diffusion models. / More- 
over, the estimates of the angular momentum correlation 
times Ts which give best agreement between the observed 
spectra and those calculated from the semiclassical 
diffusion models are in close agreement with the correla- 
tion times required to obtain a satisfactory fit of the 
symmetrical part of the observed band shape with the 
spectral densities obtained from the classical extended 
diffusion models. 

The calculated band shapes discussed above have 
been compared to experimental spectra assuming that the 


experiments were performed under conditions in which 


instrumental broadening of the bands was insignificant. 
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Since the exact nature of the slit function of the 
Spectrometer was not known, the effects of convoluting 
theoretical absorbance spectra with a triangular slit 


Sanctions Gmeideh (455. cn 


were studied. The band shapes 
computed in this way were compared with observed spectra, 
and with the exception of the low density methane-helium 
band shape, the slit function had very little effect on 
the band shape, and had no effect on the assessment of 
the angular momentum correlation times obtained in these 
comparisons. For the spectrum of CH, in helium at 96 
amagats, the slit function had a significant effect on 
the shape of the band, and the value of Ts which yielded 


best agreement between theory and experiment was 12 x 10 


S$@C pein, contrast.to;thervalue.o£f.9.0_x aie sec obtained 


tS 


neglecting instrumental broadening. Because the effects 
of the slit function on the band shape are of the same 
order of magnitude as the effects of collisional broaden- 
ing in this low density gas mixture, the angular momentum 
correlation time obtained in our analysis is subject to 
large uncertainties. 

The analysis of the gas phase spectra did not give 
a clear cut distinction between M- and J-diffusion limits. 
The reason for this is the similarity between M- and J- 
diffusion: spectra in the range of interest. The semi- 


classical correlation times of the angular momentum from 


Bb 


ih 


mri 40 gos tor? ahh: ash ae a i 
caktedlovaoo 10 mins. eee g sete 
Mise  tely QUsie? A bi ie 

pre ied amet, feibo te: . 7.2 

eorrosaw 5: vase ley (htt casas le ‘etd mb 

wnt teat xedion, yeas wor sidan Birhioe sw 6 14 


7 + ei: "ate 
ms 308705 sass tings a8 oe votes te. wth, 


Kam tnstila waa” et 8 ak inh 
agaa tials ogee oda | 
emda Oud. 24 a i fais, 

; rrr 
mt hee we DBAS “oa ee eS the ik et 
ot Fitigiie 2i akaytoma: sos if banisddo. Lia oiaewn 1 
7 do ee patil lie sca pean 
ovip Jom ba as Sasry aap Sit Se), akeptans ok: ip 
vad bin soiadeaib-0 Bike = neawtad ct gel Si ls s 
a Bik A. gacwted Rai tal Lat Bt? at olty, sh anon sant - 
ee ade eaedre, 10) eae, hy ee Mf inh aobauat hb 
a aginst ety PONT aptuigae Act's te. ‘momid nobis anit aepetonnee 


49. 


this analysis are given in Table 1 for both the M- and 
J-diffusion limits. Except in the cases of the very 

high pressure samples we see that the values are essen- 
tially identical. In the high pressure samples M-dif- 
fusion exhibited motional narrowing effects which were 


inconsistent with the experimental band shapes. 


D. DISCUSSION 
bs Comparison of Semiclassical Results with Previous 
Treatments 


The magnitudes of the angular momentum correla- 
tion times obtained by fitting the observed spectral 
band shapes with band shapes calculated from the semi- 
classical rotational diffusion models can be compared 
with those obtained from the analysis of proton magnetic 
relaxation times which are dominated by spin-rotational 
interactions. Since the only relevant experimental 
data in the literature consists of proton relaxation 
rates for methane gas at densities up to 550 amagats 2c) 
the angular momentum correlation times obtained from 
the relaxation data are not strictly comparable to those 
obtained in this work because the conditions were not 
identical. However, it is of interest to ascertain the 
reasonability of our estimates of Ty: Trappeniers 
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orientation to analyse their data. As noted above, 
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this implies that the effects of only the first 
"collision" are included and the anisotropic spin-rota- 
tion contributions to the NMR relaxation times are 
treated in an approximate fashion. We have reanalysed 
their data using the classical extended rotational dif- 


fusion theory. ° The values of tT. obtained in this way 


J 
are shown in Fig. 6 together with the present Ts values 
for CH, in dense gas mixtures. Although, as mentioned 


4 


earlier, the numbers are not directly comparable becaus 
of the difference in the perturbing species, we see 
that the angular momentum correlation times obtained by 
the different methods are very similar. 

BineecHe general conclusions of the analysis of 


the v, band of methane both in liquid and high pressure 


3 
gas systems by the classical extended diffusion and 

semiclassical diffusion models are the same, it is of 
interest to compare the band shape computed from the 


Sa! with the sym- 


classical extended diffusion model 
metric portion of the band shape computed using the 
semiclassical rotational diffusion model. Several 
representative examples are shown in Fig. 7 for the 
J-diffusion limit. It is clear that, with the excep- 
tion of the rotational fine structure which appears 


' in the semiclassical spectra at long Tye the two 


models give similar results for correlation times 
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FIGURE 6 


Density dependence of tT. determined by comparison of 


J 


observed v., band shapes with those computed using the 


3 


semiclassical J-diffusion model. 


O - CH ,/He, A-CH,/N,, o-determined from proton 


relaxation times in methane gas (ref. 56). 


FIGURE 6 


Density dependence of TS determined by comparison of 
observed V3 band shapes with those computed using the 


semiclassical J-diffusion model. 


O - CH ,/He, A-CH,/No, o-determined from proton 


relaxation times in methane gas (ref. 56). 
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FIGURE 7. 


Comparison of the symmetrical part of the V3 band shapes 
of methane calculated from the semiclassical J-diffusion 


model with the band shapes calculated from the classical 


extended J-diffusion model 


Semiclassical J-diffusion - , classical J-diffusion ----. 


The angular momentum correlation times, T are (a) Jit 


ans 
10° ¢?, * (By FORO etc) 259) 10) <>) anal cy eee 


Togs sec. The frequency scale is in reduced units of 


= 1 ; 
27C VY (I/kT) @ where I is the moment of inertia of CH, and 
c is the speed of light. The band shapes are area normal- 


ized as in Fig. 5. 
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greater than 2 x Tomes sec. For correlation times less 


ehanvoex 200 


sec the semiclassical spectra begin to 
narrow more rapidly than the classical spectra. It 

is interesting to note that the classical model produces 
the "better" fit for high pressure samples. One must 
conclude therefore that the asymmetry of the spectral 
band is due to quantum mechanical effects, but that the 
symmetric part of the band contains enough information 
about the reorientation of the molecules to permit study 
using the simple cosine Fourier transform of the band 
shape. This procedure has been applied by a number of 


4,6,38,39,46 


workers with little justification and the 


results presented here establish its validity. 


Ze Limitations of the Theory 

In our semiclassical calculations we have used 
intrinsic line amplitudes appropriate to a rigid 
spherical rotor. °° Since vibrational-rotational coupling 


40,58,60 this rigid rotor approx- 


is important in methane 
imation requires some justification. In Fig. 8 we have 


compared the v, infrared spectra of CH, determined at 


3 
a density of 0.1 amagat, with the stick spectra cal- 
culated using the rigid rotor approximation. We see 
‘that although therigid rotor approximation predicts 


the R-branch to be more intense than the P-branch the 


reverse is true in the experimental spectrum. This 
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observation has been reported previously by Welsh, 
Pashler and Dunn. °4 However, the intensities of the R- 
branch lines are greater than those of the P-branch at 
higher densities (see Fig. 5), which is more compatible 
with the rigid rotor line strengths. There is, there- 
fore, little justification in using intrinsic line ampli- 
tudes other than those of the rigid spherical rotor un- 
less a more exact treatment of the vibrational-rotational 
interaction is to be attempted. Because of this failing 
in the calculated intrinsic line amplitudes, the appli- 
cation of this theory to the band shapes of CH, in low 
density mixtures may be of questionable significance. 

We have used the extended rotational diffusion 
picture of a collision to define the collisional trans- 
fer matrix, ne rather than semiclassical scattering 
theory because it drastically reduced the computation 
time required to calculate a band shape and thus makes 
the theory practical for "every day" use. In addition, 
we felt that the application of scattering theory with 
an ill-defined intermolecular potential merely served 
to obscure any information that could be realized from 
the experiment. Gordon et al. a6 have fitted the infra- 
red band shapes of the carbon monoxide molecule using 


“the semiclassical nan approach. They found that the 


widths of the vibrational-rotational lines predicted 
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by their calculation were 23% larger than the observed 
widths of low density gas mixtures. This discrepancy 
was attributed to either inadequacies in the intermolecu- 
lar potential used in their collisional cross section 
calculations ortopossible errors in their trajectory 
calculation routines. At higher densities, the agree- 
ment between their calculated and observed band shapes 
was no better than ours. Clearly, from the success of 
our simple theory relative to the success of the more 
sophisticated treatment of Gordon, we must conclude that 
until more accurate intermolecular potentials become 
available, the simple theory is adequate for the inter- 
pretation of band shapes in terms of intermolecular 
interactions. 

The semiclassical theory developed in section II- 
B is restricted toband shapes which can be described by 
a Limited «number: (<40)' of transitions. . This restriction 
is due to the need to diagonalize a complex matrix whose 
order is the number of transitions included in the cal- 
culation. For the V3 band shape calculations we included 
all transitions arising from rotational states with 
total rotational angular momentum quantum number J<12. 
Although this included approximately 96% of the CH, 
‘molecules inthe system, the predicted band shapes were 


still inadequate in the wings of the spectrum. Coriolis 
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coupling is much more important in the V4 mode than the 


V3 mode of cH The larger Coriolis splitting in the 
a. mode invalidates the averaging of the frequencies of 
the lines within a branch arising from each J-manifold 
which we have used to simplify the V3 calculations. 
Therefore, in order to include the same number of ground 
rotational states as in the V3 study, approximately 200 
transition lines would have to be considered. This is 
far too great a number to handle efficiently by the 
present procedure and still maintain reasonable computa- 
tion time. Similarly the application of the theory to 
more massive molecules would be impractical because 
larger numbers of rotational states are populated. 

In chapter III we will develop a "true" semi- 
classical version of the extended rotational diffusion 
theory. This theory will not be so severely restricted 


to a limited number of rotational states and will be 


used to study the V4 band shape. 
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SEMICLASSICAL EXTENDED ROTATIONAL DIFFUSION THEORY 


CHAPTER III 
A. INTRODUCTION 


In Chapter II, we have developed a semiclassical 
extended rotational diffusion interpretation of Gordon's 
eae semiclassical band shape theory. This approach 
adequately described the infrared absorption band shapes 
of the V3 mode of methane in gaseous and liquid mixtures, 
however the application of the theory was restricted to 
band shapes which could be described by a limited number 
of transitions (<40). For this reason the V4 mode of 
methane, in which Coriolis coupling is important, could 
not be treated by this procedure nor could band shapes 
of larger molecules in which a greater number of rotation- 
al states have significant populations. This restric- 
tion to a minimal set of transition lines also resulted 
in very poor agreement between predicted and experimental 
results in the wings of the band shapes of the V3 mode 
of CH, and made the predicted correlation functions 
meaningless at short times. Because the classical 
theories cannot be used to interpret band shapes in 
which Coriolis coupling is significant and the previous 
semiclassical theory cannot be used to treat systems in 


which there are a large number of transitions, there 


is a need for a semiclassical model which is not 
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severely restricted by the number of distinct transi- 
tions included in the calculation. 

The classical extended rotational diffusion theory 
NEN a has been widely used to interpret spectral band 
shapes and magnetic resonance linewidths in terms of 


poate The versatility of this 


molecular reorientation. 
theory lies inthe fact that the individual steps in the 
molecular reorientation processes are followed in de- 

tail. Thus the theory is not restricted by the length 
of the rotational diffusion step as is the case for the 


Debye >* and perturbed free rotor models.°’ 1/28 


Because 
of the versatility of the extended rotational diffusion 
model one would expect that a semiclassical interpre- 
tation of this model would prove quite effective in 
treating band shapes in which Coriolis effects are 
important. 


ee have recognized the 


Recently Biioe. et al. 
potential of a memory function approach to the extended 
rotational diffusion model. Berne and Harp ei have 
shown that the reorientational correlation function 


G(t) and its associated memory function K(t) satisfy 


the Volterra equation 


ce 
G(t) = - faux ayo(e-w ' (III-1) 
re) 


where G(t) is the time derivative of G(t). Bliot 
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GY le 


were able to show that, in the J-diffusion 
limit of the extended rotational diffusion model, the 
memory function K(t) was related to Kap(t)s the memory 


function for molecules undergoing free rotation, by 


K(t) = Kap (t)exp(-t/t,) ' (1TT-2) 


wnere te is the angular momentum correlation time. DE 
the free rotation correlation function See and its 
time derivative Gap (t) are known then Kap (t) can be cal- 
culated using Eq. (III-l). Then Eq. (III-2) can be used 
to calculate K(t), the memory function for a system of 
molecules undergoing rotational diffusion. Once K(t) 
is known, G(t) can be readily calculated from Eq. (III-1). 
Therefore, using Eqs. (III-1) and (III-2), and the 
functions Gar (t) and carte), the reorientational cor- 
relation function G(t) can be calculated. All details 
of the free rotational steps are included in Gap (t) and 
Gap (t) and therefore all that is necessary to define a 
semiclassical version of the extended rotational dif- 
fusion model is to calculate the semiclassical free 
rotor correlation function and its time derivative. 

in’ section B.of this. chapter, Eq. (LII-2) will be 
derived and the necessary relationships required for 
the application of the memory function approach will 


be developed. In section C the numerical analysis 
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procedures used in the application of the theory are 
presented, while in section D the model will be applied 


to the calculation of the V3 and wre band shapes of high 


pressure methane/inert gas mixtures.>” Pinaliv, in 


section E, the t_'s obtained from this analysis are dis- 


J 


cussed and the present model is compared to the model 


presented in Chapter II. 


B. THEORY 
ice The Memory Function and Extended Rotational 
Diffusion 


The infrared absorption band shape can be expressed 


co >> >> 
* 
I(w) = Ref at exp[-tut]) 9, ilp.E(0) ls)  (£|p.B(t) |i) 
fe) Lek 


GELS) 


where 4 is the transition dipole moment vector and E(t) 

is the electric field of the applied radiation, |i) 

and £9) are the initial and final states of a vibrational 
ede. trier transition, OP; is the Boltzmann factor for 
state |i) and w is the frequency of the incident radia- 
tion. In Eq... (2iifs3) both § and E(t) are defined ina 
molecular coordinate system. The time dependence of 

E(t) arises from the reorientational motion of the mole- 


> 
cule since the direction of E is fixed in the laboratory. 


.€9 


ames yroeds any aa Es tous . 
bailqys od Lfiiw rabott ond a 58 i ea od ; 
abi i to ida casa bred 1 baw 98 Wobaatuo toy 

it yhinos ee sons yeep stead\sandtion @ 
~-ebb ous eleylette shit wout noaisiider ° x“? eds .f | 


fefonm ett od Dal tate ganna a bn | 


f fer tule [Ve ’ fot iy ’ 


Git ise, any @ wie nae * ie 
oie arR) | | 
i me “as iio a re | a 7. 
(oe bre todoaw + ind olowié noid tenant ads #i q sxoriw © z ) 


beck 


CEN aontenben | boblage ents Yo blez? viwioole ens ak 

(onotssadiv 2 te toxin foakt bhp tndthat odd eas (as on. | 

xt aos9e3 ptemet1oa edi af #9 vnokt bens feaokssso3 * 7 
~Si bs sngbiont erty te svomeupas? ott Bi we bah ee edate 
i nd bagett ots (3)8 faci ite dod seen oe at nolo _ 

26 “Spndbreget amis adn smpaage eaaniba i 2 sfuvoton 

mp Loon ant Yo aebjom: Lsoktagne sce itt sae eoekxs (3)8 | 


i 


ill ort mi BOxsD tag sh aod sexkb ed? eons olus 


neti 
i] ¢ ' : — : . ne i : hide 
u rs ws : Peep ; Se i wa Dh es oui 
4 iM ; “a f i. v i } mn ; ahh — , nt, i) = i} a ' 4 
, s - a f - & 7 ; J : ora c 0 , = . = aa 
: , je Ie ee : ‘ er te ie 


64. 


The electric field E(t), defined in the molecular coor- 

> 
dinate system, can be related to the electric field 2, 
of the laboratory coordinate system using the Wigner 


(1) 68 


rotation matrix, The properties of the {YJ mat- 


rices. are discussed in Appendix I. Using the ‘J matrix 
> lab — (1) i 
E,, (t) = En ae RG (III-4) 
m 


where E,. (t) is the k-th spherical tensor component of the 
-> 
first order spherical tensor (a vector), E(t), defined 


in the molecular coordinate system and Eee is the m-th 


> 
spherical tensor component of the spherical tensor, ey 
defined in the laboratory coordinate system. Q(t) repre- 
sents the set of Euler angles describing the transforma- 


tion from the laboratory to the molecular coordinate 


system. Introducing Eq. (III-4) to Eq. (III-3), we 


obtain 
co m 
Liu) = ref at exp[-ivt]) yin pe eee) bi ae onan 
_ if ky! 
m,m' 
(1) (2)* 
eXt(D .p t(ONE) af cl1D. , .(aieril4) (gtr 3) 


where P_, are the spherical tensor components of the 
transition dipole moment defined in the principal axis 


coordinate system (the coordinate system in which the 
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inertia tensor of the molecule is diagonal). For a spher- 
ical top, the inertia tensor is diagonal in all coordin- 
ate systems and the spherical tensor components of the 
transition dipole moment can be taken to be the transi- 
tion moments ofthe degenerate vibrational modes of the 
molecules allof which are of equal magnitude. It should 
be noted that in expressing I(w) in the form of Eq. (III- 
5) we have assumed that the transition dipole moment is 
independent of the vibrational-rotational state of the 
molecule. In Appendix I, it is demonstrated that the 
space average over the product of rotational matrices in 
Eq. (III-5) vanishes except when m =m and he = k. 


Thus Eq. (III-5) becomes 


¥ ; 2) tals, 2 
T(o)-= ref dt exp[-iwt]) 0; >, [eae Ee 
k,;m 
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The orientation of E is arbitrary so we average over 


all possible orientations to obtain 
ed eee os (1T1=7) 
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4 > 
where € is the magnitude of E Substituting Eq. 


(III-7) into Eq. (III-6) and omitting all constant multi- 
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plicative factors 


eee Ref ae exp[-twt]) (uf (2(0) 1 gq “(2120 ) 
k,m 


2) 


(TTI-8) 


and the associated correlation function is given by 


g(t) (t) > (2,00 [2(0) ] Dt" *ia1) .  (III-9) 


k,m 


In going from Eq. (III-6) to Eq. (III-8) the sum over the 
projection operator, BORGAL has been carried out and 
the statistical average indicated by the weighted sum over 
the expectation values (il |i) has been indicated by 
the angular brackets. 

In the remainder of this section, a classical 
description of reorientational motion will be used to 
predict the reorientational correlation function defined 
by Eq. (III-9), and the relationship between this cor- 
relation function and its associated memory function 
will be examined in order to derive the relationship of 
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the classical extended diffusion model has been chosen 
here, rather than a semiclassical approach, because the 


physics of the model is most clearly conveyed in clas- 
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memory function as a vehicle for simplifying the math- 
ematical apparatus used to derive the expressions for 
the reorientational correlation functions and spectral 
densities, not as a statistical mechanical function 

in its own right, we will develop the theory using the 
Classical description of reorientation and show that 
the memory function formalism leads to the same results 
as the detailed geometric approach does. The develop- 
ment of the semiclassical model will then be made using 
the memory function approach. 

The extended rotational diffusion model describes 
the reorientational motion of a molecule as a series of 
rotational steps which are terminated by "collisions". 
Two limiting types of "collisions" are considered: the 


DOI) Se pees ec ka 


M-diffusion and J-diffusion limits. 
fusion limit, the orientation of the angular momentum 
vector and its mati cude are randomized onto a Boltzmann 
distribution, while in the M-diffusion limit, only 
the orientation of the angular momentum vector is ran- 
domized by "collisions". In both cases, the reorienta- 
tional motion will be changed by the collision. 

If we divide the ensemble into groups of mole- 
cules which are in their n-th diffusive step at time t 
(i.e. have experienced exactly n-1 collisions at time 


t) then the correlation function of the ensemble can 


be expressed as 
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Ee lee pinta int), (ITT-10) 


n=l 


where 6) (nt) is the correlation function for those 
molecules which are in their n-th diffusive step at time 


er pin; ct) 2.s" given vy the Poisson distribution 


-l 
Gere Je exp [-t/t5] 
p(n,t) AL ae Sat Lease ee ea li 
(n-1)! 
(III-11) 
and Ty is the angular momentum correlation time which is 
the average time between "collisions". From Eq. (III-9), 


it is clear that the partial correlation function 


6) tnt) is given by 
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(IITT-12) 


where Q(t) describes the orientation of the molecular 
coordinate system of a molecule in its n-th diffusive 
step at time t relative toa lab-fixed coordinate sys- 
tem. 

We will now consider the classical reorienta- 


tional motion of a spherical molecule undergoing changes 


in its angular momentum due to collisions. It will be 
assumed that the collisions have occurred at times tye 
tor ae aye ead and that the molecule rotates freely 
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> 
about the direction of the angular momentum vector, J, 


with angular frequency, be given by 


> ee 
w = J/I (1Ei-13) 


where I is its moment of inertia. The matrix oN AED} 
describing the transformation from the laboratory coordin- 
ate system to the molecular coordinate system of the 
molecule in its first diffusive step at time t can be 
expressed as a product of D matrices in the following 


way 


(1) Leas! (1) 
De (he CES a a io By 
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i Ti uie,8, 088 > a) lawuie Oe 
(ITI-14) 


where D7) [24 describes the transformation from the 
laboratory to the molecular coordinate system at time 
zero, ST) Heo Vex8 51 describes the transformation from 
the molecular coordinate system, at time zero, toa 
coordinate system in which the z-axis lies along the 
direction of the angular momentum vector, D AYos, t, .010) 
describes the rotation of the molecule about this z- 
axis through an angle wit, and gy A aries ee ily % des- 
cribes the transformation back to the molecular coordin- 


ate system defined at time zero. 
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Equation (III-14) defines the transformation of a 
first rank spherical tensor from the laboratory coordin- 
ate system to the molecular coordinate system of a mole- 
cule’ which is iff its first diffusive step at time t. This 
transformation for a molecule which is in its second dif- 


fusive step at time t is given by 


(1) LP) (1) 
2" helt) ) = Do? (2, (t,)) DB" fo, 8570) 


0 DEL) pe (esevdi, ope Diaklalomeas.0) 


(1) (1) (1) 
DN 195) Di faye By OID *™” Tw, t,,0,0] 
Seva fat, Bama) NG Laer 0] 


(1) (1) -1 
x D [w, (t-t}) 0,0] D [a5,B,,0] ° 


ChIISb5) 


where Coty describes the transformation from 
the laboratory to the molecular coordinate systems at 


time t the -bimesotitchbe first coll3sion,..and 


nf 
=-1 

describes the reorientational motion during the second 

diffusive step in the period (t-t,)- Similarly, Lox 


a molecule which is in its n-th diffusive step at time 
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where TT indicates an ordered running product from 

i=l 
left to right in ascending order of i, and t; is the 
time at whichthe i-th diffusive step terminated. The 
angular momentum vector, we during the i-th diffusive 
step has an orientation (a, ,8;,0) with respect to the 
molecular coordinate system at time ea) and Was the 
angular frequency of rotation about the direction of 
ae is Given by Eqs (ILI-13)? In Eq. (Iit-16) t..and t 
are associated with 0 and t respectively. 

Equation (III-16) describes the transformation 
from the laboratory to the molecular coordinate system 
for a molecule which has undergone a particular series 


ofen.diffusive steps. The partial reorientational cor- 


relation function eo wares can then be written as 


(1) (1)* 
CL) a PD WS 
G (nt) -) (o [20] ae [20] 
k,m 


CL) on 1) 
x a [a,,B,,Olexp[tbw,t,] aoe [a5 B, 70] 


(SL+ TTR) 


sgt at Poa bas <t Qo ee ve lone 1 


ot hotsghess3 ‘qote eM, au 


avied?tsb sit ort opty «30a 7 


arid oa Sonaua’ dodo’ (Ody) ibe 
BAS». a Dee: Lage once mm i Tea taeen ea ex 
io notion ocd esti ¥ox 8 yoneuRewT side ae 


of See 41 ( a AEE), * alt OR co ‘inayat ak iF a 
acta iv 


i 


aavaye -adantsno99: BENDS Loui 9 
csiuotsiea ‘ae 101 y Luo 

“109 feaokak denice. savinaid en penn we 

as rotseaw ‘od Peed js” ea anes — roaster 


726 


(1)* (1) 
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in which advantagewas taken of the relationship cS 


Df) [6,0,0] = exp[-tad] ae (T7118) 


and 


(1) * 
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bl fs 
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In Eq. (III-17) and subsequent equations the Einstein 
convention of implied summation over repeated indices is 
used. 

The ensemble average over the initial Sree nee= 
tions i is independent of all other averages and in 
the absence of a strong applied field all orientations 


are of equal probability so that Se 
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The angles a OF Bye PRES Pe cay Bat which define the 


orn Dy 
orientation of the angular momentum vector in the 
molecular coordinate system at the beginning of each 
rotational diffusion step, are randomized at each col- 


lision, thus statistical averages can be independently 


taken over each set of these angles. This statistical 


average is given by oe 
27 1 (1) a. Gh)* 
-1 ; WD 
(4 7 ) fs ac f/ dg sin 8 Dany [a,B,0] m5™ [a,8,0] 
fe) fe) 
z 
== (IITI=-21) 
3 m,/m., 


Applying these relationships to Eq. (III-17) we obtain 
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(III-22) 
where the angular brackets imply a statistical average 
over tyreest sy and Wyre Ws 

In the extended rotational diffusion model the 
times between collisions are assumed to be independent 
of the angular momentum state of the molecule, thus the 


averages over the collision times t,,...,t,_,can be 


taken independently of the average over the angular 
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Wipreeer wo: It will be assumed that the 


probability that a molecule has experienced exactly n-1l 


frequencies, 


1! tor t. + oh 


: : 5 - 
a sctokay a dt 4 is given by (n=)! dt ,dt,...dt,_,/t 


Taking the statistical average over the collision times, 


collisions in the periods t ty dts 


t 


we obtain 
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a 2 
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G (rhe) it fa, ce: a gee f dt, 
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n 
a : 
a pees CATT? 
(III-23) 


In order to perform the statistical average over 
Wpreee OW, it is necessary to specify either M- or J- 
diffusion. In the M-diffusion limit, the magnitude of 


the angular momentum vector is not changed in a collision 


therefore Wy = Wy = veers = WL Ew and Ed= (1li-23)ebem 
comes 
te Le 
2 n 
er (n-1)! (1) s: 

(nie) vferion --f dt, af Con eet ago) 

fe) fe) i=l 
W 
(III-24) 


where the angular brackets indicate the statistical 


«Sf 


[on yloosxe covets il | 
tener + gee pdt ye ee 
aie | bens vee ‘en 


ee 


be ae, eel , 


| 


Revo SPEISVE cabot ab motred S tote | LS a 
| #h we “4 uaddto Setiad ccoeen ah a eee P 


o. 


to abusimpén oat petty odd red | ase - j 
aotail tos ‘eon bopnatto jon nd fee musitomore tect: 
od, eset) pa bike we ae Fie sk ee” 


‘ 
i . Pees cae 


13% 


average over w and the free rotor correlation function 
(1) : ; 
Gar [t,w] is given by 


(1) oe i 
Gor [t;w] = eg ree ; (TTIT=—25) 
In the J-diffusion limit, the magnitude of the angular 
momentum is randomized ontoa Boltzmann distribution at 
each collision and the statistical average over w must 
therefore be performed independently for each diffusive 


step. Therefore in the J-diffusion limit, Eq. (III-23) 


becomes 
€ t. 
(1) (n-1)! TU CL) oc 
G (n, ©) | J dt. soo) A Goepplt t; i! 
re) fo) 
(III-26) 
where the free rotor correlation function G (1) (4) is 


J:FR 
given by 


G,, aS te) ED) Ds, exp [tawt] » 2 (III-27) 
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The reorientational correlation functions are 
then obtained by substituting Eq. (III-24) or (III-26) 
intosEq.. (TI1=10) sand.applying Eq:, (TiPs11). .Fhus 


t 


= =r 
ea (t) 3 expt-t/%51 3. ei os eeee 
n=l 


Oo 


relpas att to wicsledaan we: ae he " 
os nokstud tages di gl ctr Seeatees “y 


i ~ a) =< 


Jetin wo seo SBA TSVE, | 
ovieva2eb dose a2, | 
(OS ohne)! ope, gems 


$ 


CL): eee 
Te ae ee Cho 


(9S—UET} 
ei (4 ed ae? scsi 


(VS=2TTT i ( [idistigxo EE ar 


pe eT ngid storia teqels Bs #8 
(a8—-1Zz)’ x0 (sget) po! Me cdi, yet 
Be: tein Pa ype xh si el s112) .p3 edad ” 


16. 


i 

x (1) * 

Tes! TU op bey taal] ) (III-28) 
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G (t) = exp[-t/t,])) Ts (ices Geka torare 


n=1 (e) 
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The memory function K(t) can be conveniently expres- 
sed in terms of the correlation function using Laplace 
transforms. By taking the Laplace transform of Eq. (III- 
1) and applying the convolution theorem 6? 8 Laplace 


transform theory we obtain 
(III-30) 


where K(s) and G(s) are the Laplace transforms of K(t) 
and G(t) respectively and s is the Laplace variable. 
Using the extended rotational diffusion expression for 
G(t) and Eq. (III-30) we can obtain expressions for 
K(s) in the extended rotational diffusion model. First 
consider the J-diffusion limit. Taking the Laplace 


transformof Eq. (III-29) we obtain 
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where Sema, is the Laplace transform of one, and 
hap ema g + t, 7) is the Laplace transform of 

1 (1) . 
exp ( t/t 5)G5. pr Cesc: The convolution theorem of 


Laplace transform theory was used in deriving Eq. (III- 
31) and the infinite sum rule for geometric progressions 
was used to obtain the second expression on the right 
hand side. “Substituting Eq. (Lit-31) into Eq. (Iii-30) 


we obtain upon rearrangement 
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G(s + lg 
= (1) , C1) 
where Ks (s) is the Laplace transform of Ks (t), the 


memory function in the J-diffusion limit. Using Eq. (III- 
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30) the memory function associated with CsePR (tf) tis 
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Comparison of Eqs. (III-32) and (III-33) we see that if 
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we associate s' with s + T.~ and thus neded (eh with 


J 

KA) (ty exp[-t/t5] then the equations are the same. 
This is the relationship given in Eq. (III-2). Bliot 
65,66 obtained this relationship [Eq. (III-2) ] 
using statistical mechanical arguments and they demon- 
strated that it defined the extended J-diffusion model. 

In the M-diffusion limit the treatment is more 
involved because only the orientation of the angular 
momentum vector and not its magnitude is changed in a 
"collision". Thus molecules with different angular fre- 
quencies are treated as separate species and there will 


be a separate memory function equation for each species. 


The development in the M-diffusion limit is analogous 
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order to calculate the ensemble function Gu 


instead of G Ct). and Ku 
(i). or 
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From the above development it is clear that all 
that is required to calculate the classical reorienta- 


tional correlation function for the extended rotational 


diffusion model is the free rotor correlation function 
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the expressions are determined by the description of the 
collision process and not by the description of the free 
rotation step. Therefore the theory can be modified to 
yield a semiclassical model simply by replacing the 


classical free rotor correlation functions Gipe eebt) and 
Bott) 


FR 


tions. In the next section, quantum mechanical free rotor 


(t,W) by the corresponding quantum mechanical func- 


correlation functions will be considered and in later 
sections these functions and their time derivatives. will 
be used to calculate spectral densities and reorienta- 
tional correlation functions for a semiclassical extended 


rotational diffusion model. 


ae Quantum Mechanical Free Rotor Correlation Functions 

The relationships developed in the previous section 
indicate that a semiclassical extended rotational dif- 
fusion model is readily developed provided one can cal- 
culate a quantum mechanical free rotor reorientational 
correlation function and its first time derivative. In 
this section, these functions will be considered. 

From Eq. (III-9) we see that the quantum mechanical 


free rotor correlation function is given by 
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where p, is the Boltzmann factor for the state |i), Shs 
represents the set of Euler angles describing the trans- 
formation from the laboratory to the molecular coordinate 
system ss of a freely rotating molecule at time t. The 
sums over i and f are over the initial and final transition 
states of themolecule and the sum over k is over the degen- 
erate vibrational modes of the vibrational transition of 
interest while the sum over m is over arbitrary directions 
in the laboratory coordinate system and can be ignored since 
it results simply in multiplication by a constant factor. 
The time dependence of D ) 10(t) ppl is determined by 
the motion of the free non-interacting molecule. Since 
AO) tngeeess is a Heisenberg operator, its explicit time 


dependence is given by 
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ide, D (1) =tent 
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D)) factrep) = expt Lt (2(0) lexp (=<) 


(Tit=35) 
where? is the Hamiltonian of the system. In the case of 
a rigid rotor, +/ is the rotational Hamiltonian. The 
free rotor correlation function can therefore be written 


in the form 
gS 
e'2) (e) pet ad, Kil D. 120) [£> | exp 2u, MM cet 


where Oss is the frequency of the transition i > f. 


gq t2)2's (i) etete eas tok we 
+EnER2 eric pnldinoae nen 


off -.3 eaiz ts otwonla'g 


2m moisiansis. snaoicsesdiy par " Rey 
enoltua tt 11639 SES 190. sm so im sot 


. yahoos 2 laiatheaitibal & or 
«eed bonterssas ak ae 


atin: siuoabont nye ) 2 


oy clei ane 


| Ae hee » ome | 


“bio, Opes i? : i a, 


(eé-r00) 


ae 5 ied wodtionns: ae to ‘eniopend ots at aes ey 


’ ¥ 
_ . i as > ; A 
@ be al A 7 Bas 
J o] a¢ & % i - uy Lice. ee vi 
- ee F f we ae, vee = 


81. 


It is clear from Eq. (III-36) that one requires only 
a knowledge of the ground vibration-rotation state Boltz- 
mann factor Pir the transition moments <i] ta(0)1/£> 
and the transition frequencies Wee of the molecule in order 
to determine the quantum mechanical free rotor correla- 
tion function from which I(w) and G(t) can be determined 
using the methods suggested in the previous section and 
which will be discussed in detail in the following sec- 
tions. .In the next section abrief description of the rota- 


tional states of molecules is presented and the semiclas- 


Sical extended rotational diffusion model is discussed. 


Sy Description of the Rotational Motion 

In a classical rotor, the angular momentum vector 
can take on any one of an infinite number of orientations 
with respect to a molecule-fixed coordinate system and 
with respect to a space-fixed coordinate system. The 
Magnitude of the angular momentum vector can vary con- 
tinuously from 0 to ~», with the probability of a given 
Magnitude and orientation following a Boltzmann distri- 
bution. ® In the quantum mechanical description of 
rotating molecules, the orientation and magnitude of the 
angular momentum vector are no longer continuous 
variables as they are in a classical top. There exist 


distinct rotational states, which we shall denote 
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| K M>. The magnitude of the angular momentum vector 
for a molecule in the state |J K m> is /g(g + 1)f, 
while the projections of this vector onto the space-fixed 
z-axis and onto the molecular symmetry axis are Mh and 
Kfi respectively. The quantum number J is a non-negative 
integer and the quantum numbers M and K can take on in- 
tegral’ values from +J to -J. It is clear that K and M 
define the spatial orientation of the angular momentum 
vector and J determines its magnitude. The equilibrium 
populations of the rotational states follow a Boltzmann 
distribution just as the orientation and magnitude of 
the angular momentum vector in the classical top does. 
In the absence of electric and magnetic fields, all 
rotational states with the same J and K, but different 
M, are degenerate. In linear molecules, the angular 
momentum vector is perpendicular to the molecular axis 
so that K = 0. For a symmetric top molecule, the rota- 
tional states with different values of |K]| are non- 
“degenerate, and their energy differences depend on the 
relative magnitudes of the moments of inertia about the 
symmetry axis and about axes perpendicular to it. The 
rotation states of a spherical top molecule, with the 
same value of J but different K values are degenerate, 
because of the inertial isotropy of the molecule. 


In the simplest version of the classical extended 
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ore. S749 a "collision" was 


rotational diffusion theory 
interpreted as an event during which the angular mom- 
entum vector of a molecule was changed. In the J-dif- 
fusion limit, both the orientation and magnitude of the 
angular momentum vector are randomized onto a Boltzmann 
distribution at every collisional event. In M-diffusion, 
the magnitude of the angular momentum vector is un- 
changed by collisions, and only its orientation is 
randomized. By analogy to this classical extended ro- 
tational diffusion picture, we suggest that, in a semi- 
classical version of this model, "collisions" be con- 
sidered as events wherein molecules undergo transitions 
from one rotational state to another due to perturba- 
tions from neighboring molecules. In the semiclassical 
J-diffusion model, "collisions" will be assumed to 
randomize all three quantum numbers J, K and M onto a 
Boltzmann distribution of rotational states, i.e. the 
probability of a transition to a state |J K M> 

brought about by a "collision" will be determined by 

the Boltzmann paseeuwisy that state and will be inde- 
pendent of the initial rotational state. Only transfers 
between levels with the same value of J will be allowed 
in the M-diffusion limit, where the quantum numbers K 
and M will be randomized by collisions. It is clear 
that in the M-diffusion limit, the molecules in the 


ensemble are subdivided according to their total 
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angular momentum quantum numbers J and remain within such 
a subdivision throughout the diffusion process. The 
molecules within a given J-manifold will therefore be 
characterized in the M-diffusion limit by a unique re- 
orientational correlation function and spectral band 
shape. The ensemble correlation function and band shape 


will, of course, be weighted sums of these components. 


4. Spectral Densities and Correlation Functions in the 


M- and J-Diffusion Limits. 

The Fourier relationship [Eqs. (1-5) and (I-6)] 
between the spectral density I(w) and the correlation 
function G(t) is conveniently expressed using Laplace 


transforms as 
ie Re { G (<u) } : (III-37) 


where G (tw) is the Laplace transform of G(t) and tw 

is the Laplace variable. Equation (III-37) takes advant- 
age of the time symmetry of G(t), i.e. G(t)* = G(-t). 
Equation (III-37) will prove a convenient means of 
calculating I(w). 

The computation of spectral densities and reorienta- 
tional correlation functions in the two limiting cases 
‘of the extended rotational diffusion model differ only 
in the manner in which the ensemble average over the 


rotational states of the molecule is performed. In 
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the J-diffusion limit, where collisional events effect 
transfer of molecules between all possible rotational 


states, the appropriate free rotor correlation function, 


(1) 
CS:FR 


over i encompasses all rotational states, | A: JKM > 


(t), is given by Eq. (III-36) in which the sum 


of the lower vibrational level. The Laplace transform 


(1) ; 4 
of GCsePR (t) is given by 


(1) 2 
e(t } - : Sa f»| 
gO) Gu tot, 2) = } rm) > Neate a [2(0) J |£ >| 
JKM f m,k 


see eh re ha Wye, ¢) Ty! > 


(Tet=36 


where ® sem is the Boltzmann factor for the state 

|A:JKM> , SK, € is the transition frequency of the 
transition |A:JKM>—~ |f> and A is a vibrational 
quantum number. The band shape in the J-diffusion 


limit [see Eqs. (III-31) and (III-37)] is 
I (w) = Re {rsas(w)/t1-a; (wo) I} ' 


(ITI-39) 


where 
: (III-40) 
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by the memory function procedures to be described in 
the next section. 

In the M-diffusion limit, as indicated in the pre- 
vious section, all molecules with a common total angu- 
lar momentum quantum number J are characterized by a 
unigue reorientational correlation function 6) (ts) 
and a unique spectral density I(w;J). The M-diffusion 
correlation function, G alee ye and spectral density 


M 
Ty 6) are related to these J-manifold functions by 


(1) (1) (III-41) 
and Ou ie > ca At 
5 
1 @) =) 0, 1 3) 
J (III-42) 
where 
ay > PIKM 
K,M 
(III-43) 


The free rotor correlation function appropriate for 
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6) (233) =} > 1 ¢ ase] DY taco le >|? 
K,M a fo ek He 

x exp (Wy, et) , 

(III-44) 
where 

a = PD talp 

KM JKM "J 
(IITI-45) 
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The Laplace transform of Gar 
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(IITI-46) 


and the M-diffusion spectral density, Ty) « is given by 


(III-47) 
ACD) oe pytyRe { By(u)/ 12 - B,(w)}} ; 
a 
where 
Bi (w) = T, age 6 Cu cf ce 3J) 


(III-48) 
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One could calculate the M-diffusion correlation function 


(1) 
FR 


the corresponding free rotor memory functions ee EE 


for each J-manifold, but this would be rather impractical 


from the free rotor correlation functions G (t:;J) via 


if the number of J-manifolds were large. On the other 
hand, one can compute Ty) from Eq. (III-47) with little 
difficulty, even for a large number of J-manifolds, and 
the correlation function can be computed by Fourier 


transformation: 


gi) (t) = fi 1, @) exp (-iwt) dw 


—c 


(III-49) 


C. NUMERICAL ANALYSIS PROCEDURES FOR IMPLEMENTATION 
OF THE MEMORY FUNCTION EQUATION 
The computation of a reorientational correlation 
function using the Memory Function approach consists of 
two steps: calculation of Kar (t) from Gar (t) and 
ome and calculation of G(t) from K(t). The first 
calculation is carried out using Eq. (III-1). Once 
Kap (t) is known, K(t) can be calculated for a given Ts 


using Eq. (III-2). Equation (III-1) can then be used 


to calculate G(t) from K(t). In this section the 
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numerical analysis procedures used to perform these cal- 
culations are described. Throughout this section, no 
distinction will be made between M- and J-diffusion since 
they differ only in the definition of Gap (t)- In the 
M-diffusion limit the G(t) of the ensemble can be cal- 


culated using Eq. (III-41). 


Le Calculation o£ Kap (t) from Gap (t) 


Garit) and Gop (t) can be readily calculated from 
Eq. (III-36) or Eq. (III-44) depending on which limit is 
desired. The numerical procedure nN for generating 
Kar (t) for t = 0, A,---,nA,----- , is based on a trapezoid- 
al approximation to the integral in Eq. (III-1) and uses 


the relations 


Re TO) suse etl. (ITT 50) 


nx 
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I 


Pig Yo ¢ 
ACpRlAl = Kap l0jG., [A] ; (ITI—51) 


and 


K p(n] = F Gpglnd] - Kyp 0169 {nd 


n-l 
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(Lit=—52) 
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2 Calculation of G(t) from K(t) 

The procedure for computing G(t) at t = 0, A, 2A,--- 
---nA,----- , given K(t) at these values of t, is based on 
trapezoidal approximations for the integrations in Eq. 


(III-1). given by, 


EipAlis -a #016 nt) Pek (AVGl n=l + Seek nel)G1CiAl 
2 
, K{nalctol 
2 (III-53) 
and in the equation 
nA 
G[nA] = G[(n-1)A] + du G(u) 
(n-1)A 
= 'Gl (nella he es }Gtna} * ét(n-2) a) : (ITI-54) 
2 
In addition, the conditions 
Gi0)/"=1"", (III-55) 
and 
Glo] = 0 (III-56) 


are required. The first is just a normalization condition, 
while the second determines the phase of G(t) and is 
required by the theory [see Eq. (III-1)]. 


Substituting Eq. (111-53) into™Eq. (I1I-54) and 
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imposing conditions (III-55) and (III-56) we obtain 


he A2 
G[A] = 1 - “xia | 1+ — stor} ; (TII=57) 
4 4 


and from Eq. (III-54) we obtain 


G[(n-1)A] = -G[(n-2)A] + £(G[(n-1) A] - G{(n-2)A)) 


CE LE-53) 


Substituting Eq. (III-53) into Eq. (III-54) and solving 


for G[nA] we obtain, 


sete 2 
G[nA] ={el(a-1)4) A = Gite ie : K[nd]- 5 K{(n-£)A}c [24] } 


(TLit=59 ) 


Therefore given the functions Grp (t) and ee(G yy G(t) 
can be calculated by simple numerical procedures involv- 
ing the intermediate memory functions Kap (t) anu) LE 
should be pointed out that Kap (t) need be calculated 
only once and can be employed to compute K(t) and G(t) 
.for the required values of Ty in later calculations. A 
listing of the Fortran programs used to calculate 


Kap (t) from Gap (t) and G(t) from Kap (t) are given in 


Appendix IIB and IIC respectively. 
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ae Investigation of the Stability of the Computation 


Procedures. 


We have employed the trapezoidal approximation in 
the numerical solutions to Eq. (III-1), rather than a 
more sophisticated numerical integration formula. In 
test calculations using Simpson's rule formulae and- 
higher order Newton-Cotes approximations, it was found 
that the calculated correlation functions and memory 
functions were much less stable than in calculations 
performed with Eqs. (III-50) - (III-59) above. This 
numerical integration procedure is also much more versa- 
tile than the power series expansion procedures used 


G27 S5 This procedure consisted of expanding 


previously. 
K(t) and G(t) as power series in t and relating the 
coefficients in these series using the Laplace trans- 
form relationships. This procedure broke down at times 
greater than 4/YI/kT seconds even when 100 or more terms 
were included in the expansion. The present procedure 

is stable at times as large as 20/ I/kT seconds and, if 
longer times are required, it is necessary to reduce only 
the size of A, the time increment to increase the range 
of stability. In our calculations, we were interested 

‘in times in the range 0 to 10/ I/kT. For ti is range of 
times, we found that a time increment of 0.02/ I/kT was 


more than adequate to assure numerical precision and 


accuracy. 
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The stability of the numerical methods were evalu- 
ated using the following procedures. The memory function 
calculated from the free rotor correlation function 
using the methods developed in section III-C-1 was 
compared to that calculated using the power series ex- 


pansion technique oye 


(which is exact for small t). 
This test indicated that the trapezoidal procedure is 
accurate over the range in which the power series 
approach is reliable. The generation of Kap (t) from 
Gap (t) was also tested by applying the numerical pro- 
cedure to the "correlation function" defined by 

(1 + t/2)exp(-t/2) whose associated "memory function" 
is defined by .25exp(-t). The values of K(t) calculated 
from this G(t) agree with those calculated from the 
closed expression .25[exp(-t)] to four significant 
figures at t = 20 y1/k?t. 

The numerical procedure for calculating G(t) 
from K(t) was tested by calculating G(t) using the 
methods of section III-C-2 when, Ty = o, The 
agreement between the function calculated in this way 
and Grp (t) indicated that this procedure is very 


accurate over the range 0 £ t < 20 /I/kT. 
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DD. APPLICATION TO V3 AND Va BAND SHAPES OF CH, 
13 Calculation of Band Shapes of CH, 


In section III-B-3, we described the properties of 
the rotational states of molecules, pertinent to con- 
Siderations inthe rotational diffusion model, and devel- 
oped these ideas in section III-B-4 into formulae for 
the spectral densities and reorientational correlation 
functions for molecules whose rotational behaviour was 
approximated by the J- and M-diffusion models. For 


spherical top molecules like CH certain complicating 


4’ 
features arise which affect the populations of the 
rotational levels and the vibration-rotation coupling in 
the upper vibrational state, and these factors must be 
taken into account. The first complication arises be- 
cause the nuclear spin wavefunctions of the equivalent 
nuclei (e.g. the protons in CH,) combine with the 
rotational wavefunctions in only certain ways, in order 
that the total molecular wavefunction satisfies the 
Pauli principle. Molecules which have rotational wave- 
functions which combine with a given set of symmetry- 
adapted nuclear spinwavefunctions are often designated 
-nuclear spin isomers. In the case of CH, where the 
four protons each have nuclear spin $i there are three 


nuclear spin isomers: ortho-CH, (rotational and nuclear 


spin wavefunctions of symmetry F) , para-CH, (symmetry E) 
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and meta-CH, (symmetry A), which correspond to nuclear 
spin triplet, singlet and quintet states respectively. 
The expected 2J + 1-fold degeneracy of the J-manifolds 

in the spherical top due tothe 2J + 1 allowed values of 

K must be modified since a different number of nuclear 
spin wavefunctions can be combined with the various K- 
states within the J-manifold. The resulting degeneracy 
varies with symmetry type and has been calculated by 
Wilson oc using group theoretic methods. The symmetry- 
adapted rotational wavefunctions [OT M>, where I desig- 
nates the irreducible representation and s is a serial 
index to distinguish linearly independent functions be- 
longing to the same irreducible representation, are 
linear combinations of the functions |JKM> with different 
values of K, but with the same values of the quantum num- 
bers J and M. 

Since the rate of interconversion of nuclear spin 
isomers is believed tobe slow compared to the molecular 
reorientation process oy, the spectral band shape is 
really a superposition of the contours for each nuclear 
spin species. This implies that calculations for each 
of the isomers must be performed separately and the 
“results combined at the end to produce the composite 


correlation function. 


If one neglects very small effects due to centri- 
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fugal distortion, all rotational states [oT .M> in a J- 
manifold of a spherical molecule in its ground vibra- 
tion state will be degenerate, but in an excited vibra- 
tional state in which one or more degenerate vibrational 
modes are excited, the degeneracies are lifted by the 
interaction between rotational and vibrational angular 


oe This Coriolis interaction causes the degener- 


momenta. 
acy of the states with different quantum numbers rs but 
identical quantum numbers J and M to be removed, and 

each of these rotational wavefunctions combine with each 
of the vibrational wavefunctions to produce as many 
distinct vibration-rotation states with a common symmetry 
label ar as there are components of the degenerate 
vibration. The magnitude of the Coriolis coupling there- 
fore varies with both rotational and vibrational quantum 
numbers. If all rotational transitions, AJ = 0, +1 
accompanying vibrational excitation were allowed, we 
would expect several lines to arise from each ground 
state molecule with rotational wavefunction [OT M>. 
However, a detailed analysis of the selection rules 20,72 
indicates that only one transition is allowed in each of 
ptherP, QO and R branches of the infrared spectrum of a 
tetrahedral top with ground state rotational wavefunction 


|JT.M>, and the magnitudes of the transition moments 


of the transitions ineach branch arising from states 
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within a given J-manifold are equal. The transition 
moments are given by Childs and ganinf2 
In section III-B-3 and III-B-4, we have expressed 

the free-rotor correlation functions and their Laplace 
transforms in terms of the transition moments 

| <A;JKM| PD aarais(0). 12>) 4 From the above discussion, 
conclude that all transitions arising from states 
|A;JKM>, or more correctly |A;oT .M> , Within. aysingle. J- 
manifold have transition moments ao which depend only on 


the value of AJ for that transition. For AJ = +1 (R- 


branch) 


S_[<asarga] DB © aco |e>|° = Gs + /QI +0 
™m 


m,k 
(III-60) 
for AJ = 0 (Q-branch) 
(1) 2 
>, [Casa Pe (MCOIlEP] =, 
m,k 
(III-61) 
and “for “AU = —Pr(P=branch) 
(1) é 
> I< Asst M|D 3 [aco] |£ >| = (25 -1)/(25 + 1) 
m 
m,k (III-62) 
’ big 59,78 
With these transition moments and the observed and 
Caieataces transition frequencies in each vibration- 


rotation band, one can compute theoretical reorienta- 
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tional correlation functions and spectral band shapes 
which can be compared to those obtained from experimental 


39 
measurements. 


a. Vv Band 


3 
The transitions in the V3 band of CH, are not strongly 
influenced by Coriolis effects because the Coriolis' 
coupling constant is very small 20 so that the transi- 
tions within a branch of the spectrum which arise from 
states within a given J-manifold are not widely separated 


in frequency (<2 cm? 


). The frequencies of the lines 
originating instates in the J = 0 through J = 16 mani- 
folds were obtained from the experimental results of 
Pivylervet®aly > eh and it was found that calculated cor- 
relation functions and spectral densities based on fre- 
quencies calculated from theoretical formulae 20 which 
do not include splitting of the lines from | oT M> LOE 
different r. were essentially identical to those based 
on the exact experimental frequencies. Since the 
experimental frequencies were known only to J = 16, we 
have used calculated frequencies for lines originating 
in the J = 17 through J = 25 rotational levels. As 

-we found in Chapter II (see Fig. 2) the spectral band 
shapes and reorientational correlation functions computed 


as superpositions of the functions for A, E and F nuclear 


spin isomers, were virtually indistinguishable from the 
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band shapes and correlation functions calculated by 
allowing "collisions" to bring about nuclear spin isomer- 
ization as well as rotational state transitions. The 
results for V3 reported below are from calculations in 
which the transition frequencies involving the rotational 
states of all symmetry species are included in a single 
computation, rather than separate calculations, each 
including only the transitions arising from rotation 
levels of a single symmetry. 

The free rotor memory function for the V3 band of 
CH, is shown in Fig. 9. It is to be noted that this 
free rotor memory function is complex because the transi- 
tion frequencies and intensities are not symmetrically 
disposed about the frequency of the band center and 
that the real part of the memory function closely re- 
sembles the classical free rotor memory function for a 
spherical rotor \ The zero of frequency for the tran- 


Sition frequencies w ¢ was chosen so that the free 


ap laa 
Ss 
rotor correlation function at t = 0, Care oe was real 
and equal to unity as required by Eq. (III-1l). The 
(1) 


free rotor correlation function, CsePR (t)))o whieh 1s 
also complex, is also shown in Fig. 9. 

The J-diffusion reorientational correlation func- 
tions eos at 295K for various values of the angular 
momentum correlation time, Tyr shown in Fig. 10 were 


calculated from Kap (t) weingenos. (Lil=53), = (1Tiri—59). 
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FIGURE 9. Free rotor correlation and memory functions 
for the V3 band of CH, at 295K. Time scale 

: 5 1, 
is inreduced units of t(kT/I)* = 8.7 x 10?t. 


K(0) = 1.95. 
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FIGURE 10 


J-Diffusion reorientational correlation functions 


band of CH, at 295K for various 


(1) 
G5 (t) for the V3 4 


values of Ts: 


-re{c!") (t)}, ----- imc}? (t)}. For Clarity, 
successive correlation functions have been displayed 
vertically by 0.2 units. Time is in reduced units 


as an thig.. 9. 
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The spectral densities L,() at 295K for the same values 
of T, are shown in Fig. 1l. 

The spectral densities Ty () at 295K, for various 
values of the angular momentum correlation time Tye in 
the M-diffusion limit were calculated using Eqs. (III- 
46) - (III-48) and are shown in Fig. 12. The correspond- 
ing reorientational correlation functions Gere) were 


obtained by Fourier transformation of the spectral den- 


sities and are given in Fig. 13. 


Db. V4 Band 


Coriolis coupling is very significant in the V4 band 
of methane, and the splitting of the lines within a given 
branch arising from molecules with ground rotational 
states [ST M> within a given J-manifold is very large, 
often large enough that the lines from different J-mani- 
folds are interspersed: /7 The frequencies of the lines 
in the v, branch were calculated from the formulas given 


4 
by Moret-Bailly toa. 


although the off-diagonal ele- 

ments of the vibration-rotation coupling were neglected. 

The frequencies calculated in this way compare well with 

the observed frequencies. /3 In the calculations reported 
below, all transitions arising from the rotational levels 
with J = 0 to J = 21 have been included. In calculations 

of the v, spectral densities or reorientational correla- 


4 
tion functions, we calculated the function for each of 
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FIGURE-11 .""'J-pi Ff Fusion spectral band shapes I, (w) 
tor the V3 band of CH, at 295K for various 
values of Tz: Intensities are normalized so 


that the integrated intensity is constant for 


all values of Ty. 
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FIGURE 12. M-Diffusion spectral band shapes Ty (w) 
for the V3 band of CH, at 295K for various 
values of Intensities are normalized 
so that the integrated intensity is constant 


for all values of Ty: 
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FIGURE 13 
M-Diffusion reorientational correlation functions 


Geet for the V3 band of CH, at 295K for various 


values of Ts: 


-re{Gi") (t)}, ----- -Im{Gy") (t) }. For clarity, 
successive correlation functions have been displayed 
vertically by 0.2 units. Time is in reduced units 


as iniFig. 79. 
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the nuclear spin isomers and computed the composite func- 
tion as a weighted average with weightings A:E:T = 5:2:9. 
Calculations of a spectral densities and correlation 
functions based on average transition frequencies (the 
method employed in section III-D-l-a.and in section 
II-C-2 for the V3 band) gave results which were signifi- 
cantly different from those obtained in calculations 
based on all of the accurately known transition fre- 
quencies for each nuclear spin isomer as described above. 
These spectral densities are compared in Fig. 14 for 

Ts = 0.3, 0.5 and 1.3 picoseconds. 

The free rotor correlation function and memory func- 
tion for the V4 band of meta-CH, (A symmetry) at 295K 
are shown in Fig. 15. The corresponding functions for 
ortho- and para-CH, closely resemble those for meta-CH , 
at short times, but differ significantly at long times. 
This is not too surprising since the wings of the 
spectra for the three nuclear spin isomers are very 
similar, but the spectral intensities near the band cen- 
ters of eae teonee are quite different. Since the 
time-dependence of the correlation function at short 
times reflects the intensity in the wings of the 
- spectrum, and the long-time behaviour reflects the 
intensities near the band center, one would expect 


significant differences in the free rotor correlation 


and memory functions of the different nuclear spin 
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FIGURE 14 


J-Diffusion spectral band shape T, (w) for the vy 


eee omer erie a a 


band of CH, at 297K. 


All transition frequencies considered ; average 
transition frequencies considered:-----; T, as 


indicated. 
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FIGURE 15. Free rotor correlation and memory 


function for the V4 band of CH, at 297K. Time 


Li 
scale is in reduced units of t(kT/I)* = 


B.7 10k. 810) = 0.60. 
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isomers only at long times. 
The J-diffusion reorientational correlation functions 
ae te at 297K, for various values of the angular mom- 


entum correlation time Ts shown in Fig. 16 were calcu- 
lated from the correlation functions for each of the 
nuclear spin isomers, which were in turn calculated from 
the corresponding memory functions. The spectral densi- 
ties I,(w), at 295K, for the same values of T, are shown 
in Fig. 7. 

The spectral densities, Ty (#) + in the M-diffusion 
limit, for various values of the angular momentum cor- 
relation time T, were calculated from the spectral densi- 
ties for each nuclear spin isomer, which were calculated 
using Eqs. (III-46) - (III-48) and are shown in Fig. 18. 
The corresponding reorientational correlation functions, 
Go ake ys, were obtained by Fourier transformation of 
the spectral densities and are shown in Fig. 19. 


2. Comparison with Experiment 

In order to fit observed band contours of CH, an 
dense gas solutions e 4 with those calculated from the 
semiclassical M- and J-diffusion models, we have nor- 
“malized observed and calculated spectra so that they 
have the same integrated intensity (area normalization). 


The value of the single variable parameter of the 


models, t., was varied until the agreement between the 
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FIGURE 16 


J-Diffusion reorientational correlation functions 


Ge) for the V4 band of CH, at 297K for various 


valu £ B 
lues o Ts 


-re{c!") (t)}, ----- -1m{ce!?) (t)}. For clarity, 
successive correlation functions have been displaced 
vertically by 0.2 units. Time is in reduced units 
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FIGURE 17. J-Diffusion spectral band shapes I, (w) 


for the V4 band of CH, at 297K for various 


values of TS. Intensities are normalized so 


that the integrated intensity is constant for 


all values of Ts: 
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FIGURE 18. M-Diffusion spectral band shapes Ty fw) 


for the V4 band of CH, at 297K for various 


values of Ts: Intensities are normalized so 


that the integrated intensity is constant for 


all values of Ty: 
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FIGURE 19 


M-Diffusion reorientational correlation functions 


6!) (ty for the v. band of CH, at 29/K for various 


M 3 4 
values of T_. 
PEE ET OAS | 
—— Boicua. (t)}, ----- -Imicy”) (4) 3. “for 


clarity, successive correlation functions have 
been displaced vertically by 0.2 units. Time is 


in reduced units as in Fig. 15. 
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calculated and observed spectral densities was satis- 
factory as judged by visual comparison. Representative 
comparisons of calculated and observed V3 band contours 
are shown in Fig. 20, while representative comparisons 
for the Vy band contours are given in Fig. 21. The 
band shapes calculated with the M-diffusion model were 
incompatible with those observed in the low pressure 
CH ,/He mixtures because the intensity of the Q-branch 
was too low (relative to the P- and R-branches) in the 
calculated spectra for values of Ty where the broadening 
of the rotational structure in the P- and R-branches of 
calculated and observed spectra were comparable. At 
higher pressures, however, the band shapes calculated 
with the M-diffusion model, with suitable values of Tyr 
agreed with the observed band contours as well as the 
band shapes calculated with the J-diffusion model did. 
The spectra shown in Figs. 20 and 21 are compared with 
band shapes calculated with the J-diffusion model, since 
this model appears tobe more compatible with the 
observed spectra over the complete range of pressures. 
The calculated band shapes discussed above have 
been compared to experimental spectra assuming that 
' the experiments were performed under conditions in 
which instrumental effects on the band shapes were in- 
significant. As pointed out in Chapter II, the con- 
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Ty =0.20 psec 


CH4/He p=437 ; Ty =0.35 psec 


CH,/N2 p=186 Ty =0.40psec 


2900 3000 3100 3200 
w (cm!) 
FIGURE 20. Comparison of observed band shapes 

for the V3 band of CH, in dense gas mixtures 
at 295K and band shapes computed with the J- 
diffusion model. The intensities are 
normalized as in Fig. 18 and the indicated 
densities, p, are in amagats. 
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FIGURE 21. Comparison of observed band shapes 
for the V4 band of CH, in dense gas mixtures 


at 297K and band shapes computed with the J- 
diffusion model. The densities are in amagats 
and the band shapes are area normalized as in 


Fig. .18). 
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calculated band shapes produced significant effects on 
the band contours for the low pressure mixtures only. 
It may be that the apparent incompatibility of the M- 
diffusion model with low pressure CH, /He mixtures is 
due only to our neglect of instrumental contributions 
to the band shapes. Because the effects of the slit 
function on the band shape are of the same order of 
Magnitude as the effects of collisional broadening in 
this low pressure region, and because the exact nature 
of the slit function is not known, the angular momentum 
correlation time obtained in our J-diffusion analysis 
of these cases has a large uncertainty. 

The values of Ts which gave the best agreement 
between observed >? ana calculated spectral densities and 
reorientational correlation functions for the V3 and V4 


bands of CH, in a large number of fluid media are given 


4 
in Table 2. Since the quality of the fits of observed 
spectra with those calculated with the M-diffusion and 
with the J-diffusion model was not significantly dif- 
ferent except at low pressures, we have included the 
values of Ts which gave the best fits in each model. 
The values of Ts obtained by fitting the V3 data are 
“not significantly different from those obtained in 
Chapter II. 

It should be noted that the agreement between the 


calculated and observed V3 band contours of CH,» parti- 
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TABLE 2 


Correlation Times Assigned Using the Semiclassical J-diffusion (SCF 


and Semiclassical M-diffusion (SCM) Models.* 


System Band Temperature Density Ty psec Ty psec 
(°K) (amagats) (SCF) (SCM) 
CH, /He V3 295 96 143 b 
437 0.35 0.35 
544 0.30 0.30 
745 0.20 0.20 
vy 297 110 1.1 
134 1.0 
156 0.90 0.80 
178 0.80 0.80 
200 0.78 0.80 
208 0.70 0.75 
222 O13 0.78 
242 0.65 / Oso 
250 0.60 0.65 
290 0.55 0.60 
328 0.50 0.58 
365 0.48 0.53 
434 0.40 0.48 
491 0.35 0.45 
541 0.33 0.43 
588 0.30 0.38 
682 0.25 0.35 
687 0.28 0.35 
747 0.23 0,33 
CH,/N., V5 295 186 0.40 0.45 
243 0.33 0.35 


(continued....) 
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TABLE 2 (continued ) 
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System Band Temperature Density tT, psec psec 
(°K) (amagats) (SCF) (SCM) 

CH,/N, V5 295 308 0.28 0.28 
463 0.28 0.23 

a 297 82 0.65 0.70 

493 0.65 0.70 

104 0.60 0.65 

114 0.58 0.65 

138 0.50 0.58 

161 0.45 0.53 

183 0.40 0.45 

203 0.35 0.44 

222 0.33 0.38 

239 0.30 0.36 

276 0.30 0.35 

303 0.28 0.33 

CH, Ar V3 29 281 0.33 0.35 


2 Estimated error in value of Ty is + 0.05 psec. 


p No satisfactory fit with M-diffusion model. 
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cularly in the wings of the spectrum is better than the 
agreement between the observed V3 band contours and those 
calculated with the semiclassical model of Chapter II 
(see Fig. 5 and 20). This is due, in large part, to the 
inclusion of all lines arising from rotational manifolds 
with J = 0 to J = 25 in the calculations reported in the 
present work, while computational efficiency demanded 
that our previous calculations included only those tran- 
sitions arising from the J = 0 to J = 12 manifolds. 

We have also compared the correlation functions cal- 
culated by Fourier transforming the observed band con- 
tours oe with the correlation functions calculated from 
the M- and J-diffusion models. In making such compari- 
sons, one must be careful to choose the zero of frequency 
for the Fourier transformation of the spectral densities 
in a consistent manner, since the correlation functions 
are complex functions and the relative magnitudes of 
the real and imaginary parts will depend on the choice 
the zero of frequency. In numerical Fourier transforma- 
tion of experimental band shapes, it is most convenient 
to choose the frequency at maximum spectral intensity 
as the zero of frequency. However, in calculating the 
J-diffusion correlation functions using the memory 
eraneti on technique described in section III-C above, 
one requires that the correlation function be real and 


equal to unity at t = 0, and that the time-derivative 
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of the correlation function be zero at t = 0. This 
implies that the zero of frequency in the calculated 


spectrum is equal to the first moment, M(1), of the 


spectrum. It has been shown a that the true 
band origin Wo is given by 
w = M(1) - 2B, (III-63) 


Oo 


where B is the rotation constant of the molecule. This 
band origin lies closer to the frequency of maximum 
intensity than does M(1), therefore choosing Ww, as the 
zero of frequency is more nearly in keeping with choos- 
ing the frequency of maximum intensity as zero frequency 


7,10,11,25,38,39 WW. have 


which is the common practice. 
chosen therefore to compare the correlation functions 
obtained by Fourier transformation of observed band 
contours with Ww, as the zero of frequency with calcu- 
lated correlation functions rephased to this frequency 
origin. Since the M-diffusion correlation functions 
were computed by Fourier transformation of the calcu- 
lated spectral densities, the choice of the frequency 
origin was made in the same manner as that for the 


-correlation functions calculated from experimental 


spectral densities. 
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The comparison of calculated correlation functions 
and Fourier transforms of observed V3 band contours gave 
essentially the same values of tT, as comparisons of cal- 
culated and observed band shapes. Unfortunately the 


39 


points in the digitized v, spectra of CH, available 


4 
were too distantly spaced for computation of meaningful 
Fourier transforms. We have therefore not attempted to 
compare calculated and experimental correlation functions 
for the V4 band. 

E. DISCUSSION 

The semiclassical rotational diffusion models 

proposed in this chapter are similar to, but not identi- 
cal to, those proposed in Chapter II which were based on 


2,34 Tf we form the matrix 


Gordon's semiclassical theory. 
products indicated in Eq. (II-14) we obtain the reorien- 


tational correlation function in the form 


G(t) = Pa we exp[\,t] (III-64) 
7: 


where A is a diagonal matrix containing the eigenvalues 
of the matrix + iw. as defined in Eq. (II-16), the 


weights, Wee are given by 


7. 
W, = > ash > (S Di gFoede 
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(III-65) 
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and the columns of S are the eigenvectors of “1 + twos 
It is clear from Eq. (III-64) and (II-15) that the re- 
orientational correlation function is a superposition 
of exponential functions with complex arguments, and 
that the spectral density is a superposition of Lorentz 
functions centred at frequencies given by the imaginary 
part of Ay with widths given by the real part of tore 
The reorientational correlation functions in the 
present models appear, at the outset, to be significantly 
different in form from those of the previous models. 
However, if one employs the Laplace transform relation 
in Eq. (III-31) using the free rotor expressions Eqs. 
(III-38) or (III-46), one can easily show that the 


Laplace transforms of the correlation functions in the 


J- and M-diffusion limits are respectively 


p | 2 | 1 ¢ 
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where the indices k and £ in Eq. (III-66) represent the 
collection of labels {IK M/k, £m} in Eq. (III-38), and 
these indices in Eq. (III-67) represent the collection 


{Mk £,m } in Eq. (III-46); 


Pala semen el (III-68) 
(1) 
d, = (AraKM| D_, (2(0)]|£) , (III-69) 
and 
af ey 
PhiwaPlee & (III-70) 


The sums over k and £2 in Eq. (III-67) are restricted to 
transitions associated with a particular J-manifold. 
Since the Laplace transforms given in Eqs. (III-66) and 
(III-67) are of the form of ratios f(s)/g(s) of poly- 
nomials in the Laplace variable s, with the polynomial 
g(s) in the denominator being of higher degree than the 
polynomial f(s) in the numerator, the correlation func- 
tions themselves must be sums of exponential functions 
with arguments {At}, with the coefficients A, being 
the solutions to the characteristic equation g(s) = 0. 
Thus the correlation functions in our previous and 
present semiclassical representations of the M- and J- 
diffusion models can be written as superpositions of 
exponential functions of t. 


In order to investigate the similarities of the 
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two models in greater detail, we have compared the expon- 
ential parameters in the two approaches to the M-diffu- 
sion limit of the extended rotational diffusion model. It 
was found that both approaches gave identical exponential 
parameters if the ie matrix of the previous M-diffusion 
model were specified by 

(1-4?) /y i= j 

i J 
Oi; = : | | (III-71) 
~ds/ts 1 # J 


rather than by Eq. (II-13) as in Chapter II. Numerical 
comparisons of the weights associated with each of the 
exponential terms in the M-diffusion correlation functions 
in the two models indicated that, although the previous 
approach could be modified to give the same exponential 
parameters as the present model, the weights of the ex- 
ponential functions were different in the two approaches, 
and the differences between them varied significantly 
with TS: A comparison of the weights given by the two 
models, for several values of Tyr for the three expon- 
ential terms in the J = 1 manifold correlation function 
for the V3 band of CH, is given in Table 3. 

The identity of the characteristic polynomials, 
whose roots define the exponential parameters in the 
correlation functions, in the two versions of the semi- 


classical M-diffusion model for the case of three spectral 


lines led us to investigate the characteristic polynomi- 
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als in both M- and J-diffusion limits of the semiclassical 
rotational diffusional models when a large number of lines 
were included. The characteristic equation in the present 


J-diffusion model is 


k R#k 
(III-72) 
and can be rewritten in the determinantal form 
+ (s-i0,) ore fy. Sle ee i (oe 
at rk Je 1 4 Al J di. 
5 (s-iw, WO Te eae) 0 ame 0 
Poo! 2 J 2 
s(s-iw, ) 0 (eet ) e (@> %6 0 
P4334 3 J 3 
P ats ee ) 0 0 ws (st tu ) 
nnon n n 
= 0 
(1?1-73) 


Since the value of a determinant is unchanged when a 
linear combination of any number of columns is added to 


: 80 F ‘ 
any given column, one can rearrange the determinant in 


: Fe ee 2 
Eq. (III-73) by adding -P55d, times column 2, ~P334, 
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times column 3, ... and -P a,” times column n to col- 


umn 1 to obtain 


ee Ay: -1. -1 -1 
s+(1 Pit) /T to, —(stt, ~t.u),) —(stt, -i0,) mets —(stt tu, ) 
2 -l1 . 

~P,od,/Ty (stt -tu)5) 0 eee 0 
“p..a2/ -1. 

33°3 Ty 0 (s+t, “tW4) ee 0 
2p daft 0 ene: 

natgtes 0 iy Pa (stt ~w) 

= 0 (III-74 


in’ the derivation“o£ fq. (111-74) from (LLI=-73) itchas 


been noted that 


(ITII-75) 
which is implicit in Eqs. (III-66) and (III-67) also. 
The determinant in Eq. (III-74) can be modified by 


adding the first column to each of the other columns 


to obtain 
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2 F 2 ‘ 9 
s + (1 Pid /ty ww, ~P94/T; eae Pi 49;/ 5 
op d“/t a (1-P. dy / 1 tue Ce dor 
a0 ea 22 2) Syne ‘ 1 ay a 
2 2 2 : 
P antn! ty P an’n’ “3 eo e e st+(1-P id) [tT 7=% aS 
= 0 


(ITTE=76) 


which can be written in the form 


Js +0-ea] = [0+ a) - sil =0 


ai 
(III-77) 

where 1 is a unit matrix, Wo is a diagonal matrix con- 

taining the frequencies of the transitions, and the 


matrix [J] has elements given by 


~ 


2 F=~< 
(1-P, ,4,)/T, 5 ie ee 

Lr gin 2 ue 
-P,4d5/T; att tag 


(III-78) 


The latter form of Eq. (III-77) is precisely the char- 
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acteristic equation whose roots are determined in the diag- 
onalization of the matrix (- I] + iw,)- Thus, provided we 
describe the matrix elements of I] by Eq. (III-78), rather 
than by Eq. (II-12) of the previous chapter, the charact- 
eristic equations of the present model [Eq. (III-72) ] 

and the previous model [Eq. (III-77)] are identical. 

This implies that the two approaches yield correlation 
functions which are linear combinations of exponential 
functions of time [Eq. (III-64)] with the exponential 
parameters of the two approaches being identical. As 
pointed out above, the weighting factors of the exponen- 
tial functions in the correlation function are not identi- 
cal. The weights Ws in Eq. (III-64) are given by Eq. 
(III-65) for the matrix diagonalization technique employed 
previously, and by 


' Zz Gg Lamy Pyreceg) 
Poedte A OA, tLe) ae g, 


' 
ae g'#R 


4 
A CQ; - AO (TTI=79} 


k#i 

with the present approach. The weighting factors in 
Eq. (III-79) reduce to those given in Eq. (III-65) only 
for large values of TS where the off-diagonal elements 
of I are negligible. It should be noted that for the 
case when all lines have identical intrinsic line ampli- 
tudes d, = 1/3. , the model proposed in Chapter II will 
have the same [J -matrix as the model given in the pre- 
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sent work. In the classical limit, (large rotational 
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quantum numbers) this case is indeed realized. 

The major differences between the present theory 
and that given in Chapter II is the method of describing 
the "collision" process. In this chapter, we have viewed 
"collisions" as events which result in transfer of 
molecules between the available rotational states, while 
in the previous theory, "collisions" were viewed as events 
wherein line amplitude was transferred between spectral 
lines. This concept of transfer of line amplitude 
between lines has some relevance in a classical interpre- 
tation of rotational transitions ehiee but its relation 
to transfer of molecules between quantum states is ob- 
scured by the fact that molecules in a given state can 
undergo more than a single transition within the rotation- 
vibration band of interest. The approach in Chapter II 
is based on the intuitive assumption that molecules with- 
in a given quantum state could be subdivided so that 
equal numbers of molecules were associated with each of 
the transitions originating in that quantum state. The 
present theory does not require such assumptions. 


Sac" has pointed out that the line amplitude 


Gordon 
transfer approach is expected to describe the time de- 

pendence of the correlation functions at long times only. 
However, with the semiclassical approach presented here, 


no such limitation is apparent in the mathematical de- 


velopment of the model. 
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The most important improvement in the semiclassi- 
cal rotational diffusion models presented here is the 
simplification of the mathematical description of the 
correlation functions and spectral densities. The treat- 
ment in Chapter II was of limited applicability due to 
its reliance on the diagonalization of the complex non- 
Hermitian matrix, - Tl + wo, whose order is equal to the 
number of spectral lines included in the calculations. 
This diagonalization procedure becomes prohibitively 
expensive when a very large number of spectral lines 
(>40) are included. In the application of the theory 


in Chapter II to the v., band of CH computation costs 


3 4’ 
required us to limit our consideration to the lines 
arising from the J = 0 to J = 12 rotational manifolds, 
which means that ca. 4% of the methane molecules which 
occupy higher rotational levels were completely neglected. 
In spite of this deficiency, which led to significant 
discrepancies between calculated and observed spectral 
densities in the wings of the band, we were able to 
obtain satisfactory agreement between theory and experi- 
ment. The application of the previous theory to the V4 
band of CH, was impractical since one must include a 
much larger number of transitions in these calculations. 
With the present theory, no matrix diagonalization is 
required since the results are derived from a free rotor 


memory function which is in turn derived from a free 


rotor correlation function, or from the Laplace transform 
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of the free rotor correlation function. We have therefore 
been able to apply the models to the calculation of the 


V4 band of CH, where 363 transitions are included for 


ortho-CH 120 for para-CH, and 121 for meta-CH 


4’ 4° 


The angular momentum correlation times, obtained 


Tse 
from comparisons of calculated J-diffusion and observed 


v., and v, band shapes of CH, in dense gas mixtures, given 


3 4 4 
in Table 2, are plotted against gas density in Fig. 22. 


The values obtained from analysis of the V3 and Vy bands 


do not differ significantly, and as was the case in Chapter 
II (see Fig. 6), the values are comparable to those 
obtained by analysis of nuclear relaxation times in high 


56 


pressure gaseous methane. It is felt therefore that 


the parameter tT, is indeed a measure of the time interval 


J 
during which the rotational state of the molecules remains 
fixed. The 20 vs. p plots in Fig. 22 appear to have 


non-zero intercepts. This can be attributed to our 
neglect of instrumental broadening since the values of Ts 
which we obtained by comparing calculated and observed 
band shapes would be lower than the "true" Ty values in 
the low density mixtures where instrumental effects could 
contribute significantly to the band shape. 

It has been found in the analysis of spectral band 


7, 22-28 and nuclear relaxation times, 28-33 that 


shapes 
the classical extended rotational diffusion models 


account for a very large number of observations and are of 
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0 200 400 600 800 
p(amagats) 

FIGURE 22. Density dependence of jatar oO- CH,/He 
at 295K from Va On CH, /He at 297K from Var 
A - CH,/N, at 295K from Var + - CH, /He at 
297K from Var All values of T, are taken 
from comparison of observed band shapes with 
those calculated with the semiclassical J- 


diffusion model. 
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value in interrelating data from different kinds of 
measurements. The model we have developed inthis chapter 
is a semiclassical version of the extended-rotational 
diffusion model. This model can be applied efficiently 
to the study of band shapes composed of a large number 

of transitions such as the V4 mode of methane and should 
prove useful in the analysis of spectra where the quantum 
mechanical effects are large and rotational-vibrational 


coupling is important. The classical model would be 


invalid and of little value under these conditions. 
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THE KIVELSON LINEWIDTH THEORY 
CHAPTER IV 
A. INTRODUCTION 


In Chapters II and III, we have developed models 
for the description of reorientational motion in systems 
in which the angular momentum correlation time Ts and 
the reorientational correlation times Ty are of the 
same order of magnitude. Under these conditions, it 
is necessary to follow, in detail, the reorientational 
motion in individual diffusive steps in order to des- 
cribe the reorientation of a molecule in terms of a 
rotational diffusion process. In this chapter re- 
orientational motion in liguid systems where Ty is much 


shorter than Tp) will be studied. Under these condi- 
tions, molecules reorient by very short rotational 
steps and they sample a large number of rotational 
states in the period of time Tp: For this reason it 

is not necessary to follow the rotational process in 
each individual step and the reorientational motion can 
be described using a rotational diffusion equation. 


In the limit where t. is much shorter than Tos the 


J 

extended J-diffusion model reduces to the model des- 
P : 8 

cribed by the diffusion equation approach. 


If W[Q_;t,2] defines the conditional probability 
re) 


that a molecule has reoriented from an initial orienta- 
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tional 2, to.an orientation 2 in time t, the classical 


rotation diffusion theory 52 


leads to a description of 
the time dependence of WIL Ft, 2] in the form of a rota- 


tional diffusion equation 


3 2 : = 
- =¢ WIQ 7 t,21 -) ToPapZa (Boi 2, (IV-2) 


a,B 


where Jy is the dimensionless operator for the angular 
momentum about the a-th axis and meas is the a8-th com- 
ponent of the diffusion tensor D given by 


’ (iVri2) 


D see bi el. dée € p(e,At) 
= = a 8B 
ap At+0 2 a 


At 

where At is a time increment which is much shorter than 
the times of interest Tos and p(e,At) gives the prob- 
ability that the molecule reorients through a solid 
angle € in time At. In order to obtain Eq. (IV-1), one 
applies the rotation operator exp[-ie.J], representing 
an infinitesimal reorientation €, to the function 
WIQ Ft, 2], expresses the result as a power series in 
é and neglects all terms of order greater than a4 

The solutions to the rotational diffusion equation 
(IV-1) for a molecule with an isotropic diffusion 


tensor, Dl, can bewritten in the form 
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(3%) (2):* 
W[Q Ft, 2] -) ore » [20] rach [Q] 
T km 
&,k,m 


x- exp [-2(%4+1) Dt] (iV=3) 


68 which des- 


where eo) is the Wigner D -matrix 
cribes the reorientation of a spherical tensor of rank 
& through the Euler angles Q,. The diffusion tensor D1 
is usually defined in a coordinate system in which it 
is diagonal and is generally interpreted in terms of 
the retarding torques acting on the molecule due to 
the viscous drag of the solvent using the Stokes- 


a3 


Einstein hydrodynamic relationship. Debye first 


proposed this type of relationship in his classic 
treatment : of dielectric relaxation in polar liquid 
and it was later used by Bloembergen, Purcell and 


Pound 20 


to interpret the nuclear spin relaxation times 
of the molecules in liquid water. 

In addition to the study of Raman and infrared 
band shapes described in Chapters II and III, and the 
study of dielectric and nuclear relaxation phenomena 
as noted above, one can also gather information about 
molecular reorientation in liguids by studying electron 
spin relaxation of paramagnetic species in dilute 


liquid solutions. The interpretation of electron spin 


relaxation times is a particularly useful technique 
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provided that a suitable paramagnetic probe is available. 

If the electronic Zeeman and nuclear-electronic 
hyperfine interactions are anisotropic then the tumbling 
motion of the molecule modulates the strength of the inter- 
actions between the electron spin, nuclear spin and 
applied magnetic field and gives rise therefore to a time 
dependent perturbation on the electron spin levels which 
leads to transitions among them or spin relaxation. For 
dilute solutions of paramagnetic species, containing a 
Single unpaired electron and one or more magnetic nuclei, 
the transverse spin relaxation time, To, is determined 
by this rotational modulation of the spin levels. Since 
TS is inversely proportional to the observed spectral 
linewidths, information about the reorientational motions 
of the paramagnetic species is easily obtained from ex- 
periment provided a suitable description of the effect 
of the reorientational motion on the anisotropic inter- 
actions is available. 

The anisotropic electronic Zeeman and electronic- 
nuclear hyperfine interactions are characterized by 
second rank spherical tensors, so that the reorienta- 
tional correlation time which describes the time cor- 
relation of these tensors is the correlation time wget) 


associated with the correlation function 


*x 
(99) 2°) ial) . The rotational diffusion 
km m 


equation theory gives the result 
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(2)* (2) (2) * 
( D [21D (21) = an, fand [20] 
km km km 


(2) (2) 
Z 

iW ise 0} = 
xD [2] W [2 7€,9] (1 t9511”) 


x exp-6Dt (IV-4) 


so that (2) 


9 = 1/6D ° (IV-5) 


v 
19 
and Bloembergen, Purcell and 


Pound ss approach, we assume the Stokes-Einstein a re- 


Following the Debye 


lationship 


Bae kT/8nr > n (Iv-6) 


where k is the Boltzmann constant, T the temperature, ry 
the hydrodynamic radius of the molecule and n the co- 


efficient of shear viscosity of the solvent: so that 


pe z amr. °n/3kT (Iv-7) 


The ESR spectrum of a paramagnetic species with 
one unpaired electron and a single magnetic nucleus of 
spin I in a dilute liquid solution consists of 2I + 1 
lines each with a different linewidth. The widths of 
these lines AHy! depend on the nuclear spin quantum 


number, M, and can be expressed by the cubic function 


z 3 


AHy =a’ +a" + BM + yM + OM (IV-8) 
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where the parameters a', a", 8, y, and 6 depend on the 


species, temperature and solvent. Kivelson and Wilson 4+ 


have obtained theoretical expressions for the contribu- 
tion to the parameters a‘, 8, y, and 6 from the relaxa- 
tion process due to the rotational modulation of the 
electronic Zeeman and hyperfine interaction using the 


81,82 


Kubo-Tomita formalism for relaxation. Other workers 


eae have verified these expressions by deriving them 


with Redfield's theory of magnetic relaxation. 2°" 87 
The linewidth parameters are related to the magnetic 
parameters of the paramagnetic species and the electron 


resonance frequency by 
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m Q 4 _=ay/}8 2 
8 = —— 3 B, (bAy + AcdSy) ong (AyB.) 
573g 8 = 
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+ ame . ort + v}) 
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a 2 2 
+ a( Btbar + 4c5y) - ae Boas (1 + £) 
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(IV-11) 
and 
fas aa 2 
6 = b (a/w,) [1 + u(l + £)] (IV-12) 
10¥3g, 8. 
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By 
A = 2 <i ay ’ IvV-14 
Y Ra ( J. o. a.) ( ) 
§ Po ( ) 15 
a al ’ IV- 
Y I, ~ Fy ( ) 
Piet AD ks (IvV-16) 
is lime 6 y z " 
CLAS eee ey (Iv-17) 
3 Zz x y f 
c = =n - A.) (Iv-18) 
4°"X yen 


’ (IvV-19) 
and 
sig 2 2 
f = W, Te, U (IV-20) 
Tye Sy! op and AL, ther A. are the principal values of 


the electronic Zeeman and hyperfine interaction tensors 
respectively, Bo is the Bohr magneton, Wo the frequency 

of the applied radiation and By is the applied magnetic 
field. All of these parameters, with the exception of 
oie can be determined from the ESR spectra of the 
species in liquid and solid solutions. Therefore, a 
study of the linewidth contributions due to motional 
modulation of anisotropic interactions can yield measure- 
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The residual linewidth, a", found in Eq. (IV-8) 
can be determined by subtracting the a' term predicted 


by the theory of Kivelson and Wilson oe from the ex- 


(2) 
v) 


a careful investigation of possible relaxation mechan- 
43 


perimentally determined a' + a" if T is known. From 


isms, Atkins and Kivelson concluded that the a" 
contribution to the linewidth was dominated by motional 
modulation of the spin-rotational coupling. A theoreti- 
cal expression for a" was obtained by a procedure anal- 
gous to that used to obtain a, 8, y, 6. The magnitude 
of the spin-rotational contribution to the linewidth 

is determined by the magnitude of the components of the 
spin-rotational coupling tensor and the angular momentum 
correlation time ae Since the elements of the spin- 


rotational interaction tensor are not usually available, 


42,43 have used second order 


Kivelson and co-workers 
perturbation theory to estimate these elements from the 
elements of the g-tensor. Furthermore, they have used 


Hubbard's relationship 


ey i5 ae z 
53 la eh eA at ly (IV-21) 


where I is the moment of inertia of the molecule k is 
the Boltzmann constant, and t the temperature, to 
estimate Ts from ie With these assumptions, the spin- 


rotational contribution to the linewidth is 
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(IV-22) 
where 6g is defined by 


2 2 2,3 
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(IV=23) 
and Se is the free electron g-value (2.00231). 


a applied the linewidth theory 


Kivelson et al. 
described above to interpret the observed ESR line- 
widths of bis (2,4-pentanedionato) oxovanadium(IV) [VOAA] 
in a series of solvents. They obtained excellent agree- 
ment between theory and experiment provided that Eq. 
(IV-7) was modified to the form 


3 
ANKE | n 


uae 5 ph ere (IV-24) 
3kT 


here lo is the hydrodynamic radius of the solute mole- 
cule determined from translational diffusion experiments 
and k is a dimensionless parameter which they found to 
be temperature independent, but which varied with the 
nature of the solvent and solute molecules. After care- 
ful consideration of reorientation in a viscous fluid, 


42,90 


Kivelson et al. have defined k to be a tensor 


given by 
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where f¢ and FP, are the i-th components of the torque 
exerted on the molecule, and the intermolecular force 
experienced by the molecule respectively. So as not to 


introduce additional parameters into the model the tensor 


k is generally assumed to be a scalar.+/’+8-4? 


17,18,44,45 


Recent investigations of the ESR line- 


widths of VOAA and other vanadyl complexes in solutions 
at low temperatures and also in solutions in hydrogen 

bonding solvents have indicated apparent weaknesses in 
the Kivelson linewidth theory. Angerman and Jordan 34 


studied the linewidths of VOAA in a number of alcohols 


and noted that the reorientational correlation time, 


(2) 
"9 
different from that determined from the y parameter. 


, assigned from the analysis of the 8 parameter was 


They attributed this discrepancy to the high sensitivity 
of 8 to experimental error in the linewidths relative to 
the sensitivity of y, rather than to weaknesses in the 
theory itself. They also observed that the spin-rotation- 
al contribution, a", was much larger than that predicted 
by Eq. (IV-22) in the region of large n/T. Similar 
observations were made by Chasteen and Hanna a ina study 
of a series of vanadyl a-hydroxycarboxylate complexes in 
aqueous solutions, by Huang and Kivelson af for tetra- 
n-butylammonium bis (dimercaptomaleonitrile) nickel (ITI) 
[(n-C,H 9) .N.NiS,C, (CN) 4] in n-butanol and in CHCl,, and 


by Hwang, Kivelson and Plachy Ae for VOAA in a number of 
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non-hydrogen bonding solvents at low temperatures and at 
high pressures. 

In order to further investigate the validity of 
the Kivelson theory, we have studied a series of vanadyl 
complexes of varying molecular size in toluene, carbon 
disulfide, chloroform and diphenylmethane (DPM) solutions. 
The complexes studied were bis (3,5-heptanedionato) oxovan- 
adium(IV) [VODE], bis(1-phenyl-1,3-butanedionato) oxo- 
vanadium(IV) [VOPM], bis(2,2,6,6-tetramethyl-3,5-heptane- 
dionato) oxovanadium(IV) [VOTB], and bis (1,3-dipheny1l-1,3- 
propanedionato) oxovanadium (IV) [VODP]. 

In sections B-E of this chapter, the experimental 
methods employed in obtaining the data are described. 
In section F, the analysis of the data is presented 
while in section G the procedures geen to interpret the 


data are assessed and the implications of the results 


considered. 
B. PREPARATION AND CHARACTERIZATION OF VANADYL COM- 
PLEXES 


The four vanadyl complexes were prepared using 


variations of the standard procedures given in Inorganic 


91 


Synthesis for the preparation of VOAA. 


Ly Bis (3,5-heptanedionato) oxovanadium(IV) [VODE]. 


The VODE complex was prepared by refluxing a mix- 
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ture of 0.5 g vanadium pentoxide (V,0.) (Fisher Scient- 
ific) and 12.5 g of 3,5-heptanedione (Eastman Organic) 
for 1.5 hours. The hot reaction mixture was filtered to 
remove unreacted V0. and the excess diketone was re- 
moved under vacuum. The resulting product was recrys- 
tallized twice from acetone and Paeeed for microanaly- 
sis. 

Analysis: | calculated for Cy gH520,V: 

C7 52633 3 {Ny 16385 


Found: Gy o2.08 3) “He Gere 


aa Bis (2,2,6,6-tetramethy1l-3,5-heptanedionato) oxo- 
vanadium(IV) [VOTB] 


The VOTB complex was prepared by refluxing a mix- 


ture or 0so"9 V and 18 g of 2,2,6,6-tetramethyl- 


2°5 
heptane-3,5-dione (Pierce Chemical) for 4 hours. The 
hot reaction mixture was filtered to remove V,0. and 
the excess diketone was pumped off. The resulting 
product was recrystallized twice from acetone and sub- 
mitted for microanalysis. 

Analysis: calculated for CoH 0.V: 
Gy 02.0005 EH 2.280 


Found: Cprereue 7 © HH, 8.75 


Sea Bis (l-phenyl-1,3-butanedionato) oxovanadium(IV) [VOPM] 


The VOPM complex was prepared by refluxing a mix- 
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ture of 0.5 9 V,0. and of 1.0 g 1-phenyl-1,3-butane- 
dione (J. T. Baker) in 25 mls of benzene for 1 day. The 
hot reaction mixture was filtered and the benzene re- 
moved under vacuum. The resulting product was recrystal- 


lized three times from acetone and submitted for micro- 


analysis. 


Analysis: calculated for CoH gOKV: 
Cy, 61.70 ; H, 4.63 


Found: OL 2.67) seule Cae Oe 


4. Bis (1,3-dipheny1-1,3-propanedionato) oxovanadium (IV) 

[VODP ] 

The VODP complex was prepared by refluxing 0.5 g 
of V,0, and 0.5 g of 1,3-diphenyl-1,3-propanedione 
(Eastman Organic) in 5 mls benzene for 1 day. The hot 
reaction mixture was filtered and the benzene removed 


under vacuum. The resulting product was recrystallized 


twice from acetone and submitted for microanalysis. 


Analysis: calculated for C39 Ho 00. V: 
C,e7O0018s; Hy eCSV29 


Found: CFa70 . SOSF4H Pp A425 


PC. SAMPLE PREPARATION 
The solvents toluene (Allied Chemical reagent), 
carbon disulfide (Matheson, Coleman and Bell), chloro- 


form (Mallinckrodt Analytical reagent), and diphenyl 


net i 1 a0 oe 


~onntunlely E> teenie ts eos te ah 

ent web £, x03 sisened to simi e ab 
“as atenned ody Bes bowed 22 pew 
~Latayross. ew Poehose~ edt La, ae 
crate +6 bods Endep: one p 


aiey 


digit, (itou Bere estes 
inpodg@th ‘big 4aaeaes (id 


* 


vite coe 
VAL. 


156. 


Methane (Fisher Scientific) were all purified by vacuum 
distillation and were stored over molecular sieves (BDH 
type 4A) under vacuum. Solutions of the vanadyl com- 
plexes were prepared by weighing the complex into a 
vessel with a sidearm fitted with a 4 mm ESR Pyrex sample 
tube. The solvent was then distilled into the vessel 

On a grease free portion of a vacuum line. Any traces 
of oxygen remaining in the samples were removed by 
pumping on these solutions while they were maintained 

at -95°C in a toluene slush. The solutions were then 
poured into the sample tubes and sealed off under vacuum. 


All solutions were o7 M. 


De INSTRUMENTATION AND MEASUREMENT TECHNIQUES 

The ESR spectra were recorded on a Varian V-4502 
spectrometer equipped with a Varian V-4557 temperature 
control accessory, and the temperature was measured 
with a copper-constantan thermocouple The magnetic 
field was calibrated with an Alpha model 3093 digital 
NMR gaussmeter. The cavity frequency was determined 
with a Hewlett Packard model X532B frequency meter. 
The frequency meter was found to be accurate to 0.05% 
using an aqueous solution of Fremy's salt (g = 
2.00550) .°7 

The widths of the individual spectral lines were 


determined by measuring the heights of the lines on a 
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1 kG sweep and the width of the M = -3/2 line ona 
narrower sweep (25 - 50 G). The lines were assumed to 
be of Lorentzian shape so that the height, H, and the 
width, W, were related through the familiar relationship, 
C= HW, where C is constant for a given spectrum and 
can be determined from the height and width of the M = 
-3/2 line. Preliminary investigation indicated that 
the influence of line overlap on the measured line- 
widths was negligible even at low temperatures and 
thus no correction was made for overlap. In order to 
decrease the experimental error in the height the 
average of four sweeps was used in each calculation. 
Tables of the observed spectral linewidths are given 


in Appendix III for each complex/solvent system. 


E. ANALYSIS OF THE ESR SPECTRA OF VANADYL COMPLEXES 
The spin Hamiltonian,Jj, for a paramagnetic mole- 
cule containing a single unpaired electron and a single 


magnetic nucleus is 


+> 


H = 6 8-9-8, + fi I.A. (Iv-26) 


where Ss is the electron spin angular momentum operator, 
“ft the nuclear spin angular momentum operator and By 
the applied magnetic field. g is the electronic Zeeman 
interaction tensor and A is the nuclear-electronic 


hyperfine interaction tensor. Since the magnetic field 
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vector Bo is fixed in the laboratory coordinate system 
and the tensors g and A are diagonal in a molecular co- 
ordinate system, the interactions are dependent on the 
orientation of the molecular coordinate system relative 
to the laboratory coordinate system and thus the ESR 
transition frequencies are orientation dependent. If 
we define Os ey WO, respectively as the ESR transition 
frequencies when the x, y or z axis of the molecular 


coordinate system lies along the direction of the 


applied magnetic field then these frequencies are given 
me 


by 
2 2 
g_8 B NC (EE +e.) =e} 
Te es Aaa al x 
TA 4 GoPoa 
(IV-27) 
where 
2% 2 A nee AK (IV-28) 
g x oe e q 


and gq = x, y or z. Ina solid solution the molecules 
are randomly oriented, however the position of the 
individual transitions can be identified and these 
transition frequencies can be used in conjunction with 
Eq. (IV-27) to determine the individual components of 
the g and A tensors. In Fig. 23, the spectrum of VOTP 
in a toluene glass is given and some of the transition 


lines identified. 


In a liquid solution where the molecule is tumbl- 
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ing freely, the x, y, z components of each spectral 


line are averaged to give lines at the frequencies 


8 B Aa 
oo = Sa + aM + 5 ae ee (IV-29) 
where By is the Bohr magneton, B the magnetic field and 
So and a have been defined in Eqs. (IV-13) and (IV-16) 
respectively. The liquid spectrum consists of 2I + 1 
lines whose positions and widths are determined by the 
averaging produced by the tumbling motion of the mole- 
cule. It is this averaging process which is monitored 
in studying the reorientational motion of molecules 
using the linewidth of the ESR spectra. From an itera- 
tive procedure. using Eq. (IV-29) and the position of 
the 8 lines of the liquid spectrum, So and a can be 
calculated. Although the g and A tensors could be 
determined from the glass spectra, Wilson and Kivelson'? 
suggest that determining qe and a from the liquid 


spectrum, Go g Sur Aw and A, = A, from the glass 


x “y 

spectrum and combining these data to obtain the re- 

maining parameters is the best approach. g, can be 

determined from the positions of the + M, lines in the 
- i - A. - A_ can be determined 
z-spectra while g,. dy and " y ine 


from the separation of the x- and y-spectra of a given 


M, transition (see Figs. 23,24). 
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In Fig. 24 a section of the VODE/toluene glass 


spectrum, in the region of the M_ = 3/2 x and y tran- 


I 
sitions, is shown and the positions of the x and y 
transitions are indicated. 

The components of the g and hyperfine interaction 
tensors determined from the glass and liquid spectra 
using the procedures outlined above are given in Table 
4. These determinations were all carried outon tolu- 
ene solutions of the complexes. 


In Fig. 25, the region of the M_ = 3/2 x and y 


zx 
transitions in the VOPM/toluene glass spectrum is 
shown. Comparing this to the same region of the VODE 
spectrum shown in Fig. 24, we see that there appears 

to be an additional transition in the region where 

the y-transition is expected. We have analysed this 
spectrum in terms of the two possible y-transitions 

and obtained two sets of parameters. VOPM can exist 

as either the cis or trans isomer. X-ray diffraction 
studies indicate that the cis isomer is the predominant 
species in the solid State." However ESR studies by 
Belford et ale on single crystals of bis(l-phenyl- 
1,3-butanedionato) palladium which has a trans geometry 
_doped with VOPM indicate the possible existence of 


VOPM as the trans isomer, while the VOPM complex in 


bis (1-phenyl-1,3-butanedionato) zinc crystals (cis 
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Magnetic Parameters of Vanadyl Complexes in Toluene 


Parameter 


Estimated error: 


TABLE 4 
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10 Gauss 


3380G 


FIGURE 25. ESR Spectrum Of VOPM in toluene art r= 90°K 


in the region of M = 3/2 x, y-transition. 
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geometry) exists as the cis isomer. The g values 
assigned by Belford etic-valbeme to theccis:(Zn matrix) 
and trans (Pd matrix) isomer of VOPM are similar to 
those assigned from our analysis of the glass spectra 
indicating the simultaneous existence of both the cis 
and trans isomers in the toluene glass. The values 
reported in Table 4 for VOPM are an average of our two 
sets of results and are sufficiently accurate for the 


purpose of our linewidth analysis. 


B. ANALYSIS OF ESR LINEWIDTHS 

The theoretical expressions for the linewidth 
parameters [Eqs. (IV-9) - (IV-12)] can be combined 
with the expression for «fle [Eq. (IV-24)] to obtain 
expressions for a', 8, y and 6 in terms of the ratio 
n/T and the adjustable parameter x. Provided n/T is 
known, K can be determined from a nonlinear least 
squares fit of the theoretical linewidth parameters to 
those determined from experimental linewidths. 

In previous applications of the Kivelson linewidth 
theory, 17,18,41,43-45,89 | Jas determined from the n/T 
dependence of the y parameter. The procedure had cer- 
tain advantages, y is nearly a linear function of n/T, 


and kK can be determined directly from the slope of 


a plot of y vs n/T. The values of y are usually more 
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precise than the values of 8 since 8 is determined 
essentially from the differences in observed linewidths. 
Consideration of y is also more pleasing in that the 
agreement between theory and experiment is more satis- 
factory than for 8. On the other hand, consideration 
of only the y parameter is not a stringent test of the 
complete theory since only a portion of the theory is 
being utilized. Secondly, the deviations between 
experimental and predicted 8's cannot be entirely 
attributed to experimental uncertainties. Finally, 
fitting one's data for only the parameter which gives 
the best agreement is a questionable procedure especi- 
ally in the case of y, which is less sensitive to changes 
in n/T than is 8. It is for these reasons that we have 
adopted the procedure of analysing the M-dependent parts 
of linewidths. We defined the M-dependent part of the 
linewidths as the experimental linewidths minus the 

a' + a" contribution obtained by a least squares fit 

of the experimental linewidths to Eq. (IV-8). It would 
be more satisfactory to fit the entire linewidth using 
all the linewidth parameters, however the theory of 

the residual linewidths is still open to question in 
view of recently reported results in regions of large 


18,44 


n/T A nonlinear least squares routine was used 


to fit the M-dependent part of the linewidth, obtained 
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at several different temperatures, to the theoretical 
expression for 8M + ym? + 6M? employing Eqs. (IV-10)- 
(IV-12) and (IV-24). Portran IV listings of the com- 
puter routines used in both the linear and nonlinear 
least squares analysis are given in Appendix IID. The 
values of the product si obtained from the analyses of 
the M-dependent part of the linewidth are given in 

Table 5. In Table 5, we have also included the values 
of ce obtained from a nonlinear least squares fit of 
the n/T dependence of experimental 8 and y's [deter- 
mined by fitting the experimental linewidths to Eq. (IV- 
8)] to Eq. (IV-10) or Eq. (IV-11). A similar analysis 
was not carried out on 6 because of the large experi- 
mental errors in this parameter. A study was also made 
in which the experimental 8 and y's were simultaneously 
Fitted to Eq. (IV=-10) and Eq. (IV-11). As expected, 
this latter approach yielded values of Kee between 
those obtained from the separate 8 and y fits and are 
not included in Table 5. In these nonlinear least 
squares analyses convergence was determined by mini-~ 
mizing therelative residuals for the 8, yand 8£y fits 
and by minimizing the absolute residual for the fit 

_to the M-dependent part of the linewidths. The con- 
vergence of the linear least squares routine in fitting 


Eq. (IV-8) was also determined by minimizing absolute 
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residuals. In these least squares analyses several 
methods of weighting the data based on widths of the 
spectral lines were examined however, these different 
weighting procedures had an insignificant effect on 
the final result. In general, the results obtained by 
analysis of the y vs n/T data are similar to those 
obtained from the analyses of the M-dependent part of 
the linewidths. The VODP/DPM system was not analysed 
because the very broad lines in these spectra made the 
measurement of the spectra difficult since overlap 
between successive lines was significant. 

In Figs. 26-33 the 8 and y parameters calculated 
from Eqs. (IV-10) and (IV-11) using the values of Kee 
obtained by fitting the M-dependent part of the line- 
widths are compared with the experimental values for 
the complexes in toluene and carbon disulfide. It 
should be noted that at high temperatures these com- 
plexes begin to decompose and that the scatter of 
points in this region is possibly attributable to this 
fact. No effort was made to overcome this problem 
or to try to obtain data at higher temperatures be- 
cause viscosity data was not available in the region 
above the standard boiling points of these liquids. 
Data were not obtained in the very low temperature 
region because of the lack of viscosity data. The vis- 


cosities were calculated from standard foxenelae. > ° 
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8, y (GAUSS) 


0 ] 3 5 
n/T x 10° P/°K 


FIGURE 26. 8 and y vs n/T for VODE/toluene. 
Experimental, + - B, X - Y- Predicted 
by Kivelson theory, - with Stokes- 
Einstein, O - without Stokes-Einstein. 
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n/t = 10° P/°K 


B and y vs n/T for VODE/CS,. 


FIGURE 27. 
Predicted 


Experimental, + - By * — Y- 
- with Stokes- 


by Kivelson theory, 
Einstein, O - without Stokes-Einstein. 
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8, y (GAUSS) 


0 1 3 


n/t 10° P/°K 


FIGURE 28. 8 and y vs n/T for VOPM/toluene. 
Experimental, + - B, X - ¥. Predicted 

- with Stokes- 

Einstein, O - without Stokes-Einstein. 


by Kivelson theory, 


Stl 


fo 
Or 
f 


——— 


tay ® 


i | | 
eee Bomar here 0's Seka 
; = 


ar\@ Bia Ait. 
Ae aot som ios 


~seite ; he as Xe 


) castes 


173. 


8, Y (GAUSS) 


1 Z 
n/t x he P/°K 


FIGURE 29. 8 and y vs n/T for VOPM/CS,,. 
Experimental, + - 6, * = y. Predicted 
by Kivelson theory, - with Stokes- 
Einstein, O - without Stokes-Einstein. 
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8, Y (GAUSS) 


0 1 3 


nvr 10> Py -k 
FIGURE 30. 8 and y vs n/T for VOTB/toluene. 
Experimental, +°- By," - Y« Predicted 
by Kivelson theory, = with Stokes- 


Einstein, O - without Stokes-Einstein. 
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8, (GAUSS) 


0 1 a ae | 3 
n/t x 10° P/°K 


ELCURE, 32. 6 and ¥, WS nZT yior VOTB/CS,- 
Experimental, + - By x - Y- Predicted 
- with Stokes- 


Einstein, O - without Stokes-Einstein. 
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8, y (GAUSS) 


0 1 z--—L 5 
Tie 10° P/°K 


FIGURE 32. 8 and y vs n/T for VODP/toluene. 
Experimental, + - 8, x - Y. Predicted 

- with Stokes- 

Einstein, O - without Stokes-Einstein. 


by Kivelson theory, 
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Bry (GAUSS) 


0 re 2 | 3 
n/t x 10> P/°K 
FIGURE 33. g and y vs n/T for VODP/CS,. 


Experimental, + = 8, X - y. Predicted 
- with Stokes- 


by Kivelson theory, 
Einstein, O - without Stokes-Einstein. 
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The region of temperature and n/T over which the vis- 
cosity formulae are valid are given in Table 6. 

From a consideration of Figs. 26-33, it is clear 
that, in general, the agreement between theory and ex- 
periment is very good at high temperatures and that it 
becomes very poor at low temperatures. From examination 
of Figs. 26-33 and similar figures for the chloroform °? 
and DPM solvent systems and the VOAA results of Kivel- 


18,41 we conclude that in order of decreas- 


son et - al. 
ing "goodness" of fit the complexes are VOAA, VODE, 
VOPM, VODP, VOTB. With the exception of VODP and VOTB 
this order is also the order of increasing molecular 
weight and increasing molecular size. 

The analyses of the experimental data using the 
Kivelson theory yielded values for the product Kee 
Any assignment of the parameter K is dependent on the 
magnitude of the effective hydrodynamic radius, ro: 
The determination of Lo is somewhat uncertain, however 
reasonable estimates can be made from translational 
diffusion experiments interpreted using Stoke's Lawes 76 
Such experiments were not performed because we are 
primarily interested in the relative values of k fora 
given complex in different solvents. In order that 
our K's bear some resemblance to those already appear- 


ing in the literature oe we have adopted the following 


procedure. An approximate hydrodynamic radius re for 
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TABLE 6 


Sources of Solvent Viscosities 


Temperature range 


(°C) 
0 = 107 
0 - 46 
=. 5) 60 
60 - 260 


n/T range 
(107? P/eK) 
28 = HOT 
1.8 - 0.8 
3.2 = 1.2 
4,5 7=- 0.5 
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Source 
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Ref. 95 


Ref. 95 
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each complex/solvent system was calculated from the 
experimentally determined product ah for that system 
by assuming that the value of K was equal to the k for 
VOAA in that solvent. ?? An average hydrodynamic radius 
Lo for a given complex was obtained by taking the aver- 
age of the values of ro! for that complex. The values 
of K determined from the average hydrodynamic radii and 
the values of beat obtained from the analysis of the M- 
dependent part of the linewidths, along with the values 
for VOAA in the same solvents are reported in Table 7. 
We see that the k's for the different complexes in 

a given solvent are similar - a result which is indeed 
expected in view of the interpretation of K given by 
Eq. (IV-20) and the similar nature of the complexes 
studied. 

In the vanadyl systems studied here, the spin- 
rotational contribution to the linewidth, oq", cannot 
be calculated independently of the larger contribution 
a' from the modulation of the g and hyperfine inter- 
action tensors. An "experimental" q" is calculated by 
subtracting the value of q' predicted by the Kivelson 
theory from the q' + q" term obtained by fitting Eq. 
(IV-8) to the experimental linewidths. Unfortunately 
this procedure results in the inclusion of all the 
experimental error in q' + q" into the much smaller q" 


term. Secondly, any deficiency in q' predicted by 
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TABLE 7 


& 
Values of k Assigned to Complex-Solvent Systems 


Complex Toluene CHC. cs DPM 


* 
Values indicate best fit but probably significant to 


103%. 


kk 
From references 41,42,89. 
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the Kivelson theory is included in the experimental a" 
term. These two possible sources of error must be con- 
sidered when comparing experimental and predicted val- 
wes of a”. 

a" predicted by Eq. (IV-17) is compared with the 
experimental a" in Figs. 34-35 for the cS, and toluene 
systems. As was the case for 8 and y for these systems 
there is excellent agreement between predicted and ex- 
perimental values in the region of small n/T (high 
temperature) and the agreement becomes very poor in the 
region of large n/T (low temperature). We noted that the 
deviation between predicted and experimental values 
occurred at larger values of n/T for the small complexes 
and for solvents with small k's as was the case for the 
8 and y parameters. This consistency suggests the pos- 
sibility that part of the discrepancies between predicted 
and experimental values of a" may be due to the failure 
of the Kivelson theory to correctly predict a' at low 


temperature. We will have more to say about this 


point in the next section. 


G. DISCUSSION 


In the previous sections we have suggested that a 
fit of the M-dependent part of the linewidths to the 


corresponding theoretical expression is a more repre- 
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FIGURE 34 


a" vs n/T for toluene systems 


Experimental - +, predicted by Eq. (IV-22) - 


a" (GAUSS) 


E 
n/T x Ton P/ ik 


ta 


i= re ns v a y 


win 


> tae E : 


ea 


- ee 


a ar 


ta ave 
— —_ | 


+ xa ebHAiwaib aodzse 20% 


7 aa 


FIGURE 35 


a" vs n/T for carbon disulfide systems 


Experimental - +, predicted by Eq. (IV-22) - ———. 


a" (GAUSS) 
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sentative test of the theory than is the previous method 
of fitting only the parameter y. Previous workers have 
justified the neglect of the 8 parameter in terms of 

the sensitivity of 8 to experimental error. We have 
performed the following computer experiments in order to 
evaluate the sensitivity of each parameter to experiment- 
al error. Eq. (IV-8) was used to generate a set of 
linewidths from a given set of parameters. An error 

of + 10% was then added to these linewidths and these 
modified linewidths were fit to Eq. (IV-8). The sign 

of the error added to each linewidth was determined 


ey The firsti$ digits 


from a table of random numbers. 
of the table were associated with the 8 lines of the 
first spectrum, the next 8 with the lines of the 
second spectrum and so on. When a digit in the range 
0-4 was encountered, an increment of -10% was added 

to the width of the line associated with this random 
digit, while an increment of +10% was added if the 
digit encountered was in the range 5-9. The results 
of the least squares analyses of the modified line- 
widths are shown in Table 8. The three sets of para- 
meters selected are representative of the parameters 
in the low temperature region. Also included in Table 
8 are the results for systems in which the 10% incre- 
ments were all assigned the same sign and the cases 


where the sign of the error was determined by the sign 
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TABLE 8 (continued) 


Numbers in brackets indicate the percent contribution 


of the parameter to the sum defined by 


) Ie + g[M| + yM? + sin? 


M 


Numbers in the table represent the percent relative 


error in the parameter resulting from the application 


of the indicated random error. 
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of M. The system in which the sign of the error was 
assigned according to the sign of M gave relatively 
large errors in @ and relatively small errors in y com- 
pared to those systems in which the sign of the error 
was determined randomly. Although the results in Table 
8 represent only a few of the very large number of 
possible assignments of error it isclear from an exam- 
ination of the results that no binding generalization 
can be made about the significance of experimental 
error in 8 and y. Although y certainly appears to be 
less sensitive to error than 8, the errors in 8 are 
not sufficient to warrant ignoring the parameter com- 
pletely. It is interesting to note the apparent 
stability of a and also that parameters which make 
larger contributions to the linewidth appear to be less 
sensitive to the error than those that make smaller 
contributions. The stability of a is particularly grati- 
fying since the procedure of subtracting a from the 
linewidths to obtain the M-dependent component adds 

any error in a to the error in the M-dependent compon- 
ent. Thus it would appear that the procedure of fitting 
the M-dependent part of linewidths provides the best 
compromise between a 8 or y fit in that the importance 
'of each term is considered. 

In order to study the correctness of Kivelson and 


Wilson's expressions [Eqs. (IV-10) - (IV-12)] independ- 
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dently of Eq. (IV-24), we have fit the M-dependent 
parts of the linewidths from each spectrum to the cor- 
responding theoretical expressions to obtain a value 
(2) 
6 e 


fOr “1 This procedure is equivalent to allowing «Kk 


to be temperature dependent. The ee 


s.determined by 
this procedure are plotted against n/T in Fig. 36 for 
three of the complexes in carbon disulfide, also in- 


cluded in this figure are the values of eae 


pre- 
dicted by Eq. (IV-24) using the appropriate k from 
Table 7. These plots are very similar in appearance 
to the Y vs n/T plots given in Figs. 26-33. We see 
that Eq. (IV-24) is compatible with the observed re- 
mite at high temperatures, but, at lower temperatures, 
marked deviations occur. The fact that the values of 


a determined directly from a fit of the M-depen- 


3 
dent part of the linewidths for a given spectrum to 

the theoretical expressions differ from the value of 
rue obtained when Eq. (IV-24) is employed suggests 
that the requirement that the validity of Eq. (IV-24) 

be assumed throughout gives rise to a discrepancy bet- 
ween predicted and experimental values of the para- 
meters. The values of 8 and y obtained from these has 
have also been plotted in Figs. 26-33. From Figs, 26- 
“33 we see that the quality of the fit has been improved 


by neglecting the Stokes-Einstein relationship. This 
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FIGURE 36 


T vs /T for carbon disulfide systems 


Predicted by Kivelson theory with the Stoke- 


Einstein relations - , predicted by Kivelson 
theory without Stokes-Einstein relation, GJ -VODP, 
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observation could be accounted for in one of two ways: 
either Eq. (IV-24) is not correct over the whole range 
of n/T or that the temperature and viscosity data is in 
error. Since the deviations generally occur in the 
regions outside the region where the viscosity formulae 
are valid the viscosity data is certainly open to 
question, however errors in the viscosity data would 
not account for the dependence of this deviation on 

the complex itself. For example consider the toluene 
systems shown in Figs. 26,28,30,32. For the VOTB and 
VODP systems the theory appears to break down at n/T 
gope,derpd, ip Tg-* poise/°K while for the VODE and VOPM 
systems the theory appears to be valid as far out as 
Netra e 2 ee to° poise/°K. Similarly for the car- 
bon disulfide systems, the VODP and VOTB results are 
compatible with the theory up to n/T = 1.3 - 1.7 x ri ee 
poise/°K while the VODE and VOPM are good out to 

N/ie= Let =" 1. o x to™ poise/°K. Note that the upper 
limits for a good fit nearly coincide with the upper 
limit of the viscosity formulae (see Table 6), this 

was also the case for the chloroform and DPM systems. 
This is possibly pure coincidence since the VOAA/tolu- 
ene data already appearing in the literature na tee 
produce good fits in regions outside these limits. 


The values of eis determined by analysing the 
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data independent of Eq. (IV-24) were used with Eq. (IV- 

22) to predict new values of &". The overall agreement 
between these "'s and the experimental a"'s obtained 

7 (2) 
) 


using the a" predicted by Eq. IV-S and the new T is 


not significantly better than the agreement between pre- 


dicted and "experimental" values when Pe 


was re- 
stricted by the Stokes-Einstein relationship. 

Kivelson et arate have considered anisotropic 
rotational diffusion as a’ possible explanation of the 
discrepancy between predicted and experimental values 
of a". However, for complexes such as VOAA and those 
considered in this study the value of pee for re- 
orientation about the V-O axis should not differ by 
more than a factor of 3 or 4 re from that about axes 
perpendicular to the V-O axis, and, since the contri- 
butions to the linewidths from rotation about the V-O 
axis involve only quadratic and bilinear terms in the 
very small parameters Sy and c, it is very unlikely 
that anisotropic diffusion could account for the 
observed discrepancies. 

In this study we have attempted to pinpoint the 
source of discrepancy between experimental linewidths 
and those calculated from the Kivelson linewidth theory 
assuming a modified Stokes-Einstein relationship. Our 


results indicate that the Stokes-Einstein relationship 


is the main source of the discrepancy between the pre- 
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dicted and experimental parameters 8 and Y. However, 
this observation does not account for the failure of 
the theory to predict the spin-rotational contribution 
to the linewidth, a". Thus we conclude that for large 
values of n/T, the Stokes-Einstein relationship is not 
correct and that the simple description of the spin- 
rotational contributions to the linewidths given in 


Eq. (IV-20) is inadequate. 
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CHAPTER V 


Spectral band shapes contain a gréat deal of in- 
formation about molecular dynamics and intermolecular 


Snteractions.*'" ~° 


In order to extract this information 
from the band shape, one moe predict the influence of 
these factors on the observed spectral band shape. ii 
this thesis, we have considered two different sources 
of information about molecular reorientation. Chapters 
II and III have been concerned with the extraction of 
reorientational information from infrared band shapes, 
while Chapter IV has been concerned with the extraction 
of this information from electron spin resonance spectra. 
The Fourier transform of the infrared band shape 
is the reorientational correlation function of the 


2,4 


transition dipole moment of the molecule. Thais cor- 
relation function contains information about how the 
orientation of the molecule changes in time. The 
angular momentum of a molecule is changed by strong 
torques exerted on the molecule due to the strong 
intermolecular interactions of the colliding species, 
during a collision. In principle, if the. incermoteculer 
potential for the molecule in a fluid is known, then 

the average rotational motion of the molecules and 


thus the reorientational correlation function can be 


predicted. This predicted correlation function can then 
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be compared to the corresponding experimental correlation 
function and the accuracy of the intermolecular potential 
assessed. 24737 
The aim of the work considered in Chapters II and 
III was to introduce a model of reorientational motion 
which could predict the infrared band shape of a reorient- 
ing molecule when that band shape was distorted by quantum 
mechanical interactions. The previously existing theor- 
ies fue Orb 2030 were either only applicable to systems 
in which the band shapes were nearly symmetrical or were 
highly mathematical theories, which were not easily 
applied. We have developed two semiclassical models: the 
first a modification of Gordon's more sophisticated 


semiclassical theory ARTES the second a semiclassical 


extended rotational diffusion theory. 759 


Gordon's semiclassical band shape theory Se? 
required explicit knowledge of the intermolecular pot- 
ential of the molecule in order to predict the outcome 
of a collisional event using semiclassical scattering 
theory. In Chapter II, Gordon's semiclassical theory 
was modified by describing collisional events in terms 
of the M- and J-diffusion limits of extended rotational 
diffusion theory. This model was successfully used to 
interpret the V3 infrared band shapes of methane in 


gaseous and liquid mixtures. This model is practical 


only for band shapes which could be represented by a 
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small number of individual transitions (<40) and for this 
reason could not be applied to band shapes of larger 
molecules where a significant number of rotational states 
were populated, nor could it be applied to the V4 infra- 
red band of methane in which Coriolis coupling caused 
significant splitting of the transition lines. 

In Chapter III we developed a "true" semiclassical 
extended rotational diffusion model. This model is not 
severely restricted by the number of individual transi- 
tions which could be included and was applied to the Vy 
infrared band shape of methane in gaseous mixtures. 

The angular momentum correlation times obtained 


from the analyses of the V3 and V4 band shapes 32 and 


those obtained from an independent NMR-T) study 36 show 


a similar density dependence (see Figs. 6,12), which 
indicates that the values of Ts determined from these 


studies are indeed a measure of the time between strong 


30 


collisional events. Recently Maryott et al. have 


zs predicted from NMR-T, studies for several 


liquids including methane, with reat the time between 


compared T 


collisions predicted using either a cell like or gas like 
liguid model. They reported that the ratio ba ae 
varied between 1 and 2 for both models indicating that 

Ty is a reasonable measure of the time between collisions, 


not a meaningless adjustable parameter of the theory. 


In introducing the extended J-diffusion picture of a 
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collision. we have ignored the details of the inter- 
molecular potential in the liquid, however the agreement 
between experimental and predicted band shapes for the 
CO/He system obtained by Gordon and McGinnis a6) employ- 
ing an assumed intermolecular potential in their molecu- 
lar dynamics calculations, is no better than ours using 
the simpler J-diffusion picture. Thus we must conclude 
that until more accurate intermolecular potentials become 
available our simple theory is adequate for the interpre- 
tation of band shapes in terms of intermolecular inter- 
actions. 

Because we have relied on the picture of a collision 
proposed by classical extended rotational diffusion in 
Chapters II and III, it is relevant at this time to very 
briefly mention some of the experimental evidence which 
supports the J-diffusion model as a reasonable interpre- 
tation of reorientational motions. Campbell, Seymour 
and Jonas Fe have investigated the relationship between 
T gute and Ts predicted by the classical J-diffusion 
model. They studied the CF, in the gas phase (density 
range 9-120 amagats) and determined Ts from ie nuclear 
relaxation times and measured meet “2 independently from 
the Raman band shape. They found that the classical 
| J-diffusion theory correctly related these two correla- 
tion times throughout the entire density range studied. 


Gillen et ye have obtained similar results in 
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the liquid region using CCl The spin lattice relaxa- 


4° 
tion time of ys is dominated by spin-rotational relaxa- 


£3 


tion and t. can be calculated from the C relaxation 


J 


: P 3 
times. The transverse relaxation time of the 01 


nucleus is dominated by quadrupole relaxation which 


(2) (2) 
8 v) 


depends on T so that values of T were determined 


35 


from the Cl relaxation times. They found that 


classical-extended J-diffusion correctly predicted the 


(2) 


relationship between To and Ts throughout the entire 


diquid region,of ccl Similar results for less sym- 


4° 
metrical molecules have been obtained in the liquid 


EN aie In addition, several researchers have 


region. 
obtained excellent agreement between experimental and 
predicted infrared and Raman band shapes using the 
elagadeal inctational-diffustonsmodels ici” 
The energy levels of the spin states of the para- 
Magnetic species considered in Chapter IV are dependent 
on the orientation of the molecule relative to the 
magnetic field, thus the reorientational motion of the 
molecule modulates the energies of these spin states 
and the observed ESR spectral linewidths can be inter- 
preted in terms of the tumbling time Cee of the mol- 
ecule. The spin magnetic moment of the molecule is 
also weakly coupled to the rotational magnetic moment 


of the molecule and this spin-rotational interaction 


is modulated rapidly in the liquid systems considered 
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here since the rotational state of the molecule is 
rapidly changed by intermolecular interactions. This 
gives rise to a linewidth contribution (a") governed by 
the angular momentum correlation time TS: 

The ESR linewidth approach does not yield informa- 
tion about the intermolecular potentials directly although 


ad have indicated that the parameter kK 


Kivelson et ‘al. 
does contain information about the anisotropy of the 
intermolecular potential. The reorientational correla- 
tion time te (7), which can be obtained from the linewidths, 
is of direct interest since it indicates the nature of 


(2), 
6 


the reorientational motion. T s determined from ESR 


are often used to obtain information about the active 


aeer spin label (a molecule containing 


sites of enzymes. 
a paramagnetic center) is allowed to bind to the active 
Site. From the T° 78 of the bound and free spin label 
and from other changes in the observed spectra, informa- 
tion about the nature of the binding to the active site 
and about the active site itself may be obtained. For 
example it is possible to distinguish a binding site 
that firmly binds the spin label to the larger molecule 
from a binding site in which the bound spin label is 
allowed considerable freedom for independent motion. 

The aim of the ESR work described in Chapter IV 


: L ; 41-4 
was to further evaluate the Kivelson linewidth theory = 
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which interprets the linewidths of an observed ESR 
spectrum in terms of the modulation of the anisotropic 
electronic Zeeman, nuclear-electronic hyperfine and spin- 
rotational interactions by molecular motion. Several 


tg) 84a AS have noted discrepancies between 


workers 
experimental linewidths and those predicted by the 
Kivelson theory. WE have introduced the procedure of 
analyzing the M-dependent part of the linewidths in 

terms of the Kivelson theory rather than simply consider- 
ing the contribution tothe linewidth quadratic in 

M which was the previous practice. We have justified 
this procedure using computer experiments on test cases 
and feel that the results from studying the M-dependent 
part of the linewidth are much more representative of 
the theory because each parameter is considered in terms 
of its contribution tothe linewidth. We have also 
analyzed single spectra and obtained values of ee 
independent of the Stokes-Einstein hydrodynamic relation- 
ship. °* These results indicate that if the Stokes- 
Einstein relationship is ignored, then the agreement 
between predicted and experimental results is improved. 
Thus it appears that the assumption of the Stokes- 
Einstein relationship is the main cause of the discrepancy 
pecueee predicted and experimental results, however it 


does not explain the discrepancy between the predicted 


and experimental a" which is much greater than that 
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expected from experimental error. Therefore, we conclude 
that although the deficiency in the 8 and y parameters 
could be accounted for in terms of the assumptions of 

the Stokes-Einstein relationship, the simple theory pre- 
dicting the linewidth contribution due to the spin- 
rotational coupling must be reconsidered. 

It would be interesting to study the 636u diketo- 
nate complexes corresponding to those studied in Chapter 
IV. The spin-rotational coupling for copper is much 
larger than for vanadium and thus the a" contribution 
to the linewidth would be proportionately larger for 
the copper complexes. The a"“'s determined from such a 
study, being larger, would have a smaller relative 
experimental error and thus would provide a more stringent 


test of the validity of the expressions used to predict 
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THE WIGNER MATRIX AND ANGULAR MOMENTUM THEORY 


In order to describe the effects of rotational 
motion of molecules on the interactions between the mole- 
cule and an external field, it is often necessary to 
define a tensor in a given coordinate system in terms of 
its components in another coordinate system. For example 
we are often interested in defining an electric field in 
a molecular coordinate system in terms of its components 
in the laboratory coordinate system. Throughout this 
thesis, rotations have been described using the conven- 
tions of Bosal?” In this convention, the coordinate 
systems are rotated and not the functions. The sense 
of a rotation follows the right-hand screw rule i.e. 
roations of a right-handed coordinate system +90° about 
the z-axis rotates the x-axis onto the original positive 
y-axis. The rotation of a coordinate system by an 
amount a is equivalent to the rotation of a function 
by an amount -a thus if the operator R, (a) represents 
the rotation of the coordinate system by an angle a 


- about the z-axis then 


R, (a) £ (4) = f($-a) (A-1) 


where f£(¢) is the function and ¢ is the azimuthal angle 


of polar coordinates. Performing a Taylor series 
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expansion of the function £(¢-a) we obtain 


5 n 
£(o-a) = 1 -a = £ (o>) (A-2) 
Ds 2. ( is) 
= A) 
= €xp (-s2-) Coys 
Ko) 


thus the rotation operator R, (a) is given by exp (—« &). 
In spherical polar coordinates, the reduced angular 


momentum operator J. is given by 


ft of 3 
tT te Rane 
thus the rotation operator is 
R, (a) = exp (-taJ_) “ (A-3) 
Similar considerations yield 
R, (a) = exp (-toJ,) , (A-4) 
and / 
Ry (a) = exp (“tad ' (A-5) 


where J, and Jy are the reduced angular momentum oper- 
ators for the x and y components of the angular 


momentum. 


Any general transformation between two coordin- 


ate systems can be represented by three successive 


rotations. The Eulerian transformation is a most useful 
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way of carrying out such a coordinate transformation. It 
is defined in terms of the following three rotations: 

(1) the coordinate system (x, y, z) is rotated Coroman an 
angle @ about.the z-axis to become (x', y', z'), (2) the 
coordinate system (x', y', z') is then rotated through 

an angle 8 about the y'-axis to become (x", y", 2"), 

(3), the coordinate system (x", y", 2") is then rotated 
through an angle y about the z"-axis to yield the desired 
ecoordvnate-systems (x, vy" 4 2"s)s .Thesang les) a,is6, 

Y = 2 are referred to as the Euler .angles of the trans- 
formation. In terms of the rotation operators defined 


in Eqs. (A-3) - (A-5) the Eulerian transformation is 


“a nan 


R(a,8,Y) = Rou Cy) Ry (B)R, (a) 


exp (-77I.,4) exp (-78I,,) exp(-taJ,). 
(A-6) 


gee that the rotations, described in Eq. 


It can be shown 
(A-6) about axes in different coordinate systems, can 
be expressed as rotations about axes in the initial 


coordinate system, provided that the rotations are 


done in the reverse order. Thus, 


R(a,B,y) = exp (-iaJ,)exp(-t8J,) exp ¢ tyJ,). (A-7) 


The Wigner rotation matrix or {Q-~-matrix is the 
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Matrix representation of the rotation operator in angu- 


lar momentum space and is defined by 


“a (J) - 
R(a,B,y) |T-M) -) Full: ap ayer (a,B,y). (A-8) 
: 7 


where [J.M) is an angular momentum wavefunction. Sub- 
Stituting Eq. (A-7) into Eq. (A-8) and premultiplying 
by (o,M'| , we obtain 


(J) 


D M'M (a,8,Y) = (5 M'|R(a,B,¥) |T-M) 


= exp(-7aM') CoM! | exp(-165)) |5,M) exp(-tyM), 


(A-9) 


where the operator J, has been applied to CoM" | and 
|5,M) . 
From Eqs. (A-6) - (A-9) the relationships [Egqs. 


(III-18) and (III-19)] used in Chapter III can be 


obtained: 
(J) : 
D (a,0,0) = exp(-t) Sy ue (A-10) 
M'M 
and 
(J) ee (J) 
iss (a,B,Y) ‘ = AO eK) 
M'M MM 
(J), 
ees (a,BrY¥) (A-11) 
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f 100 
also it can be shown that 


(J) * M‘'-M (J) 
LB (a,8,y) = (-1) es (a,B,Y) (ASE?) 
M'M -M',-M ee 


The orthogonality relationships [Eq. (III-20) - 
101 
(III-21)] are developed in detail by Rose and will not 
be given here since the development is extensive. 
; 102 : 
Wigner has shown that the ')-matrices are propor- 
tional to the rotational wavefunction for a symmetric 


top, ac¥ This relationship is given by 


JKM* 


OR eee (J) 
¥oxm(OrBry) = ( | Diym (4-B-Y¥) (Ani3) 


where the multiplicative term is a normalizing constant. 
The coupling of two angular moments is described 


by the Clebsch-Gordon series defined by, 


|a,™) =» C(J4I557M, MM |) [7, mM, ) | J, ,M-M, ) 


1 (A-14) 


where 

1 528 >| Saad (A-15) 
and C is the Clebsch-Gordon coefficient. The function 
|5,M) is the angular momentum wavefunction of the 


coupled representation while |, -™,) and |3,.M-M, ) 
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are wavefunctions of the uncoupled representation. 


Using Eqs. (A-8) and (A-14) one can obtain L032 


a (J) ie (J,) 
: eee a (Q) = > C(J,5,5;K,K,)C(JI,I,5;M)M,) 
a o2 5 
(J) . 
yi: (Qn (A-16) 
K,+K,,M,+M, 


This relationship is often useful in simplifying expres- 
sions containing products of -matrices. 

If we define tensors so that they transform under 
rotations in the same way as the angular momentum wave- 
functions do, then all the relationships derived to 
manipulate angular momentum wavefunctions using the 

-matrix can be used to manipulate these tensors. 
Tensors defined in this way are referred to as spheri- 


“a 


cal tensors. The components of a tensor operator T, 


of rank & satisfy the following comutator relation- 


ships! 


ye (£41) -u(utl) Te els (A-18) 
[3., ¥. iy = y2 (241) -u(u-1) La onary 4 (A-19) 


where J, and J are the raising and lowering operators. 
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A first rank spherical tensor operator has components 


wh = qT, (A-20) 


= ae 1 A a 
ay 44 =+ — 2 ORS ae Aa (A-21) 


A 


where ak gq = X,y,Z are the components of the correspond- 
ing cartesian tensor operator. The angular momentum 
operators used in Eqs. (A-17) - (A-19) can also be de- 
fined in terms of the cartesian components. 

In obtaining Eq. (III-6) in Chapter III, we have 
assumed that the mixed products of  -matrices in 
Eq. (III-5) vanished when ensemble averages were per- 
formed. This will now be demonstrated using the rela- 
tionships given above and the orthogonality relation- 

101 


ships given by Rose. The expression of interest is 


the integral over all orientations, 


(2)* (2) 
K =/ dQ (aKM[D = [Q] |T'K'M') CItK'M'|D [2] |IKM) . 
km k'm' 


(A-22) 


Henceforth the argument 2 of the -matrices will not 
be shown explicitly since it is 2 in all cases. Sub- 
stituting Eq. (A-13) into Eq. (A-22) and using Eq. (A-12) 


we obtain 
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' 1_qtwz! (7 )* (Uy* 
Z - faa (254+1) (25'4+1) Pe shied aia » 
(3n2) 2 -K-M km 
pee) (2) (J) 
D D D D : (A-23) 
K'M! -K'-M! k'm'  -K-M 


After applying the coupling expression [Eq. (A-16)] we 


obtain 


K = (23+1) (23'+1) (-1) M'-K! C(IRI;-Kk) “C(T2RI;—Mm) * 


(807) 2 Jae okt 


KC oO 8 hk pk’) C(I os" eM MY) C (lds: KK) Cas em aM) 


(J) $s (J") 
[es D D JB) (A-24) 
KM 00 K"M" 
where K = -Ktk K" = -K+k' 
(A-25) 
M = -M+m M" = -Mtm' . 


The orthogonality relation for a triple product of D- 
101 


matrix components given by Rose reduces the integral 
in Eq. (A-18) to 
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Thus the integral vanishes unless K' = K and M' = M. 
Combining these conditions with (Eq. (A-25) we find the 
integral vanishes unless k = k' and m = m' which is the 


required condition assumed in Chapter III. 
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This Appendix gives the FORTRAN computer pro- 


grams used in the course of this work. 


= 


A. Semiclassical Band Shape Program (Gordon method). 


B. Program to calculate the free rotor memory 


function from the free rotor correlation function. 


C. Program to calculate the correlation function 


from the free rotor memory function. 


D. Linear and non-linear least squares programs 


used in Chapter IV. 


1. Program to fit linewidths to a cubic 


function in M. 


2. Program to fit M-dependent part of the 


linewidth. 
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SECTION A 223. 


THIS PROGRAM CALCULATES SPECTRAL DENSITIES AND REORIENTATIONAL 
CORRELATION FUNCTIONS USING EITHER M OR J DIFFUSION 


CARD #1 : TEMPERATURE; TAUJ; HEIGHT OF HEIGHT NORMALIZED 
SPECTRA; MAXIMUM J QUANTUM NUMBER (DEFAULT=9 
: MAX=12); 
IPLOT ( CONTROL PARAMETER); FN M OR J ? 
TEMP <0.0 STOP. 
FM=1 100% M-DIP 
FM=0 100% J-DIPF i 
IPLOT 0=-YES;""~te- NO 
1. HEIGHT ? 
2. G(T) ? 
3. AREA PLOTS ? 


CARD #2: TITLE TO APPEAR ON PLOTS AND OUTPUT 
COMPLEX*8 SM(3,3) ,MLAM(3) ,DUMP3 (3,3) 
REAL OMO (37) ,D (37) ,P(37) ,RI (345) ,OMOW (345) ,WH(191) ,BUFFER(1024), © 
# TITLE (22) ,XX (191) 
COMPLEX*8 S (37,37) , SINV(37,37) ,LAM (37) ,C1,SUMA,SUMB, A(37) .B(37), 
# RC (345) ,DUMMY(37,37) 
REAL NORM, DEG (13) ,XD(2) 
LOGICAL INTH (37) 
INTEGER INT (37) ,IPLOT (3) 
DATA OMO/3028. 740, 3038.500,3048.180,3057.720,3067. 240, 3076.650, 
# 3086.410,3095.190,3104.450, 3113. 430, 3018.810, 3018. 640,30 18.380, 
# 3018.620, 3017.470,3016.940, 3016.510, 3016.070, 3015.610,3009.000, 
# 2999.010, 2988.920,2978.880,2968-630,2958.340,2948.060,2937.830, 
# 2927.200,3122.75, 3014. 13, 2916.72, 3131.87, 3013.29, 2906.20, 
# 3140.311, 3012.44, 2895.62/ 
DATA TITLE/*INTE',*NSIT',*Y  *,*FREQ*',*UENC',"Y °,16*! 7 
DATA D/1.0000000, 0. 7453560, 0.6831300,0.6546536,0.6 382847, 
0. 6276459,0.6201736,0.6186362,0.6103679,0.6069769,0.5773502, 
0.5773502,0. 5773502, 0. 5773502, 0. 5773502, 0.5773502,.5773502, 
0.5773502, 0. 5773502, 0. 3333333,0.4472135,0.4879500,0.5091750, 
0.5222329,0. 5310850, 0. 5374838, 0. 5423261, 0.5461186,0. 6042180, 
0.5773503,0. 5491696, 0.6019293,0.5773503, 0.551677 3,0. 6,0. 5773503, 
0.5537749/ 
DATAs DEG/5~ 0,90, 25.0,.77-0,117.0,121.0,273.0,285.0,357.0,513.0, 
# 609.0,621.0,925.0/ 
DATA NORM/*NORM!/ 


th SS te Me Se Gh 


INITIATION OF PLOTS 


REAL YD(2)/0.0,2-0/ 


Ic=0 
YORIG=- 10.0 
XORIG=17.0 


XSIZE=14.9606 
YSIZE=7. 874 

CALL PLOTS (BUFFER, 4096) 
CALL PLOT(2.0,2-0,-3) 

1 READ (5,901,END=29) TEMP,TAUJ,HEIGHT,J, (IPLOT(I) ,I=1,3) ,FM 
IP (TEMP.LT.0.0) GO TO 29 
READ (5,902) (TITLE(I) -I=7,22) 
N=28 
IF (J-NE.0) N=JtJ+3+1 
JJI=3 


CALCULATION OF POPULATIONS 
DEG=((2J+ 1) **2) (NUC.SPIN.FACTOR) 
P=DEG* EXP (-E/KT) 
~E/KT=-7.51702*3* (J+1) /TEBP 


EPAC=-7.51702/TEMP 
P (1) =DEG (1) 
SUM=DEG (1) _ 
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DO 2 J=2,10 
TERM=DEG (J) *EXP (EPAC*FLOAT ( (J-1) *J)) 
K=J+9 
M=J+18 
P (J) =TERN 
P(K) =TERS 
P (M) =TERS 
2 SUM=SUM+TERN 
IP (N.EQ.28) GO TO 4 
NJ=N/3 
KKK=NJ+ 1 
K2=28 
DO 3 J=115KKK. 
K=K2¢1 
K1=K+1 
K2=K+2 
TERM=DEG (J) *EXP (EFAC*PLOAT ( (J-1) *J) ) 
P(K) =TERM 
P(K1)=TERM 
P (K2) =TERM 
3 SUM=SUM+TERM 
4 CONTINUE 
DO 5 J=1,N 
5 P(J)=P(J)/SUM 


SETTING UP S-MATRIX 


DEM= (1.88496E11) *TAUJ 

P (1) =3.0*P (1) 

DO 6 I=1,H8 

DO 6 J=1,¥8 

DEL=0.0 

IF (I.EQ.J) DEL=1.0 

REALS= ( (DEL-0. 3333333*P (J) *D(J) /D(I)) /DEM) * (FM-1.0) 

IF (P(J3).EQ.P(I)) REALS=REALS-FM* (DEL- (D(J)/(3-0*D (I))))/DEM 

COMPS=DEL*OMO (1) “7 
6 S(J, I) =CMPLX (REALS, COMPS) 

P (1) =P(1)/3.0 


DIAGONALIZATION 
IF (FM.NE.1.0) GO TO 15 
4-DIFFUSION 


DO 10 IDIG=1,9 
KK=IDIG+1 
DO 7 I=1,3 
K=KK+9* (I-1) 
SM (1,1) =S(K,KR) 
Sm (1,2) =S(K,KK+9) 
7 SM (1,3) =S(K,KK+18) 
CALL ALLMAT (SM,MLAM,3,3,NCAL, DUMP3, DUNMY,A,B, INTH, INT) 
WRITE (6,903) NCAL 


PLACING INFORMATION BACK INTO S-MATRIX 


DO 8 I=1,3 

K=KK ¢9* (I-1) 

S (K, KK) =S(I, 1) 

S (K, KK+9)=SM(I,2) 

S (K, KK+ 18) =SM (I, 3) 
8 LAM(K)=MLAM (I) 

CALL CHINVS (SH, 3, 3) 

Do 9 I=1,3 

K=KK+9* (I-1) 

SINV (K, KK) =SM (I, 1) 

SINV (K, KK¢9) =SB(I,2) 
9 SINV (K,KK¢18) =SM (I, 3) 
10 CONTINUE 
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LAM (1)=S(1,1) 
S(1, 1)=(1.0,0. 0) 
SINV (1,1)=(1.0,0.0) 
IP (N.EQ.28) GO TO 17 
KK=26 
K=28 
DO 14 IDIG=10,JJJ 
KK=KK+3 
DO 11 I=1,3 
K=K+1 
SM (I, 1) =S(K,KK) : 
SM (I,2)=S(K,KK+1) 
11 SM (I, 3) =S,(K,KK+2) 
CALL ALLMAT(SN,MLAN,3,3,NCAL, DUNP3, DUMMY,A,B, INTH, INT) 
WRITE (6,903) NCAL 
K=K-3 
DO 12 I=1,3 
K=K+1 
S (K, KK) =SM(I,1) 
S (K, KK+1)=SM (1,2) 
S (K, KK+ 2) =SM (I, 3) 
12 LAM (K)=MLAM (I) 
CALL CMINVS(SN,3, 3) 
K=K-3 
DO 13 I=1,3 
K=K+1 
SINV (K, KK) =SM (I, 1) 
SINV (K,KK+1) =SH(I, 2) 
13 SINV(K,KKt2)=SM(I,3) 
14 CONTINUE 
GO TO 17 
15 CONTINUE 


J-DIFFUSION 
CALL ALLMAT (S,LAM,N,37,NCAL,SINV, DUMMY,A,B,INTH, INT) 
WRITE (6,903) NCAL 


THIS SECTION INVERTS THE S-MATRIX 


DO 16 I=1,8 

DO 16 J=1,N 
16 SINV (I,J) =S (I,J) 

CALL CMINVS (SINV,N,37) 
17 CONTINUE 


CALCULATION OF LINE SHAPES 
C1= (0.0, 1.0) . 
DOT PRODUCTS 


DO 19 J=1,H 
SUMA=(0.0,0.0) 
SUMB=(0.0,0.0) 
LAM (J) =LAM (J) *C1 
DO 18 I=1,N 
SUMA=SUMA+D (I) *S(I,J) 

18 SUMB=SUMB+SINV (J, I) *P(I) *D (I) 
A (J) =SUMA 

19 B(J) =SUSB 


ACTUAL CALCULATION OF BANDSHAPE 


ON=2844.0 
DO 21 I=1,385 
OM=ON+1.0 
SUMA=(0.0,0.0) 
DO 20 J=1,8 
20 SUMA=SUMNA+A (J) *B(J) / (OM*#LAB (J) ) 


22555 


(Pwd HT OAT ne pent 


*« 
zz. .e 
oO hy as 
Ne , 
A = ; 


BSA a 


Noha 
; oe 


io 
he 


p »@ 7 : ite i 
s ~ ' )-O2A Thy ih 7) 4 
Ss ® i f iu > ‘ 


- - tar - ; Pr ' ¥ bo a 
i ree 


es os 


aaAA 


SUMA=SUMA/(5.92176211) 226. 


RI (I) =AIMAG (SUMA) 
RC (I) =SUMA 
21 OMOW (I) =OM 
WRITE (6,904) (TITLE (I) ,I=7,22) 
WRITE (6,905) FM 
WRITE (6,906) TAUJ,TEMP 
WRITE (6,907) (OMOW(I) ,RC(I) ,I=1,345) 
RMAX=RI (69) 
DO 22 I=100,250 
IF (RMAX.GT.RI(I)) GO TO 22 . 
RMAX=RI (1) . 
WMAX=ONOW (I) 
22 CONTINUE 
RAT=HEIGHT/RMAX 
DO 23 J=1,345 
23 RI(J)=RI (J) *RAT 
WRITE (6,908) HEIGHT 
WRITE (6,909) (OMOW(I),RI(I) ,I=1,345) 


PLOTS 


DO 24 J=1,191 
K=J+83 
WW (J) =OMOW (K) 
24 XX (J) =RI (K) 
XD (1) ="W (1) 
XD (2) =WW (191) 
IF (IPLOT(1).£0Q.1) GO TO 25 
CALL PLOTME (XSIZE,YSIZE,2,XD,YD,TITLE,. FALSE. ,0. 0, 4) 
CALL PLOTME (XSIZE,YSIZE, 191, WW, XX, TITLE,. TRUE. ,0.0,4) 
CALL RDORIG (IC, YORIG, XORIG) 
25 CONTINUE 
RATNOR=RAT 
CALL DPCOR(OMOW,RI,WMAX, TITLE, TEMP, RAT ,IPLOT) 
ARAT=31.41593/RAT 
REALA=RAT/RATNOR 
TP=ARAT/ (HEIGHT*1.884956E12) 
WRITE (6,910) REALA,RAT, TP 
DO 26 J=1,191 
26 XX (J)=XX (J) *ARAT. 
WRITE (6,911) 
WRITE (6,909) (WwW(I),XX(I) ,I=1,191) 
IF (IPLOT(2).EQ.1) GO TO 27 
CALL RDORIG(IC,YORIG, XORIG) 
27 CONTINUE 
TITLE (22) =NORM 
IP (IPLOT(3)-EQ.1) GO TO 28 
CALL PLOTME (XSIZE,YSIZE,2,XD,YD, TITLE,. FALSE. ,0.0, 4) 
CALL PLOTME(XSIZE, YSIZE, 191, WW,XX,TITLE,. TRUE. ,0.0,4) 
CALL RDORIG (IC,YORIG,XORIG) 
28 CONTINUE 


GO TO 1 
29 CALL PLOT(0.0,0.0,999) 
STOP 


901 FORMAT (F10.0,E10.3,F10.0,13,511,F10.0) 


902 FORMAT (* *,16A4) 


903 FORMAT (* NUMBER OF EIGENVALUES FOUND=',13) 


904 FORMAT (* 1, 16A4) 


905 FORMAT (* FRACTION OF M-DIFFUSION=*,P5.2) 
906 FORMAT (* CALCULATED INTENSITIES ','J-DIFFUSION TAUJ='",£13.3, 
#°SECONDS, TEMP=',F6.1//*FREQUENCY INTENSITIES (REAL, IMAGINARY) °) 


907 FORMAT (2(F10. 1,2E20.6)) 


908 PORMAT ('1SPECTRA NORMALIZED TO MAXIMUM HEIGHT OF',F10. 5) 


909 FORMAT (5(F10.1,F10.3)) 


910 FORMAT (‘AREA UNDER ORIGINAL SPECTRAL DENSITY',E20. 6/" AREA UNDER B 
#ORMALIZED(HEIGHT) SPECTRAL DENSITY',E20. 6/*REORTENTATION CORRELATI 


#ON TIME',E20.6,'/SEC") 
911 FORMAT (‘1 SPECTRA NORMALIZED TO AREA OF 10 oe 


END 
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SUBROUTINE DPCOR(W,INT,WMAX, TITLE, TEMP,RAT,IPLOT) 


LCULATION OF DIPOLE CORRELATION PUNCTIONS 


WB= (K*T/I) **0.5 = ((0.2565E24) * TEMP) **0.5 


K 
I 


901 
# 
# 

902 


= 1.38E16 GRAM (CM/SEC) **2 
= 5.38E-40 GRAM(CM) **2 


REAL INT (345) ,W¥(345),G(101),TT (101) 
REAL SCALE(22) , TITLE (22) 

REAL SG (101) 

DATA SCALE/! ',"G(T)t,? ', *,*TIME', 17%! 
REAL XD(2)/0.0,7.6/ 

REAL YD(2)/0.0,1.0/ 

INTEGER IPLOT(5) 
WB=SQRT ( (0. 2565E24) *TEMP) 

fol ts? 22 

SCALE (I) =TITLE (I) 

DO, 3 IT=1,101 

RT=0. 1*PLOAT (IT-1) 

TT(IT)=RT 

RT=RT/WB 

SUM=0.0 

SSUM=0.0 

FPRT= 1.884956E11*RT 

DO 2 IW=1,345 
SSUM=SSUM+INT (IW) *SIN(FRT* (W(IW) -WMAX) ) 
SUM=SUM+INT (IW) *COS (FRT* (W (IW) -WMAX)) 


PI*SPEED LIGHT = 1.884955592E11 


IF (IT.EQ.1) RAT=SUM 
SG (IT) =SSUM/RAT 

G (IT) =SUM/RAT 

WRITE (6,901) WMAX, WB 

WRITE (6,902) (TT (I) ,G(I),SG(I) ,I=1,101) 
IF (IPLOT(2).EQ.1) RETURN 
XSIZE=38.0*0. 3937 

YSIZE=20.0*0. 3937 


CALL PLOTME (XSIZE,YSIZE,2,XD,YD,SCALE,. FALSE. ,0.0,4) 


CALL PLOTME(XSIZE,YSIZE,77,TT,G,TITLE,. TRUE. ,0.0,4) 


CALL PLOTME (XSIZE,YSIZE,77,TT,SG,TITLE,.TRUE.,0.0, 4) 


RETORN 
FORMAT (*1 DIPOLAR CORRELATION FUNCTION CALCULATED 


ABOUT*,F10.1," CM-1'/* TIME IN REDUCED UNITS OF 1/OMEGABAR, 


OMEGABAR=",E16.6) 
FORMAT (4(F5.1,2F7.4)) 
END 
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228. 
SECTION B 


THIS PROGRAM READS THE TIME(REDUCED UNITS) AND 

THE FREE ROTOR CORRELATION FUNCTION AND ITS FIRST 
DERIVATIVE ON 3. IT WRITES TIME AND THE ASSOCIATED 
MEMORY FPUNCTION ON 4. . 
NUCLEAR SPIN SYMMETRY ORDER A,E,P 


TIME (502) =OMEGABAR 

TIME (503) =TEMP 
G(502)=AVERAGE FREQUENCY 
GP (502)=K (0) 


ON 5 
CARD #1 3: NT; NUMBER OF INCREMENTS OF DEL T TO 
BE CONSIDERED 


COMPLEX*8 KA (502) ,KE(502) , KF (502) ,G (502) ,GP(502) 
REAL TIME(503) 

READ (5,901) NT 

READ (3) TIME 

WRITE (4) TIME 

DELT=TIMNE(2) 


READ (3) G 

READ (3) GP 

CALL GFKF(G,GP,NT,KA,DELT) 
KA (502) =G (502) 

WRITE (4) KA 


E 


READ (3) G 

READ (3) GP 

CALL GPKF(G,GP,NT,KE,DELT) 
KE (502) =G(502) 

WRITE (4) KE 


READ. (3) G 

READ (3) GP 

CALL GFKF(G,GP,NT,KF,DELT) 

KF (502) =G (502) 

WRITE (4) KP 

RAW=REAL (KA (502)) 

REW=REAL (KE (502) ) 

RFW=REAL (KF (502) ) 

WRITE (6,902) TIME(503) ,RAW,REW,RPW 
WRITE (6,903) (TIME (I) -KA(I) ,KE(I) ,KF(1) ,I=1,NT) 
STOP 


901 FORMAT (I4) 
902 FORMAT ("1 FREE ROTOR MEMORY FUNCTIONS FOR NU-4 IR BAND OF*, 


€* METHANE : TEMP=',P10.1,* DEG K'/'- CALCULATED ABOUT ', "AVERAGE 


# FREQUENCY OF THE STICK SPECTRA'/*' A ‘',F10.3/' E *,F10.3/ 
#° F *, F10.3) 


903 FORMAT ('- TIME(RED)*,10X,"A-SYMMETRY',19X,"E-SYMMETRY',19X, 


#*P-SYMMETRY'/(F8.2, 3X, 3(2E13.4,3X))) 
END 
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THIS SUBROUTINE CALCULATES THE FREE ROTOR 
MEMORY FONCTION FROM THE FREE ROTOR FUNCTIONS 


SUBROUTINE GPKF(GT,GP,ITINE,K,DELT) 
COMPLEX*8 GT (502) ,GP (502) ,K(502),C1 
K (1) =-GP (502) 

DELT2=DELT/2.0 

K (2) =-(GP(2)-GP(1)) /DELT2-GT(2) *K (1) 
DO 2 N=3,ITIME 

NM1=N-1 | 

c1=0.0 : 

NP1=Ne1 

DO 1 M=2,NM1 

C1=C 14#GT (NP 1-4) *K (M) 

CONTINUE 

K (N) =-C1-C1- (GP(N) -GP(1) ) /DELT2-GT (N) *K (1) 
CONTINUE 

RETURN 

END 
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SECTION C 230. 


CALCULATION OF REORIENTATIONAL CORRELATION FUNCTION FROM FREE 
ROTOR MEMORY FUNCTION AND THE MEMOBY PUNCTION EQUATION 
K (T) =K (T) FR*EXP (-T/TAUJ) 
PREE ROTOR MEMORY FUNCTION READ IN ON 4 
K (502) =AVERAGE FREQUENCY OF STICK SPECTRA 
TINE (502) =OMEGABAR; TIME(503) =TEMP 
COMPLEX*8 CPHASE(501) ,CEXP,CMPLX,CA,CE,CF,CPHA (501) ,CPHF (501), 
# KAT(501),KET (501) ,KFT (501) ,GA(501) ,GE(501) ,GP (501) 
COMPLEX*8 K (502),C1,GT (501) ,KP (501) ,C3,KA (502) yKE(502) , KP (502) 
REAL TIME(503) ,BUPPER(1024) ,KR (501) ,KI (501) , TITLE(22) , TITLEK (22) 


INITIATION OF PLOT ROUTINES 


CALL PLOTS (BUFFER,4096) 

DATA TITLEK/* K(',"T)/K*,*(0) *,* TI","ME(R*,"ED) *, 16¥*° W4 
DATA TITLE/* *,"G(T)*,! '," TI',*ME(R',*ED) *,16*° 7 
REAL YD(2)/0.0,1.0/ 

REAL YDK (2) /-0.5,1.0/ 

REAL XD(2)/0.0,7.6/ 

REAL XC(3)/0.0,4.0,7.6/ 

REAL YC(3)/0-.0,0.0,0.0/ 


Ic=0 
YORIG=- 10.0 
XORIG=17.0 


XSIZE=14.9606 
YSIZE=7. 874 
CALL PLOT (2.0, 2.0,-3) 
READ (4) TIME 
READ (4) KA 
READ (4) KE 
READ (8) KP 
WB=TIME (502) 
WA=R EAL (KA (502)) 
WE=REAL (KE (502)) 
WF=REAL (KF (502)) 
WPB= (5.0*WA+2. OFWE+3.0*WF) /10.0 
WA=WPB-WA 
WE=WPB-WE 
WF=WPB-WP 
TEMP=TIME (503) 
DELT=TIME (2) 
B1=5.194 
CPHASE(1)=1.0 
CPHA (1) =1.0 
CPHF (1) =1.0 
CA=C EXP (CMPLX (0.0, (B1#B1+WA) *DELT*1.885E11/WB) ) 
CE=CEXP (CMPLX (0.0, (B1+B1+WE) *DELT*1.885E11/WB) ) 
CF=C EXP (CMPLX (0.0, (B1#B1+WF} *DELT* 1.885E11/WB) ) 
DO 1 I=2,501 
CPHA (I) =CPHA (I-1) *CA 
CPHASE (I) =CPHASE (I- 1) *CE 
CPHF (I) =CPHF (I-1) *CP 
1 CONTINUE 
DO 2 I=1,501 
CPHA (I) =CPHA (I) *5.0 
CPHASE(I)=CPHASE(I) *2.0 
CPHF (I) =CPHF (I) *3.0 


2 CONTINUE 
TAUJ (SEC) 

TMAX (RED) MAXIMUM TIME TO BE CONSIDERED 
IPLTK PLOT K(T) ? 

1¢ G(T) ? 

IPLTG PLOT G(T) ? 

PAR=0 YES 

PAR=1 ne) 
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3 READ (5,901) TAUJ,ITNAX,IPLTK,IG,IPLTG 
IF (TAUJ.LT.0.0) GO TO 8 
TMAX=TIME(ITMAX) 

READ (5,902) (TITLE(I) ,I=7,22) 
DO 4 I=7,22 
TITLEK (I) =TITLE(I) 

4 CONTINUE 

TJRED=TAUJ*WB 


MULTIPLYING K(T)-FR BY EXP(-T/TAUJ) 


RATK=REAL (K (1) ) 
R1=EXP (-DELT/TJRED) 
R2=1.0 ¥ 
KAT (1)=KA(1) 
KET(1)=KE(1) 
KPT (1) =KF(1) 
DO 5 I=2,501 
R2=R2*R1 
KAT (I)=KA(I) *R2 
KET (I)=KE(I) *R2 
KFT (I) =KF (I) *R2 
5 CONTINUE 
TF’ (IG.£OQ.1) GO To 3 
WRITE (6,903) (TITLE(I) ,I=7,22) 
WRITE (6,904) TEMP, TAUJ, TIRED 
WRITE (6,905) 


WRITE (6,906) 
CALL KG (KAT,GA,DELT,1ITMAX, TINE) 


E 4 


WRITE (6,907) 
CALL KG(KET,GE,DELT, ITMAX, TIME) 


WRITE (6,908) 
CALL KG(KFT,GP,DELT, ITMAX, TIME) 
GT (1) =5. O*GA(1) +2.0*GE (1) +3.0*GP (1) 


NORMALIZING G(T) AND PHASING SO FIRST SOMENT 
OF THE SPECTRAL DENSITY = 2 B1 


RRR=REAL (GT (1) ) 
GT(1)=GT (1) /RRR 
C3=0.5 
DO 6 I=2,ITMAX 
C1= (CPHA (I) *GA (I) *CPHASE (I) *GE (I) #CPHF (I) *GP(I)) /RRE 
GT (I) =C1 
C3=C34+C1 
6 CONTINUE 
TAUTHA=REAL (DELT* (C3-GT (ITMAX) /2.0))/(1.-REAL (GT (ITMAX) )) 
TWB=TAUTHA/WB , 
WRITE (6,909) 
WRITE (6,910) (TIME (I),GT(I) ,I1=1,ITMAX) 
WRITE (6,911) TAUTHA,TWB 
IP (IPLTG.£Q.1) GO TO 3 
J=IFIX (7.6/DELT) 
IF (TMAX.LT.7.6) J=ITMAX 
CALL PLOTME (XSIZE,YSIZE,2,XD,YD,TITLE,.FPALSE.,0.0,4) 
IJ=3+1 
DO 7 I=1,Td 
KR (I) =REAL (GT (I)) 
KI (1) =AIMAG (GT (I) ) 
7 COXTINUE 
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CALL PLOTME (XSIZE,YSIZE,J,TIME,KR,TITLE,.TRUE. ,0.0,4) 
CALL PLOTME(XSIZE,YSIZE,J,TIME,KI,TITLE,.TRUE. ,0.05,4) 
CALL RDORIG (IC,YORIG, XORIG) 


GO TO 3 
8 CALL PLOT(0.0,0.0,999) 
STOP 


901 FORMAT (£10.0,I4, 311) 

902 FORMAT (16A4) 

903 FORMAT (*1',16A4) 

904 PORMAT ('-', 
#*REORIENTATIONAL CORRELATION PUNCTION CALCULATED FROM',* FREE ROTO 
#R MEMORY FUNCTION'/* TENP=",F8.2," DEG K; TAUJ=",E10.3," SEC, =%, 
#£10.4," (RED) ') ¥ 

905 FORMAT (* CORRELATION FUNCTION CALCULATED FROM MEMORY',* FUNCTION 
#BY TRAPEZOIDAL RULE SQRT(G(T)*G(T)) & G(T) *G(T) *) 


906 FORMAT ('- A-SYMMETRY'") 
907 FORMAT (‘1 E-SYMMETRY*) 
908 FORMAT ("1 F-SYMNNETRY*) 


909 FORMAT ("1 GA,GE,GF: WEIGHTED SUM(5:2:3) ;PHASE SUCH THAT THE®, 

.. #* AVERAGE FREQUENCY IS WEIGHTED SUM OF WPA,WPE,WPF*) 

910 FORMAT ('- REORIENTATIONAL CORRELATION FUNCTIONS ADJUSTED so*, 

_  #* PIRST MOMENT=2B1=10.388 CM-1'/(2(F10.2,5X,2E15.7, 10X))) 

911 FORMAT (*-REORIENTATIONAL CORRELATION TIME= ',E11.3,'(RED) = °, 
#E11.3,* SEC") 
END 


THIS SUBROUTINE CALCULATES G(T) FROM K(T) 


SUBROUTINE KG(KP,GT,DELT,ITMAX, TIME) 
COMPLEX*8 KP(501) ,GT(501) ,C3,KSUM (501) ,C1 
REAL TIME(503) 
WB=T IME (502) 
RC1=DELT*DELT/2. 
RC2=REAL (1.0+RC1*KP (1) /2.) 
C3=(-1.0+#RC1* (KP (2) +KP(1) /2-) ) /RC2 
RC4=RC1/RC2 
KSUM (1) =C3 
ITMN1=ITMAX-1 
DO 1 I=2,ITH1 
KSUM (I) =RC4*(KP(I+1) +KP(I)) 
1 CONTINUE 
GT(1)=1.0 
GT (2)=(1.0-RC1*KP (2) /2.0) /RC2 
WRITE (6,901) TIME(1),GT(1) 
WRITE (6,901) TIME(2) ,GT (2) 
C3=GT (1) /2.0+GT (2) 
DO 3 N=3,ITMAX 
IF (TIME (N).LE.7.6) IJ=8 
C1=-KSUM (N-1) /2. 
NM2=N-2 
DO 2 J=1,NM2 
C1=C1-GT (N-J) *KSUM (J) 
2 CONTINUE 
GT (N)=C1 
C3=C3+C1 
AAAA=CABS (GT (N)) 
AAAAA=AAAA*AAAA 
WRITE (6,901) TIME(N) ,C1,AAAA, AAAAA 
3 CONTINUE 
TAUTHA=REAL (DELT* (C3-GT (ITMAX) /2.0))/(1.-REAL (GT (ITMAX))) 
TWB= TAUTHA/WB . 
WRITE (6,902) TAUTHA,TWB 
RETURN 
901 FORMAT (F10.4, 4E16.7) 
902 FORMAT ('-REORIENTATIONAL CORRELATION TISNE= ',E£11.3,' (RED) = ', 
#£11.3," SEC") is 
END 
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~SECTION D =] 233. 


INTEGER ISKIP(8),MWT (8) 


THIS PROGRAM FITS OBSERVED EXPERIMENTAL ESR 
LINEWIDTHS TO A CUBIC FUNCTION IB &. 


USES DR. BIRSS*'S LEAST SQUARES ROUTINE ala les crf 
TO WEIGHT POINTS. 


INPUT ON 5 
CARD #1 : IMOD,IWT 
- IMOD=1 READS ISKIP ON 3 
IWT=1 WEIGHTS POINTS AS 1/(WIDTH**3) 
=0 NO WEIGHTING 

CARD #2 : TITLE UP TO 40 CHARACTERS 

CARD #3 : PPW,TC,ETA 
PPW, PEAK TO PEAK WIDTHS IN GAUSS 
PPW(1)=0.0 STOP 
PPW(1)=100.0 READS NEW TITLE CARD 
TC=TEMPERATURE IN DEGREE C. 
ETA=VISCOSITY IN CENTIPOISE 


INPOT ON 3 

ISKIP, INTEGER VECTOR 

IF ISKIP(I)=1 I-TH LINE NOT CONSIDERED IN FIT 
OUTPOT ON 9 SUMMARY OF PARANETERS 
OUTPUT ON 6 LSQ OUTPUT 


REAL D(11),A(10),PPH(8) , PPW (8) 

REAL TITLE (10) 

REAL INC(11)/* t," NIC',Q*# ‘/ 

REAL M(8)/-3-5 ¢—2e5S ye Sy V0 57 0054 1054 20 5p 30 5/ 

REAL MM (8) /12% 25, 6.25, 2.25,0 025 40025,2025 96-25 412025/ 
REAL MMM (8) /-42. i515. 625,-3. 475-125). 125, 32375,15-625,42- 875/ 
DATA ISKIP/0,0,0,0,0,0,0,0/ 

IW (A) =IFIX (0. 25E6/ (A**3) ) 

WRITE (9,901) 

D(1) =1.0 

READ (5,902) IMOD,IWT 

GO TO 2 

READ (5,903) (PPW(I),I=1,8) ,TC,ETA 

IF (PPW(1)-EQ.0.0) GO TO 11 

IF (PPW(1).NE.100.0) GO TO 3 


2 READ (5,904) ISOL,TITLE 


WRITE (9,905) TITLE 
Go EO al las 


3 CONTINUE 


WRITE (6,906) TITLE 

IF (IMOD.EQ.0) GO TO 5 

READ (3,907) (ISKIP(I) ,I=1,8) 
GOsnO ob | 6: 

ISKIP(8) =0 

IF (PPw(8).EQ.0.0) ISKIP(8)=1 
CONTINUE 

DO 7 I=1,8 

IF (ISKIP(I).EQ.1) GO TO 7 

D (2) =4(I) 

D (3) =4M (I) 

D (4) =SMM (I) 

D (5) =PPW (I) 

MWT (I) =I8 (PPW(I)) 

IP (IWT.EQ.0) MWT(I)=0 

CALL LSQSM (4,D,MWT(I)) 
CONTINUE 
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ANS 


CALL LSQSY (A) 


A (1) =ALPRA 
A(2) =BETA 

A (3) =GAMMA 
A(4) =DELTA 


TK=TC#273.15 

GOB=A (3) /A(2) 

EOT=ETA/TK 

WRITE (6,908) GOB,EOT,ETA,TC,TK 


THEORETICAL WIDTHS AND STD(PROM EXPERIMENT) 


DO 8 I=1,8 
8 PPH(I)=A (1) +A (2)*M(I) A (3) *MM (I) +A (4) *MMM (I) 
WRITE (6,909) 
SUM=0.0 .... 
COUNT=0.0 
DO 10 J=1,8 
CAL= PPW (J) 
THE= PPH (J) 
DIP=CAL-THE 
IP (PPW(J).EQ.0.0) GO TO 9 
SUM=SUN+ DIF*DIF 
COUNT=COUNY+1.0 
GO TO 10 
9 DIF=1.0E20 
CAL= 1.0520 
10 WRITE (6,910) CAL,THE,DIP,INC(ISKIP (J) +1) 
SUM= (SUM/COUNT) **0.5 
WRITE (6,911) SUM 
WRITE (9,912) TK,EOT, (A(JJ) ,JJ=1,4) ,GOB, ETA 
IF (IMOD.EQ.0) GO TO 1 
GO TO 4 
11 WRITE (9,913) 
STOP 


234. 


901 FORMAT (* TEMP DEG K ETA*/TEMP',6X,'ALPHA',9X,*BETA*,10X, 


_  #"GAMMA', 9X,"*DELTA',4X,'GAMMA/BETA *,4X,"ETA*) 
902 FORMAT (213) 

903 FORMAT (10F7.0) 

904 FORMAT (10A4) 

905 FORMAT (* *,10A4) ; 


906 FORMAT (*'1',10A4) 


907 FORMAT (811) 

908 FORMAT (* GAMMA/BETA',E16.8/" ETA/TENP',E18.8/" ETA(CP) 
#* TEMP(C) ',F18.2/" TEMP(K) ',F18.2) 

909 FORMAT (*- CALCULATED THEORETICAL DIFFERENCE‘) 

910 FORMAT (* *,F10.3,2F13.3,A4) 


°911 FORMAT (*-RMS DIFPERENCE=',F10.5) 


912 FORMAT (*.',F8.2,E£12.3,F13.3,2P14.4,P15.4,2P11.4) 


913 PORMAT (" ALL VALUES OF ETA IN 10 TO -3 CENTIPOISE') 


END 
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SECTION D-2 FA sa FA 


THIS PROGRAM DOES A NONLINEAR LEAST SQUARES FIT TO THE 
MN-DEPENDENT PART OF THE LINEWIDTH TO KIVELSON'S EXPRESSIONS 
WITH STOKES-EINSTEIN. - 
VARIABLE PARAMETER KAPA*R**3 


USES KARL HALLINS'S DRIVER ROUTINE FOR BIRSS'S LINEAR PROGRAM. 


ABSOLUTE RESIDUALS 


INPUT 

CARD #1 :-AX,AY,AZ,GX,GY,GZ (A'S AND WO IN (1/SEC) *E-9) 
CARD #2 : WO,BO (BO IN GAUSS) 

CARD £3 : TITLE (ENDFILE STOP) 


CARD #4 : PPW,ETA,TEMP, ALPHA (EXPERIMENTAL) 
PPW(1)<O STOPS READING DADA 
( THERE WILL BE ONE OF THESE CARDS FOR EACH SPECTRUM) 
IMPLICIT REAL*8(D) 
REAL*4 TITLE(10) , TEMPV (301) 
REAL*8 D(11),DM(301) ,DA(10) ,DX (301) ,D¥ (301) , DPY 
REAL*8 F,U,TR, WTR2,PP& (8) , ALPHA, DEOT 
DIMENSION MUWT (301) 
COMMON /NAME/ NT,MB,NPT,NPAR 
COMMON /PARDPY/ P,U,TR,WTR2 
DATA DS/0.0D0/,DO0T/0.0D0/,NF/0/ 
REAL*8 MVEC (8) /-3.5D0,-2.5D0,-1.5D0,-.5D0,.5D0, 1.5D0,2.5D0,3.5D0/ 
EXTERNAL DPY 
NT=0 
NPT=0 
CALL COMPUT 
1 READ (5,901,END=8) TITLE 
NIT=0 
NPT=0 
NPAR=1 
DA (1) =64.0D0 
WRITE (6,902) TITLE 
DO 3 I=1,100 
READ (5,903) PPW,ETA,TEMP,ALPHA 
IF (PPW(1).LT.0.0D0) GO TO 4 
DEOT=DBLE (ETA/ (TEMP+273.15)) 
DO, 2. J=1,8 
IF (PPW(J).EQ.0.0D0) GO TO 2 
NPT=NPT+1 
DY (NPT) =PPW (J) -ALPHA 
DX (NPT) =DEOT 
DM (NPT) =MVEC(J) 
MWT (NPT) =1 
TEMP V (NPT)=TEMP 
CONTINUE 
CONTINUE 
WRITE (6,904) (I,DA(I) ,I=1,NPAR) 
CALL DERIV(D,DA,DX,DM,DY,4WT) 
DO 6 IB=1,NPAR 
DA(1B)=DA(IB) #D(IB) 
WRITE (6,905) NT 
WRITE (6,906) (I,DA(I) ,I=1,NPAR) 
DDDD=DA (1) **, 3333333: 
WRITE (6,907) DDDD 
IP ((NIT.GT.0) .AND.(NT.LE.NIT)) GO TO 5 
IF (NB.EQ.1) GOTO 5 
WRITE (6,908) 
DO 7 IC=1,NPT 
DPR= DPY (DX (IC) , DM (IC) ,DA, PAR) 
DEL= DY (IC) -DPR 
DS=DS+DEL*DEL 


DOT=DOT+DEL 
7 WRITE (6,909) DX(IC),DY(IC),DPR,DEL,MWT (IC), DM (IC) ,TEMPV (IC) ,TR 


an MewWwh 


= aT ve oi one 00 , Ta4 
whoreeseiXs ZS ao2TaVIN OF 


-Aspoee AASKTA 272eRIS 0% 


‘ae 


(2-8 e (EZAT OD OW Ew. Be 


(Sua rogve Hose ACY 208K SER 


1g, {F085 fily HOE) 


a?. 
" 


‘ 


\Oae.£,002.5,002st 002. 002.000 Fam 


q : = 


I ‘ 


_ 87, (DI) vanat, (21) MOG (IE) TWH USO gH (EP TE: 


i 


a 7 
| nas 


8 


901 
902 
903 
904 
905 
906 
907 
.908 
909 
-910 


236; 


DS=DSOQRT (DS/DPLOAT (NPT-NPAR) ) 
WRITE (6,910) DS,DOT 

WRITE (2,901) TITLE 

WRITE (2,911) DDDD,DS,DOT 
DS=0.0D0 

DOT=0.0D0 

NT=0 

NF=0 

NPT=0 

GO TO 1 f 
WRITE (6,912) 

STOP 0 , 

FORMAT (1044) 

FORMAT (*1',10A4) 

FORMAT (11F5.0) 


PORMAT (* LSQSM PIT : INITIAL GUESSES*/(*ODA(*,I2,*) = *,1PD12.5)) 
PORMAT (*OAFTER *,I2,! ITERATIONS THE PARAMETERS NOW ARE :' 
FORMAT ("ODA(',I2,") = *,1PD18.8) e 

FORMAT (! Ee aicarscns; ',F10.3) 


FORMATS (*=EBOT, Y (EXE) ,Y (THE) ,E-T,WT,A, TEMP, TR*) 

PORHAT  (D71.3, 3F10.2,15, 2F6.1,D12.4) 

FORMAT ('OOBSERVED-CALCULATED STANDARD DEVIATION = *,1PD18.8/ 
*OSUMMED DIFFERENCES = ',D18.8) is 


911 PORMAT (3D20.7) 
912 FORMAT (*ODATA EXHAUSTED ; ENDFILE ENCOUNTERED ON READ‘) 


+ 


END 


REAL FUNCTTON DPY*8(DX,DM,DA,NPAR) 

IMPLICIT REAL*8(A-H,0-Z) 

REAL*8 DA(1) 

COMMON /LWP/ BS,B1T,B2T,GS,G1T,G2T, DS, WO 

COMMON /PARDPY/ F,U,TR,WTR2 

TR=3.03439D-10*DX*DA (1) 

WIR2= (TR*WO) **2 

U=1.0/(1.04+WTR2) 

P=0* WIR2 

DPY= (BS+U* (B1T-P*B2T) + (GStU* (G1T#+F*G2T) +DS*DM* (1.+0* (1. +7) )) 408) * 
DM*TR 


SUBROUTINE DERIV(D,DA,DX,DM,DY,M¥T) 
IMPLICIT REAL*8(A-H,0-2Z) 

REAL*8 DA(1),D(1) ,DX(1),DM(1) ,D¥ (1) 
INTEGER MWT (1) 

COMMON /NAME/ NT,MB,NPT,NPAR 

COMMON /LWP/ BS,B1T,B2T,GS,G1T,G2T, DS, 8O 
COMMON /PARDPY/ F,U,TR,WTR2 

EXTERNAL DPY 

pO 1 IA=1,NPT 

DMM=DM (IA) 

DSDS=DS* DSS 

DXX=DX(IA) 


_ DPR= DPY (DXX,DMM, DA, NPAB) 


A= (BS+ (GS+DSDM) *DMM) 
B= (B1T+ (G1T+DSDS) *DMB) 

C= (- B2T+ (G2T+DSDM) *DMB) 

D (1) =DMM*3. 034 39D-10*DXX* (A*U* (B* (1.-F-F) #C¥F* (3.-4. *Fy)) 
D (2) = (DY (IA) -DPR) 

CALL LSQSM(0,NPAR,D,MWT (IA) ,MA) 

CALL LSQSM(1,NPAR,D, MA, SB) 

NT=NT#1 

RETURN 

END 
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SUBROUTINE. COMPUT 

IMPLICIT REAL*8(A-H,0-Z) 

COMMON /LWP/ BS,B1T,B2T,GS,G1T,G2T,DS,¥O 
READ (5,901) AX,AY,AZ,GX,GY,GZ 

READ (5,901) WO,BO 
STOG=3.9390D-7/ (GX#GY+#GZ) 

AX=AX*1.0D9 

AY=1.0D9*AY 

AZ=AZ*1.0D9 

WO=WO*6. 28319D9 

WRITE (6,902) AX,AY,AZ,GX,GY,GZ,WO,BO 
B=.6666667* (AZ~.5* (AX+AY) ) 

C=.25* (AX-AY) 

DDG= (GZ- .5* (GX +GY) ) *8. 7944D6*BO 

DG=. 5* (GX-GY) *8. 7944D6*BO*C 

BDDG=B*DDG 

DDG2=DDG*DDG 

BS=B*B 

CS=C¥*C 

AOWO= (AX+AY+AZ) /(3. *WO) 

B1T= (.2* BDDG-. 13333333*DDG2*AOWO+. RS pe sag ae 
B2T= (AOWO* (. 1333333*DDG2+5.5125*BS) ) *STOG 


BETA (NS) =U* (B1T-FB2T) 


GS=(.125*BS+. 6 6666667*CS-.2333333*AOWO*BDDG) *STOG 
G1T= (-.025*BS-. 1666667*BDDG*AOWO-. 13333333*CS) *STOG 
G2T= ((. 125*BS-.4*BDDG) *AOWO) *STOG 


GAMMA (NS) =U* (G1T#+P*G27) 


DS=(.05*BS*AOWO) *STOG 
BS= (.2666667*BDDG+1. 066667*DG-AOWO* (. 17777778*DDG2+. 8625*BS) ) *STOG 
RETURN... ; 
901 FORMAT (6F 10.0) 
902 PORMAT (*1AX=',D12.4/% AY=",D12.4/" AZ=",D12.4/8 tag): 3/ 
#* GY=',P12. ar GZ=',F12.3/" WO=",E£12.4/" BO=',F12. 2) = 
END 
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The tables in the Appendix give the peak to peak 
widths of the lines of the ESR ae ee The tempera- 
ture is given in °C and the viscosities are in poise. 
The peak to peak widths were calculated from the peak 
to peak heights of the observed ESR lines and the 


width of the M = -3/2 Lines? 


4 Linewidths which could not be determined due to the 


noise in the spectrum are indicated by a. 


Temperature quoted are accurate to + 0.5°C. 
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